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GLOBAL STRAIN ACCUMULATION AND RELEASE AS REVEALED 
BY GREAT EARTHQUAKES 


By Huco BENroFrr 


ABSTRACT 


The strain-rebound characteristic of the sequence of all great shallow earthquakes of magnitudes 8.0 and 
greater which have occurred since 1904 exhibits a saw-tooth shape with very nearly linear segments. The 
serration amplitudes and periods have been decreasing linearly with time since the beginning of the sequence. 
If a proposed interpretation is correct, the characteristic indicates the following conclusions: (1) Earth- 
quakes in this magnitude range are not independent events, but are related in some form of world-wide 
stress system. (2) From 1908 (and presumably 1904) to 1950 the rate of total secular strain generation in 
the crustal layer in which these earthquakes originate was remarkably constant. (3) The strain was released 
in five active periods of decreasing lengths separated by quiescent intervals of very small or no activity 
during which crustal strain accumulated at a constant rate. (4) During the active periods the strain release 
proceeded at approximately twice the rate of secular strain generation. 

The strain-rebound characteristic of the sequence of all great earthquakes of intermediate depth (h = 
} 70-300 km) exhibits no resemblance to that of the shallow sequence. 

The strain-rebound characteristic of the sequence of all great deep earthquakes (h > 300 km) can be 
represented by the equation S = a + b log ¢ indicating that activity at this depth has been falling off con- 


a tinuously since 1904. 

The different behavior of the three sequences is considered evidence for the existence of three layers in 
A. the earth’s crust having different relative movements in the tectonic sense. 

Heretofore evidence for layering in the earth’s crust has been derived from studies of the propagation 
rg of seismic waves and is thus limited to short-time mechanical properties (less than 5 minutes) only. The 


world strain-rebound characteristics presented herein provide the first evidence for layering on the basis of 
secular tectonic mechanical properties. From this point of view, the crust is made up of three distinct lay- 
ers each with its own proper movement, the shallow 0-70 km, intermediate 70-300 km, and deep 300-680 
km approximately. 














2. Observed periods in strain-rebound charac- 
teristic. Sequence of all great shallow earth- 
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3. Origin times of terminating shocks of active 
intervals in the sequence of great shallow 
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INTRODUCTION 


The elastic strain-rebound accompanying or 
producing an earthquake can be calculated from 
the seismic energy by means of a method 
described in an earlier paper (Benioff, 1949). 
The energy J is derived from the earthquake 
magnitude M in accordance with the following 
relation due to Gutenberg and Richter (1942). 


log J = 12+18M (1) 


Briefly, the method assumes that an earthquake 
is generated by the elastic rebound of a volume 
V of rock having an average elastic constant yu, 
and an average elastic strain e preceding the 
earthquake. The energy stored in the rock as 
elastic strain is therefore 


J = hueV (2) 


Assuming the conversion efficiency of strain 
energy into seismic waves is unity—which 
cannot be far wrong—J is then also the energy 
of the earthquake and « is also the strain re- 
bound. Taking the square root of equation (2) 
we may write 


Ji = (}uV) te (3) 


In general it is not possible to determine yu 
and V but for any given active fault it may be 
assumed that they are constant so that 


Ji=C (4) 


where C = (34V)?. The square root of the earth- 
quake energy is thus proportional to the elastic 
strain-rebound. Consequently, for a sequence of 
shocks occurring on a given fault, a graph of 
the accumulated sum S of the increments 
J;* plotted against time, represents the elastic 
strain-rebound characteristic (times C) of the 
sequence. 

A study of the elastic strain-rebound charac- 
teristics of a number of aftershocks indicated 
that aftershocks are produced by elastic after- 
working of the fault rock (Benioff, 1950). 
Applied to a sequence of earthquakes other 
than aftershocks, the method has provided evi- 
dence that oceanic deeps are the surface ex- 
pressions of great faults which dip under the 
neighboring continents (Benioff, 1949). 
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SEQUENCE OF GREAT SHALLOW EARTHQUAKES 


In this paper the method is applied to three 
world sequences of great earthquakes extending 
back to 1904 when reliable instrumental obser. 
vations first became available. Attention will 
be directed first to the sequence of great shallow 
earthquakes (h <70 km) of magnitudes 8 and 
larger. The data for this sequence are taken 
from the catalog of Gutenberg and Richter 
(1949) and represent a total of 48 shocks!. Since 
the earthquakes of this sequence occurred on 
many different faults having presumably dif- 
ferent values of u and V, it was not anticipated 
that a plot of the strain-rebound characteristic 
would be of significance. However, when the 
graph was constructed (Fig. 1), its singular 
form and regularity indicated that this pessi- 
mistic view was not justified. Presumably u and 
V or the product uV are substantially constant 
for all faults on which these great shocks occur. 
The characteristic approximates a saw-tooth 
curve with amplitudes and periods which de- 
crease linearly with the time. Thus, beginning 
in 1904, the shocks were occurring frequently 
and the magnitudes were so related that the 
strain-rebound increments formed a large, 
nearly linear curve-segment of relatively long 
duration. The activity representing this seg- 
ment of the curve ceased abruptly with the 
earthquake of 1907 October 21. Thereafter, ex- 
cept for the two shocks of 1911 January 3, and 
1912 May 23, all of the faults were relatively 
quiescent for some 94 years until the earthquake 
of 1917 May 1 initiated the next active period. 
In this period, the average strain release con- 
tinued at a nearly constant rate for 7 years and 
ended with the shock of 1924 April 14. Except 
for 3 earthquakes, the faults again remained 
relatively quiescent until the new cycle of ac- 
tivity began 1931 August 10. Similar alternate 
intervals of activity and quiescence, but with 
decreasing durations and amplitudes, have been 
continuing up to the date of this writing. 

The initial points A, B, C, and D of the active 
intervals mark the times when the accumulated 
world strain was maximum. The end points 4, 
b, c, d, and ¢ of the active intervals mark the 


1 Data for the last 3 shocks which occurred since 
the catalogue went to press were kindly given to me 
directly by Dr. Gutenberg. 
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times at which the accumulated world strain 
passed through minima. The upper dashed line 
drawn through the points of minimum strain 
must therefore represent the total secular 
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alternating intervals of activity and quiescence, 
it is necessary to assume the existence of some 
world-wide mechanism which locks all of the 
faults simultaneously, possibly through varia- 
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FicurE 1. STRAIN-REBOUND CHARACTERISTIC. SEQUENCE OF ALL GREAT SHALLOW 
EARTHQUAKES MAGNITUDES > 8.0 
Zach black rectangle represents the elastic strain-rebound increment accompanying a particular earth- 
quake sequence. The height of the rectangle is proportional to the size of the increment and is related to 
the earthquake magnitude as shown in the legend. Blank spaces between rectangles have been provided in 
the few cases where the time-interval between two shocks is so small as to make the two corresponding 


rectangles appear as one. 


strain generation (times C of equation (4)) 
of the system. The small spread of points 
on this line indicates that the rate of strain 
generation has been remarkably constant since 
1907 and presumably since 1904. For any 
given date, the vertical separation between 
corresponding points on the generation line 
and the strain-rebound characteristic is a meas- 
ure of the accumulated world strain (times C) 
available for the production of earthquakes in 
the magnitude range here studied. Since during 
the active intervals the average strain-release 
segment is very nearly linear, it may be con- 
cluded that earthquakes in this magnitude 
range are not independent events, in that the 
release of strain by one shock effects the time 
and presumably the amount of strain-release 
of following shocks. In order to explain the 





tions in the normal stress on the fault surfaces. 
Thus, although global strain is being generated 
at a constant rate given by the slope of the 
upper dotted line, it is released only during the 
intervals when the fault-locking is reduced. 

At the present date, the oscillatory amplitude 
of the characteristic is so small that there is no 
effective secular accumulation, and strain is 
being released substantially as fast as it is being 
generated. Referring again to Figure 1, it may 
be seen that the serration amplitudes of the 
strain-rebound curve decrease linearly with 
the time, since they define the two converging 
dashed lines (Fig. 1). The consecutive periods 
also decrease linearly with time as shown by the 
curve (Fig. 2). 

Figure 3 is a plot of the origin times of the 
terminating earthquakes of the active intervals. 
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Ficure 2. OBSERVED PERIODS IN STRAIN-REBOUND CHARACTERISTIC. SEQUENCE OF ALL GREAT 
SHALLOW EARTHQUAKES 


The time-intervals ab, be, cd, and de (Fig. 1) between accumulated world strain minima are plotted as 


ordinates and the times of the beginnings of the intervals as abscissas. 
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FicurE 3. OriGin Tres oF TERMINATING SHOCKS OF ACTIVE INTERVALS IN THE SEQUENCE OF GREAT 
SHALLOW EARTHQUAKES. 
The shocks are designated a, b, c, d, e, and f to correspond with Figs. 1 and 2. 


It appears that the great earthquake of 1950 
August 15, designated f on the figure, marks 
not only the end of the minor active interval 
with which it is associated, but also the end of 
the whole secular period which began prior to 
1904. 

Epicenters of the terminating earthquakes 
are plotted on a world map in Figure 4. The 
dashed circle drawn to fit the points is very 
nearly a great circle and hence defines a plane 
passing through the center of the earth. The 


pole of this circle, indicated by the two small 
concentric circles on the map, is situated ap- 
proximately at ¢ = 55° N, A = 127° W. It 
may be mere chance that the six epicenters fal 
so nearly on a plane. If not, then this plane 
may define some fundamental shearing surface 
which divides the earth into hemispheres. 

In view of the short interval during which 
observations have been available, it is not pos 
sible to determine whether or not the charactel- 
istic is cyclic in nature. If it is, the amplitude 
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and periods should begin to increase within a 
short time. On the other hand, if the charac- 
teristic represents a transient oscillation, there 
may be a long quiescent interval during which 
strain is accumulated, before the oscillatory 
activity begins again. 
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FicuRE 4, EPICENTERS OF TERMINATING EARTH- 
QUAKES OF ACTIVE INTERVALS IN THE SEQUENCE 
oF GREAT SHALLOW EARTHQUAKES 
The shocks are designated a, b, c, d, e, and f in 
accordance with Figs. 1 and 2. 


A possible cause of the simultaneous locking 
of all the great faults is a variation in the di- 
mension of the earth’s radius such as has been 
postulated by astronomers to explain the 
observed changes in the length of the day. On 
this basis, to account for the shape of the strain- 
rebound characteristic, it would be necessary 
to assume that the earth is executing some form 
of radial relaxation oscillation in which both 
the period and the amplitude have been de- 
creasing since 1904 or earlier. During the phases 
of oscillation when the radius is short, the 
normal stress on the faults is increased and 
they become quiescent. Likewise, when the 
radius is long, the normal stress is reduced and 
the faults become active. 

Another possible interpretation of the origin 
of the sawtooth-shaped characteristic is sug- 
gested by its similarity to the motion of a bowed 
string. Thus the constant speed of the bow cor- 
responds with a constant rate of strain genera- 
tion within the crust. At the conclusion of an 
active period the faults all become locked cor- 





responding with the phase in the bowed-string 
cycle when the string becomes attached to the 
bow. During the quiescent period, the strain 
accumulates at a constant rate in the same way 
that the string is displaced by the bow at a 
constant rate. Finally, when the stress becomes 
greater than the locking friction, the faults 
begin to slip and this corresponds with the 
condition in the string cycle when the bow loses 
hold of the string. The constant rate of strain 
relief during the active periods would thus be 
determined by some form of solid friction in- 
volving the whole shallow crustal layer of the 
earth. The change in period and amplitude 
with time can be attributed to a secular change 
in normal force on the faults rather than a 
short-time variation determined by the saw- 
tooth period as suggested in the preceding 
hypothesis. In the string analogy, a change in 
bow pressure changes the amplitude of string 
vibration but not the period. This results from 
the fact that the string period is determined by 
mass and elastic reactions. On the other hand, 
in the general crustal movements of the earth, 
the mass effects are negligible, and the period 
is determined by elastic and frictional forces 
only. In this type of oscillation it should be pos- 
sible to set up conditions in which period is 
proportional to amplitude. 

By either hypothesis, it is difficult to explain 
the fact that the magnitudes of the shocks 
appear to be independent of the total accumu- 
lated crustal strains. When the strain is large 
it is relieved by a large number of shocks oc- 
curring during an active interval of long dura- 
tion rather than by a small number of larger 
shocks occurring in a short active period. 

If the shallow sequence is extended to include 
all of the world class A shocks of Gutenberg 
and Richter by adding the remaining shocks of 
magnitudes 72 to 8.0, the characteristic loses 
most of its saw-tooth form and, except for the 
1904-1908 segment, becomes approximately 
linear (Fig. 5). It appears therefore that the 
great earthquakes of magnitude 8 and larger 
differ in kind from those of lesser magnitude. 
Apparently the larger earthquakes are produced 
by a primary strain pattern whereas the smaller 
shocks represent some form of auxiliary strain- 


relief. 
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SEQUENCE OF GREAT EARTHQUAKES OF 
INTERMEDIATE DEPTH 


The data for this sequence also taken from 
Gutenberg and Richter (1949, p. 133-139), 
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SEQUENCE OF GREAT DEEP EARTHQUAKES 


The data for this sequence are also taken 


from Gutenberg and Richter (1949, pp. 140- 
141). The sequence includes all earthquakes of 
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FiGurRE 5, STRAIN-REBOUND CHARACTERISTIC. SEQUENCE OF ALL SHALLOW EARTHQUAKES 
MacnituDEs > 7} 


represent all world earthquakes of magnitude > 
7 in the depth range 70-300 kilometers. There 
is some uncertainty as to the eligibility of 
some shocks in this range, since, in some regions 
such as Tonga-Kermadec, shocks in the 70-300 
km depth range belong with the deep shocks 
(Benioff, 1949). Consequently, in this region 
there should be only two classes, shallow, 
0-70 km and deep, 70-680 km. However, until 
a more extended study is made of this distribu- 
tion, it will suffice to use the present classifica- 
tion. The strain-rebound characteristic of this 
sequence (Fig. 6) bears no resemblance to that 
of the shallow sequence (Fig. 1). It suggests a 
long term periodicity, but since less than two 
cycles are represented, no definite conclusions 
can be drawn. 


magnitudes > 7.0 with focal depths of 300 
km and over. The uncertainty as to classifica- 
tion of some of the intermediate shocks men- 
tioned in the preceding paragraph affects this 
sequence also. 

The elastic strain-rebound characteristic of 
this sequence is shown in Figure 7, plotted on 
semi-logarithmic coordinates. This character- 
istic resembles neither the shallow nor the inter- 
mediate sequences. In 1904 when observations 
first began the rate of strain release was maxi- 
mum and it has been decreasing continuously 
since that date. The straight line is calculated 
from the equation shown on the figure and is 
of the form of a compressional elastic (recover 
able) creep-strain. The fact that the observed 
points fit the curve so well suggests that the 
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great deep earthquakes also are related in a 
single stress system and that in this system 
compressional elastic creep strain has been 
accumulating or releasing since 1888. In any 
case the curve shows definitely that the rate of 
elastic strain-rebound in this sequence has been 
decreasing as 1/t since 1904. 


CRUSTAL STRUCTURE OF THE EARTH 


Heretofore knowledge of the interior of the 
earth has been derived principally from the 
behavior of seismic waves. Since the largest 
observed periods in these waves do not exceed 
several minutes, the’ properties of the medium 
exhibited by these observations must be limited 
to short-time characteristics only. A medium 
can transmit shear waves and behave other- 
wise as a solid for waves of such short duration 
and yet be incapable of maintaining a shearing 
stress for as long as several hours or days. Thus 
the properties of the crust defined by waves 
do not necessarily refer to its tectonic or secular 
properties. On the other hand, elastic strain- 
rebound characteristics of earthquake sequences 
do provide information as to the tectonic cha - 
acteristics of the crust. Thus (Benioff, 1949) 
in the region covered by the Tonga-Kermadec 
sequence and the South American sequences, 
there exists a profound tectonic discontinuity 
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at the base of the shallow layer some 60-70 km 
from the surface and; in the region of the South 
American sequences, there is an additional 
secular discontinuity at approximately 300 km 
depth. The fact that the three world-strain- 
rebound characteristics of Figures 1, 4, and § 
are so strikingly different is evidence that two 
at least of the three tectonic discontinuities 
found for South America and Tonga are in gen- 
eral world-wide. Since the tectonic movements 
in these separate layers appear to be unrelated, 
it must be concluded that they actually move 
relative to each other. 
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GEOLOGY OF THE BARRILLA MOUNTAINS, TEXAS 
By G. K. Errter, Jr. 


ABSTRACT 


The Barrilla Mountains, in the northeastern part of the Davis Mountains of Trans-Pecos Texas, are 
composed of Tertiary volcanic materials. Five tuffs and five lava flows, 1500 feet thick occurring 
throughout the mountains, persist in thickness and lithologic characteristics. Their upper surfaces show 
little erosion. The lavas are chiefly Silicic and soda rich. 

The volcanic succession is underlain by a Tertiary sandstone above Upper Cretaceous marine forma- 
tions. These were slightly deformed by the Laramide revolution, subsequently beveled, and everywhere 
covered by the sands of coalescing streams. The sandstone contains well-rounded chert and quartzite 
pebbles. 

Broad folds and normal faults succeeded the extrusions of the youngest lava. 
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tains resemble the more extensive Davis Moun- 
tains, which comprise alternating tuffs and 
lavas, and overlie marine Upper Cretaceous 


August 1947 and June and July 1948 with as- 
sistance of Joseph A. Kennedy and Roger S. 
Plummer. Mrs. Kathryn O. Dickson of the 
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Ficure 1.—INpEx Map SHowrnc Location oF BARRILLA MOUNTAINS 


strata in a structurally low area. Uplifts on the 
north, east, and west expose Lower Creataceous 
and Permian rocks. Southward the Davis 
Mountains merge with well-dissected and al- 
luvium-covered tuffs and flows. No doubt some 
uplift occurred after the lavas were formed. 
The volcanic rocks of the Davis and Barrilla 
mountains are generally considered early Ter- 
tiary, but diagnostic fossils are absent. 

Field work was conducted during July and 


Bureau of Economic Geology, The University 
of Texas, described the igneous rocks petro- 
graphically. 

A grant from the Penrose Bequest of The 
Geological Society of America defrayed expense 
of part of the field work and cost of aerial photo- 
graphs and thin sections. This paper is one of a 
series reporting results of a Trans-Pecos Texas 
project sponsored by the Bureau of Economic 
Geology, The University of Texas. 
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GENERAL STRATIGRAPHY 


The oldest rocks cropping out in the Barrilla 
Mountains are Cretaceous, but wells drilled 
for oil have encountered rocks as old as Permian. 
Older Paleozoic rocks probably occur at greater 
depth. A thick Permian section comprises 
marine limestone, dolomite, subordinate shale, 
and evaporites near the top. The Leonard, 
Guadalupe, and Ochoa series have been iden- 
tified. The Triassic comprises nonmarine red 
and gray shales and sandstones a few hundred 
feet thick. The Jurassic is absent. Approxi- 
mately 1500 feet of marine sandstone, shale, 
marl, and limestone constitutes the Cretaceous. 
Although these rocks range from Trinity to 
Taylor in age, only Washita, Austin, and Taylor 
strata crop out. The Tertiary comprises thick 
tuffs and lavas with minor fresh-water limestone 
and sandstone. Thin gravels and alluvium con- 
stitute the Quaternary. 


PERMIAN SYSTEM 


In the Humble Oil and Refining Company 
No. 1 McCutcheon well east of Star Mountain, 
Skinner (1948, p. 61) identified alternating 
sandstone, limestone, and dolomite as Capitan 
(?) in age, the upper part of which possibly is 
Rustler. 

In the Humble Oil and Refining Company 
No. 2 Flores well northwest of the Barrilla 
Mountains, Creataceous rocks directly overlie 
the Permian, of which the Rustler formation, 
with possibly the Salado and Castile, the 
Capitan limestone, and the Word were identi- 
fied. However, some of the beds assigned to the 
upper part of the Capitan may belong to the 
Tessey limestone. The identification of many of 
the formations is not too certain (Skinner, 
personal communication). 

Northeast of the Barrilla Mountains in the 
Argo Oil Corporation No. 1 Dora Roberts well, 
Reeves County, W. A. Waldschmidt (personal 
communication) identified Cretaceous, possibly 
Triassic, Permian, and Pennsylvanian rocks. 
The Permian section, approximately 11,500 
feet thick, includes well-developed Ochoa, 
Guadalupe, Leonard, and Wolfcamp series. 

From available subsurface data Cretaceous, 
Triassic, Permian, and Pennsylvanian rocks 
underlie the Barrilla Mountains area. 


CRETACEOUS SYSTEM 
Comanche Series 


Georgetown group—The oldest Cretaceous 
strata crop out in the low hills of limestone 2-3 
miles east of the U ranch. On the upthrown side 
of a fault trending northwest, limestone beds, 
white to light-gray, contain abundant rudistids, 
some Toucasia sp., and one specimen of a large 
Neithea. Fossils could not be broken out. These 
beds lithologically resemble the Georgetown 
rudistid facies limestone cropping out along 
the Fort Stockton-Balmorhea highway to the 
north. The lack of identifiable fossils precludes 
a definite correlation, but these beds are tenta- 
tively assigned to the Georgetown group. 


Gulf Series 


Austin group.—-Upper Cretaceous marine 
strata crop out along the periphery of the 
mountains except to the southwest where cov- 
ered by Tertiary volcanics along the axis of 
a syncline (Pl. 1). Good exposures are rare. 
No section justified measurement and descrip- 
tion. 

Outcrops are mostly blue green marls weath- 
ering ochre yellow, but locally a white to pearl- 
gray chalky limestone alternates with poorly 
bedded white indurated marl and marly lime- 
stone. These chalky and marly limestones crop 
out in several ravines about 73 miles west of 
the U ranch. They contain fragments of a large 
Inoceramus, probably Inoceramus undulato- 
plicatus, with small cupped oysters attached. 
Parts of a large rudistid, either Durania sp. or 
Sauvagesia sp., are associated with the Inocera- 
mus. The exposed section here is only about 20 
or 25 feet thick and apparently underlies the 
yellow-weathering Taylor marls. 

In a shallow ravine 6 miles northeast of the 
Jeff ranch a white chalky limestone, inter- 
bedded with cream marl, weathers into small 
flat chips. A few beds, 1-2 feet thick, resemble 
the Austin chalk of central Texas. Gryphaea 
sp. and a large Inoceramus, possibly I. undulato- 
plicatus, occur in these beds, which apparently 
underlie Taylor marl. The contact between the 
chalky limestone and overlying marl is not 
exposed. Chalky limestone crops out about 7 
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miles west of the Jeff ranch and about 3 miles 
southeast. 

The chalky and marly limestones probably 
are Austin chalk. This correlation, however, is 
based upon the very meager evidence of a few 
fossils and uncertain field relations with the 
Taylor marl. 

Taylor group.—Taylor marl is exposed along 
ravines by undercutting and slumping; else- 
where it is covered. It is massive and moderately 
well consolidated. It weathers ochre yellow and 
erodes into rounded smooth hills. About mid- 
way between the U ranch and Jeff ranch, yellow 
to brown sandstone concretions occur at an 
unknown stratigraphic level in the Taylor marl. 
Approximately 1.1 miles southeast of Star 
Mountain 85-90 feet of Taylor marl containing 
abundant Gryphaea newberryi (?) underlies Ter- 
tiary fluviatile sandstone. The marl is dark 
gray with interbedded marly limestone and 
indurated yellow to buff marl. Similar dark, 
flaky chalky marls and thin beds of limestone 
crop out on the north slope of a hill west of the 
Fort Davis-Toyahvale highway 4.74 miles south 
of Toyahvale close to an igneous intrusion. 
These dark marls may be stratigraphically 
higher than the Taylor marl. 

Few species of fossils have been found. 
Gryphaea newberryi (?) and Exogyra ponderosa 
are abundant. East of the U ranch the zone of 
abundant Exogyra ponderosa directly overlies 
the zone of abundant Gryphaea newberryi (?). 
Since no complete section of the Taylor marl 
was found, these zones of abundance may al- 
ternate with each other. 


TERTIARY SYSTEM 
McCutcheon Volcanic Series 


Introduction—Volcanic rocks aggregating 
1500-1700 feet overlie marine Cretaceous rocks. 
A succession of five lavas alternating with five 
tuffs, some of which contain sandstones, brec- 
cias, and fresh-water limestones, is here named 
the McCutcheon volcanic series. (Willis 


McCutcheon ranch 7 miles south of Toyahvale.) 
The base is the lowest sandstone and conglomer- 
ate which here overlie marine Upper Cretaceous. 
The top, the highest lava, is exposed along the 
syncline crossed by the Fort Davis-Toyahvale 
highway north of Star Mountain. 


Apparently one unit succeeded another fairly 
rapidly. Some of the tuffs and lavas may have 
been slightly eroded, but nowhere is there 
evidence of channeling or other marked differ- 
ential erosion. However, a few sandstones and 
conglomerates with lava fragments interbedded 
in the series indicate erosion elsewhere in the 
Davis Mountains area. 

Apparently the lavas spread widely over 
fairly flat surfaces, yet they are silicic and pre- 
sumably were viscous. Perhaps they erupted 
from fissures in the Davis Mountains. Each is 
fairly uniform in thickness and lithologic char- 
acteristics. 

Age and correlation—There is very little 
evidence concerning the exact age of the series. 
The few fossils are not diagnostic. The series 
is younger than the Taylor marl and probably 
is Tertiary. E. W. Berry (1919, p. 4) correlated 
plants from the lower tuffs of the Barrilla 
Mountains with those of the Raton and Denver 
formations (lower Eocene). The tooth of an 
Oligocene rhinoceros has been reported from 
the lower tuffs of the volcanic succession of the 
Davis Mountains about 11 miles west of Bal- 
morhea. Land gastropods and fresh-water algae 
in the lower tuffs of the Barrilla Mountains 
near Barrilla Springs can be dated only as 
Tertiary, possibly early Tertiary. Goldich and 
Seward (1948, p. 15) assign Eocene, Oligocene, 
and possibly Miocene age to the volcanic series 
that overlies eroded Upper Cretaceous rocks 
in the southeastern part of the Davis Moun- 
tains. The tuffs of the Tierra Vieja Mountains, 
about 60 miles southwest of the Barrilla Moun- 
tains, contain vertebrate remains identified by 
Bryan Patterson as Oligocene (cited by Goldich 
and Seward, 1948, p. 18). No attempt is made 
to correlate the tuffs of the Barrilla Mountains 
with those of the southeastern Davis Mountains 
or of the Tierra Vieja Mountains. 

The series is divided into three formations 
(from oldest to youngest): Huelster formation, 
Star Mountain rhyolite, and Seven Springs 
formation. 

Huelster formation —The Huelster formation 
forms the lower part of the McCutcheon vol- 
canic series. The type locality is near the head 
of a deep ravine about 1} miles south of the 
ruins of the Huelster ranch house in the north- 
western part of the Barrilla Mountains. The 
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base is the lowest sandstone and conglomerate 
overlying marine Upper Cretaceous; the top 
is the highest tuff directly underlying the Star 
Mountain rhyolite. The formation, almost en- 
tirely tuff, contains thin layers of sandstone and 
conglomerate, lenses of fresh-water limestone, 
and of trachydoleritic lava. At the type locality 
it is approximately 400 feet thick, but it thins 
eastward. 

The formation crops out widely about the 
Barrilla foothills except on the southwest where 
it dips below the surface in a broad syncline. 
In the interior parts of the mountains, out- 
crops occur in deep ravines eroded along major 
faults. The formation is eroded into low rounded 
hills. It forms no distinct topographic break 
with the underlying Cretaceous, but the over- 
lying lava forms high rugged cliffs. 

The Jeff conglomerate member, a sandstone 
and conglomerate unit constituting the lowest 
beds of the Huelster formation averages 25 
feet thick. The type locality is in a tributary 
ravine on the left bank of Horse Thief Canyon 
about 33 miles upstream from its junction with 
Limpia Creek near the Jeff ranch. 

The Jeff conglomerate member lies with 
slight angularity upon Cretaceous beds, mostly 
Taylor marl, but also Austin chalk. The contact 
is sharp and marks the Mesozoic-Cenozoic 
boundary, but is not well exposed. The member 
is overlain with apparent conformity by vol- 
canic tuff, although a surface of erosion might 
separate the two. 

The sandstone is moderately well sorted, 
medium- to coarse-grained, slightly angular, 
and varies from white to light tan to yellow. 
The conglomerate is composed mainly of well- 
sorted pebbles and cobbles, but boulders occur 
at some localities. Some beds are almost entirely 
pebbles and cobbles with little sandstone matrix 
(Pl. 2, fig. 2), and others are mainly sandstone 
with scattered pebbles and cobbles. 

Bedding is poor, but in places indistinct cross- 
bedding occurs, and at a few exposures the 
sandstone is distinctly laminated. In general 
the basal beds are sandstones varying consid- 
erably in thickness and in places cross-bedded. 
Above are conglomerates with abundant peb- 
bles and cobbles in a sparse, sandy matrix. 
At some localities beds of conglomerate are 
separated by a thin sandstone or by pudding- 


stone. Above the conglomerates are sandstones 
with a few scattered pebbles. 

Most of the pebbles, cobbles, and boulders 
are ellipsoidal to subspherical (Pl. 2, fig. 1), 
but some are greatly flattened. They are chiefly 
white to light lemon-yellow structureless 
quartzite. Others are gray and dark limestone, 
white, black, pink, and gray chert, chert con- 
glomerate, and quartz. Some pebbles are com- 
posed partly of chert and partly of dark lime- 
stone. Fragmentary thick shells of Inoceramus 
(?) occur sparingly. 

At the base of Santiago Peak about 70 miles 
south-southeast of Star Mountain the writer 
(Eifler, 1943, p. 1635) found a similar sandstone 
and conglomerate overlying eroded Boquillas 
formation (Upper Cretaceous) and underlying 
Tertiary volcanic ash. Most of the pebbles in 
the conglomerate are chert, but some are lime- 
stone. South of Alpine Goldich and Seward 
(1948, p. 13) report a similar conglomerate be- 
tween the Boquillas formation and the Buck 
Hill volcanic series (Tertiary). The conglomer- 
ate is composed “. . . principally of well-rounded 
pebbles and cobbles of limestone, with some 
pebbles of marble, quartzite, chert, and chalce- 
dony in a well-cemented matrix of tuffaceous 
sandstone” (Goldich and Seward, 1948, p. 18). 
The conglomerates south of Alpine and at 
Santiago Peak are probably the same and may 
be equivalents of the Jeff conglomerate. 

The source of the pebbles and cobbles is 
certainly not the underlying Cretaceous forma- 
tions. However, fragmentary shells of Jnocera- 
mus were probably locally derived; they show 
practically no abrasion. The quartzite pebbles 
increase in size toward the Van Horn uplift, 
but no extensive quartzite is known there. 
The chert, chert conglomerate, and dark lime- 
stone resemble some Paleozoic rocks of the 
Marathon basin. In the northeastern Glass 
Mountains J. T. Lonsdale (personal communi- 
cation) found well-rounded quartzite pebbles 
in a Leonard conglomerate, which might be the 
source of the Jeff conglomerate. 

The Jeff conglomerate member is thought to 
be continuous throughout the Barrilla Moun- 
tains and the northeastern part of the Davis 
Mountains. It was deposited on a fairly smooth 
erosion surface either by coalescing streams or 
by extensive sheetwash. It is too extensive to 
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be the veneer of a pediment. The thickness is 
uniform over too great an area for the deposit 
to be a piedmont alluvial fan, although the 
sandstone and gravel might have been spread 
very thin in front of a piedmont fan. 

Sandstone, not typical of the Jeff conglomer- 
ate member, occurs between Cretaceous rocks 
and the volcanics at two localities. In a ravine 
44 miles S. 30° E. of the U ranch, sandstone 
overlies Taylor marl. The basal 6 feet of sand- 
stone, which contains no pebbles, is ochre 
yellow to brown and is cross-bedded and lam- 
inated. The upper contact is not exposed. About 
7 miles west of the Jeff ranch, 20 feet of yellow 
to white sandstone with no pebbles overlies the 
Austin chalk. The sandstone is cross-bedded and 
laminated. On the slope above are a few scat- 
tered pebbles and tuff beds, which may not be 
in place. The absence of pebbles and the yellow 
color of the basal sandstone are thought to be 
local characters, and the two sections are cor- 
related with the Jeff conglomerate member. 

At the type locality the upper contact of the 
Huelster formation is obscured. Highly gen- 
eralized section of the lower 294 feet which is 
exposed is given below. 


PARTIAL SECTION OF HUELSTER FORMATION SOUTH 
oF HvEtsteR House RuINs 


Thickness 
in feet 


Unit 
6. White to cream, fine to coarse tuffs. Petrified 
tree trunks 2-3 feet across. Upper part ends 
at top of hill covered by lava debris. ..... 
. Flinty-appearing cream tuffs. A few ane 
tree trunks. . 6 
Ledges 3-4 feet thick of hard, white | to cream 
tuffs. Brown stains along sharp joints. . 44 
. White, buff, and cream tuffs. . ree 63 
2. Cream to white tuffs. Yellow-brown stains 
along joints. Cross-bedded sandstone near 
top.. 7 ~ ae 


Description 


104 


on 


> 


we 


1. Conglomerate and sandstone. ‘Quartzite ‘and 
chert pebbles near base; gray marls at top. 
Underlain by Taylor marl................ 52 
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Approximately 100 feet of poorly exposed 
tuff with an interbedded lens of trachydolerite 
lies above unit 6 and below the Star Mountain 
rhyolite. 

The Huelster formation varies more in thick- 
ness than any other dominantly tuff section of 
the McCutcheon volcanic series. At a point 
3.5 miles north of the Jeff ranch the formation 
is 136 feet thick, and on the hill 3 miles south- 
east of the Jeff ranch it is only a few feet thick 
and locally is missing. 

The tuffs are distinctly bedded in layers a 
few inches to 1 or 2 feet thick. The flatness and 
continuity of bedding may indicate lacustrine 
deposition. 

At two observed localities, fresh-water lime- 
stone is interbedded with tuff near the base of 
the member. About 1.8 miles south-southeast 
of the Jeff ranch an unfossiliferous limestone 
layer 3-4 feet thick lies 50 feet above the Jeff 
conglomerate member. At the other locality, 
1.6 miles north of the Jeff ranch, a 5-foot lime- 
stone layer crops out conspicuously along the 
foot of several low tuff hills. The limestone, 
light lemon-yellow to light and dark gray and 
compact, contains many irregular purplish- 
biue and ochre yellow cherts of secondary origin. 

Petrified tree trunks and unidentified low- 
and high-spired gastropods occur sparingly in 
the exposed limestone. Concentric, fusiform to 
ellipsoid organic structures, 4 or 5 inches in 
longest diameter, compose a large part of the 
rock. They have been identified as algae by 
J. Harlan Johnson, whose description follows 
(personal communication): 


Megascopic: 
Rounded pellets, with irregularly banded, more or 
less concentric structure, which make up a large 
percentage of the limestone sample. 
Microscopic: 
An algal limestone formed by colonies of an un- 
described type of algae, probably belonging among 
the green algae (Chlorophyta), possibly among 
the blue-green (Cyanophyta). Shows molds of 
coarse thread-like filaments which form rounded 





PiatE 2.—JEFF CONGLOMERATE SOUTHWEST OF JEFF RANCH 


Ficure 1.—ROUNDED QUARTZITE AND CHERT PEBBLES IN SANDSTONE MATRIX 
FicurE 2.—MaAssIvE CONGLOMERATE 
Hammer gives scale. 
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Ficure 2 


JEFF CONGLOMERATE SOUTHWEST OF JEFF RANCH 
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and rosette-shaped colonies (bands and fans in 
section). Average size of filaments 0.005 mm. 


Remarks: 
Probably a fresh-water growth. 


The algae closely resemble specimens which 
S. S. Goldich obtained from a fresh-water lime- 
stone near the base of the Pruett formation 
south of Alpine. Goldich and Seward (1948, 
p. 15) also report an Eocene gastropod, Gonio- 
basis tenera carterii. These limestone beds, al- 
though not continuous, are probably of the 
same age. 

Thick lenses of a dark igneous rock, probably 
lava, occur in the upper part of the formation 
in the western part of the Barrilla Mountains 
and at Wild Rose Pass in the Davis Mountains. 
The rock weathers into smaJl brownish-yellow 
to buff chunks. Megascopically it is a black, 
fine-grained melaphyre; the microscope reveals 
a porphyritic, slightly ophitic texture. The 
rock can be called a porphyritic olivine trachy- 
dolerite. 

Star Mountain rhyolite—The thick rhyolitic 
lavas overlying the Huelster formation are 
named the Star Mountain rhyolite after Star 
Mountain. The base is the lowest lava above 
the Huelster formation and the top is the high- 
est lava below the tuffs of the Seven Springs 
formation. The number of individual flows 
cannot be ascertained. In the high vertical 
cliff on the east side of Star Mountain there ap- 
pears to be at least six distinct flows, but possi- 
bly there are more. At the type locality the 
formation is 496 feet thick, but it thins to the 
northeast. About 3.5 miles north of the Jeff 
ranch it is 262 feet thick. The rock is typically 
black to brownish black, but some parts, par- 
ticularly the lower, are green, greenish gray, 
and purplish lavender. 

Owing to its great thickness and resistance 
to erosion the Star Mountain rhyolite crops out 
widely throughout the area. The lava is exposed 
in vertical cliffs east of Star Mountain, along 
parts of Limpia Creek, and near the summit 
and northeastern part of the Barrilla Moun- 
tains. The formation is commonly surmounted 
by a distinct bench where the overlying tuff 
has been eroded. 

The flows commonly exhibit Palisade struc- 
ture, and columnar jointing is less commonly 
developed. Many joints, which essentially paral- 
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lel the Palisade structure, extend through sev- 
eral flows and must therefore be of tectonic 
origin. Flow structure is commonly developed. 
Some flows are highly vesicular at the top and 
poorly vesicular at the bottom. Blue to gray 
chalcedony fills many of the vesicles. 

The rock is porphyritic with conspicuous 
phenocrysts in an aphanitic groundmass and 
is a riebeckite soda rhyolite porphyry. 

Seven Springs formation.—The Seven Springs 
formation includes the lavas, tuffs, and inter- 
bedded sediments between the Star Mountain 
rhyolite and the top of the McCutcheon vol- 
canic series. The formation is named for Seven 
Springs near the Willis McCutcheon ranch 
house 7 miles south of Toyahvale. A composite 
section near the type locality is 682 feet thick. 
The formation is completely exposed in the 
syncline north of Star Mountain. Elsewhere the 
upper beds are eroded. Where horizontal or with 
low dip the formation erodes into distinct 
benches. 

The formation is divided into four lava mem- 
bers and four tuff members, which alternate 
with one another. 

Tuff member No. 1, at the base of the forma- 
ion, includes the tuffs and interbedded sedi- 
ments between the top of the Star Mountain 
rhyolite and the base of the overlying lava 
member No. 1. It is 102 feet thick on the lower 
northeast slopes of Star Mountain but is con- 
siderably thinner in the Barrilla Mountains. 

The tuff is evenly bedded in layers up to 2 
feet thick. The beds are white, gray, blue green, 
and lemon yellow. At the west end of Kunz 
Mesa, tuff is interbedded with a well-cemented 
conglomerate which contains igneous boulders 
up to 2 feet in diameter. Similar conglomerates 
and breccias with igneous fragments are wide- 
spread but do not occur throughout the Barrilla 
Mountains. Associated with them, but not 
present everywhere, is a blue green, highly 
vesicular lava. 

Lava member No. 1, which includes all 
lavas above tuff member No. 1 and below tuff 
member No. 2, consists of a thin, inconspicuous 
vitrophyre below and a thicker, prominent 
granophyre above, each representing a single 
flow. 

The average thickness of the vitrophyre is 
4-8 feet, but it is at least 30-40 feet thick about 
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1 mile southeast of Seven Springs. Although 
very thin, the vitrophyre must have spread 
over the entire region. 

Typically the rock is black, grading in places 
into light greenish yellow phases. It is notably 
devoid of vesicles and columnar jointing. Large 
phenocrysts are included in a glassy ground- 
mass with perlitic structure. Some parts may be 
fragmental. It is classed as a soda rhyolitic 
vitrophyre. 

The granophyre is persistently 30-40 feet 
thick with excellent columnar jointing and 
Palisade structure. The rock is purplish red, 
maroon, and reddish brown. At the top and 
bottom the flow is distinctly vesicular. The 
base is a conglomeratelike mass of pieces ap- 
parently derived from the granophyre while 
flowing. The rock can be called a spherulitic 
riebeckite soda granophyre. 

Tuff member No. 2 includes all tuffs and sedi- 
ments between the top of lava member No. 1 
and the base of the overlying lava member No. 
2. At Star Mountain the member is 105 feet 
thick, on top of the Barrilla Mountains east of 
the McCutcheon ranch it is 114 feet thick, on 
Kunz Mesa it is 90 feet thick, but in the eastern 
and northeastern part of the Barrilla Mountains 
the interval is much thinner. The tuffs are 
commonly white to gray and are well bedded 
in 1- to 2-foot layers. In places the tuff contains 
conglomerate and cross-bedded sandstone. 

Lava member No. 2 includes all lavas above 
tuff member No. 2 and below tuff member No. 
3. It ranges in thickness from 10 feet to 77 feet 
and forms no distinct topographic bench. 

The base is a red lava (possibly a well-con- 
solidated tuff) up to 2 feet thick, overlain at one 
locality by a 20-foot interval of tuff, but else- 
where by lava. The lava consists of three 
basaltic phases: a chocolate-brown with flow 
layers and blue chalcedonic amygdules, a blue- 
black, scoriaceous, conglomeratelike rock with 
chalcedonic amygdules; and a compact and 
specular phase. These phases are similar litho- 
logically and can be called glassy basalt or tra- 
chybasalt porphyry. 

Tuff member No. 3 includes all tuffs and 


sediments above lava member No. 2 and below 
the overlying lava member No. 3. It is 138 
feet thick at Star Mountain, 227 feet on the 
mesas on top of the Barrilla Mountains east of 


the McCutcheon ranch, 185 feet on Kunz Mesa, 
but it thins rapidly to the east and on Beard 
Mountain is only 25-30 feet thick. 

The interval is mainly white to gray tuff 
and contains very little sediments. 

Lava member No. 3 includes all lavas above 
tuff member No. 3 and below the tuff member 
No. 4. It forms the caprock of Star Mountain, 
of the two mesas on top of the Barrilla Moun- 
tains east of the McCutcheon ranch, and of 
Kunz Mesa. It is 42-65 feet thick. 

The lavas have moderately well defined pali- 
sade structure and abundant platy flow layers. 
The rock is prevailingly gray, but also red and 
black. The rock can be classified as granophyric 
rhyolite porphyry. 

Tuff member No. 4 includes all tuffs, sedi- 
ments, and interbedded lava above lava mem- 
ber No. 3 and below lava member No. 4. It is 
approximately 200 feet thick in the syncline 
between Star Mountain and the Barrilla Moun- 
tains. A single lava flow near the top is 54 feet 
thick. This lava, a vitric rhyolite porphyry, is 
brown at the base and lavender to light purple 
at the top. 

Lava member No. 4 includes all lavas above 
tuff member No. 4 and the top of the McCutch- 
eon volcanic series. At the type locality of 
the series a thin tuff may overlie this lava mem- 
ber, but if so it is covered by alluvium. The 
lava crops out continuously along the flatirons 
north of Star Mountain. A few scattered out- 
crops occur in the extensive alluvial flat north- 
northeast of the junction of the Jeff ranch road 
with the Fort Davis-Balmorhea highway. The 
lava member is at least 41 feet thick. The rock 
is gray and is stratiform in intervals 4—6 inches 
thick; it is a granophyric rhyolite porphyry. 


POSSIBLE INTRUSIONS 


Igneous rocks of unknown relations crop out 
at two localities along the Fort Davis-Balmor- 
hea highway, 4 and 5 miles south of Toyahvale. 
The northern occurrence, apparently a tabular 
body intrusive into the Taylor marl, is prob- 
ably a dike. The southern body crops out more 
widely, underlying two hills west of the high- 
way. Contact with the surrounding Taylor 
marl is not exposed, and it is impossible to tell 
whether the rock is intrusive or extrusive. Some 
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POSSIBLE INTRUSIONS 


marl close to the igneous rock is darkened, 
perhaps by heat. This igneous body might be a 
sill whose roof has been removed by erosion, 
but more likely it is a dike. The rocks of the 
two occurrences are soda trachyte porphyries. 


STRUCTURAL FEATURES 
Introduction 


There is only slight angularity between the 
Cretaceous rocks and the overlying McCutch- 
eon volcanic series. The structural features 
of the volcanic series are probably quite similar 
to those of the Cretaceous rocks, but Cretaceous 
formations were slightly warped and faulted 
before the Jeff conglomerate was deposited. 
Because Paleozoic rocks do not crop out in the 
area and subsurface data are scarce, no at- 
tempt is made to interpret their structural rela- 
tions with the Cretaceous or the volcanic rocks. 


Folds 


Broad folds trend northwestward in the 
region, the most conspicuous of which is a 
syncline separating the Barrilla Mountains from 
the Davis Mountains. The southwest limb with 
a dip of 18°-20° has a structural relief of ap- 
proximately 2500 feet and the northeast limb 
about 2000 feet. A few miles northwest of Seven 
Springs the syncline is mostly covered by the 
alluvium of Toyah Creek. About 2 miles north- 
northeast of the Jeff ranch road junction with 
the Fort Davis-Balmorhea highway a structural 
basin is superimposed on the northeast limb. 
The areal extent and structural relief of the 
basin is not determined because of alluvium. 
Near the Jeff ranch the northeast limb is almost 
horizontal, but the southwest limb dips steeply. 
Below the alluvium of Limpia Creek there is 
either the continuation of the syncline greatly 
narrowed or a fault upthrown to the northeast. 

An anticline trends northwestward through 
Star Mountain, from which the volcanic beds 
dip steeply northeastward into the syncline, 
and southwestward at angles of 1°-2°. This 
structural feature is probably a continuation 
of the Hovey anticline (King, 1937, p. 141). 
Another anticlinal uplift trends northwestward 
through the crest of the Barrilla Mountains. 
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The greatest structural height along this uplift 
is in the western part of the mountains. 

Folding occurred after the outpouring of lava 
member No. 4 of the Seven Springs formation 
and is not related to the Laramide movements. 
The exact date is not yet determined. 


Faults 


Longitudinal high-angle faults cut the rocks 
of the McCutcheon volcanic series into long 
narrow blocks. These faults are fairly abundant 
in the Barrilla Mountains but do not occur in 
the syncline or in the adjoining part of the 
Davis Mountains. The crest of the Barrilla 
Mountains east of the McCutcheon ranch is 
a grabenlike block (Pl. 1). The fault bounding 
the block on the southwest dropped tuff mem- 
ber No. 3 of the Seven Springs formation against 
the Star Mountain rhyolite and has a throw of 
approximately 490 feet. The fault on the north- 
east has less throw but extends for two-thirds 
the length of the Barrilla Mountains. It is a 
rotary fault, and the rocks forming the graben 
become upthrown as they extend southeast- 
ward. The northeast slopes of the Barrilla 
Mountains may contain many small faults 
downthrown to the northeast, but the thick 
mantle of igneous debris prevents an interpreta- 
tion of the structural relations. The lava-topped 
mesas at the northeast margin of the Barrilla 
Mountains are bounded on the southwest by 
faults downthrown to the northeast. This line 
of faults continues southeastward between 
Beard Mountain and the main body of the 
Barrilla Mountains. Narrow fault blocks on the 
northeastern slope of Beard Mountain have 
apparently resulted from landslide faulting. 
They are downthrown to the northeast and are 
tilted at steep angles to the southwest. A fault 
2 miles east of the U ranch separates the George- 
town and Huelster formations. The upthrown 
side of this fault is probably part of an uplift 
paralleling the anticline along the crest of the 
Barrilla Mountains. 

The structural relations are not clear in the 
foothills of the Barrilla Mountains north and 
east of the Kunz ranch. Here the Star Mountain 
rhyolite is apparently faulted against Taylor 
marl. The fault, if present, trends northwest- 
ward at the southwest foot of Kunz Mesa and 
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turns at 90° down the ravine for a short dis- 
tance and then, making another 90-degree 
turn, resumes its northwestward trend. It is 
possible that the Huelster formation here is 
not faulted out but that the Star Mountain 
rhyolite rests directly on Taylor marl. 

The faults displace beds as young as the 
tuffs and lavas of the Seven Springs formation, 
but whether the movement preceded or fol- 
lowed the folding is not clear. Faulting prob- 
ably followed folding. If the faults are genet- 
ically related to folding, they would likely 
have formed along the anticlinal uplift through 
Star Mountain and be as abundant in the north- 
eastern part of the Davis Mountains as they 
are in the Barrilla Mountains. Furthermore, if 
the faults formed in response to the anti- 
clinal uplifts, then the fault block 3 miles east 
of the McCutcheon ranch would likely have 
been a horst instead of a graben. 


GEOMORPHOLOGY 


The Barrilla Mountains have not been ex- 
tensively eroded. The highest parts essentially 
coincide with the greatest uplift. A cross sec- 
tion (Pl. 1) northeast from Star Mountain 
shows that topographic relief is not much less 
than structural relief. Nevertheless erosion has 
been sufficient to form several low fault-line 
scarps where the Star Mountain rhyolite was 
faulted against the less resistant Taylor marl 
and Huelster formation. This relation is well 
shown 2 miles north of the Jeff ranch and on 
the southwest side of the mesas which border 
the mountains on the northeast. The lavas are 
the most resistant rocks; they top the mesas 
and form the flat part of benches and scarps. 
The Star Mountain rhyolite, owing mainly to 
its great thickness, forms high vertical cliffs 
where undermined by the erosion of the 
Huelster formation and the Taylor marl. 


Economic GEOLOGY 


Subsurface water is the most important nat- 
ural resource in the Barrilla Mountains. Small 
gravity springs are fairly abundant in the 
western part of the mountains. Practically all 
of them occur at the contact between the per- 
meable Jeff conglomerate and the impermeable 


Taylor marl. Where sufficiently large the springs 
are dammed, and the water is piped and used 
for watering stock, but each of these springs 
flows very little water and is fairly sensitive 
to fluctuations in rainfall. 

The vitrophyre at the base of lava member 
No. 1 of the Seven Springs formation was 
found unsuitable as a lightweight aggregate 
because of the phenocrysts. Further explora- 
tion, however, might reveal a phase of the 
vitrophyre without phenocrysts. 

Oil has not been produced commercially in 
the Barrilla Mountains, but the thick sequence 
of Permian rocks offers excellent opportunity for 
oil exploration. If structural relations between 
these rocks and the Tertiary lava flows can be 
established, perhaps the structural features of 
the lavas will aid in the discovery of structural 
oil traps in the Permian rocks. 


PETROGRAPHY 
By Katuryn O. DICKSON 
General Statement 


The igneous rocks of the Barrilla Mountains 
belong to the same Tertiary alkalic province 
as those of the Terlingua-Solitario, Big Bend, 
and southern Davis Mountains regions of Brew- 
ster and Presidio counties, Texas (Lonsdale, 
1940; Goldich and Elms, 1949). Anorthoclase 
is the typical alkali feldspar. Plagioclase is 
negligible in the felsic members, and analcime 
is sometimes present in the mafic members. 
In the main, the lava flows are porphyritic 
quartz-poor soda rhyolites which grade to 
quartz trachytes. These lavas are characterized 
by soda pyriboles and granophyric, trachytic, 
or glassy groundmasses. The mafic types are 
more limited in occurrence. Aside from one ex- 
tensive basaltic flow, they appear as local lenses 
which have diabasic texture and are composed 
principally of titaniferous augite, olivine, labra- 
dorite, and magnetite. The important rock 
types of the area are discussed, and rocks from 
specific localities are shown in Table 1. The 
percentages represented are estimates because 
of the highly altered state of most of the mafic 
minerals and the fine-grained groundmasses. 
Where volume percentages were practicable 
and informative, Rosiwal analyses were made. 
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Probably the major portion of the ground- 
masses is anorthoclase. However, as the data 
are not critical, all alkali feldspar is included 
in the table under one heading with a notation 
when anorthoclase phenocrysts are present. 


Riebeckite Soda Rhyolite Porphyry 


Lava flows of the Star Mountain rhyolite with 
very few exceptions are composed of this rock. 
It is fine-grained light gray to reddish with 
abundant feldspar phenocrysts 1-5 mm. in 
size and small grains of dark minerals visible 
in hand specimen. The groundmass is trachytic 
to micrographic with alkali feldspar laths 0.01— 
0.3 mm. in length, averaging 0.04 mm. Quartz 
occurs interstitially or in intergrowth. A basal 
flow (Table 1, No. 6) contains no quartz and 
isa trachyte. Occasionally a trace of riebeckite 
or other soda amphibole is present, but most of 
the original soda amphibole sponges have been 
altered to an opaque black mineral and hema- 
tite. There are discrete grains of iddingsite (?), 
a trace of fluorite, and a few prisms of apatite 
as accessories. Chalcedony and calcite appear 
in vesicles with hematite. 


Spherulitic Riebeckite Soda Granophyre 


This rhyolite makes the top of lava member 
No. 1 of the Seven Springs formation. Speci- 
mens are dull gray to brick red with distinct 
flow lines and a few minute vesicles. Some speci- 
mens have a fragmental appearance both in 
hand specimen and under magnification. Rec- 
ognizable megascopically are anorthoclase phe- 
nocrysts 3-8 mm. in length oriented along flow 
lines. A few phenocrysts are feldspar and quartz 
intergrown micropegmatitically. The ground- 
mass appears aphanitic in hand specimen but 
in several slides consists largely of spherulitic 
intergrowths of quartz and alkali feldspar aver- 
aging 0.2-0.3 mm. in diameter. Rarely riebeck- 
ite is found in needles in the spherulites. Rare 
fine aegirite grains and a few grains of iddingsite 
about 0.5 mm. in diameter are present in the 
groundmass. Opal and chalcedony are abun- 
dant in amygdules together with hematite 
and limonite. In the reddish fragmental type, 
there are spherulites but to a rather limited 
extent. The cryptocrystalline groundmass en- 
closes euhedral feldspar phenocrysts, some of 
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which are much embayed. Patches and stringers 
of hematite are common. The spherulitic frag- 
ments as well as those of micrographic or 
trachytic texture apparently are inclusions in- 
corporated in the rock during flow. 


Granophyric Rhyolite Porphyry 


The lava members No. 3 and No. 4 of the 
Seven Springs formation and one bed in the 
base of the Star Mountain rhyolite section fall 
under this heading. The rhyolite is purple- 
to red gray, with a spotty appearance because 
of numerous light-colored feldspar phenocrysts 
and amygdules. Some flows or facies of flows 
are highly vesicular; others are dense and glassy. 
The groundmass ranges from granophyric 
quartz and alkali feldspar containing altered 
poikilitic ferromagnesian sponges to glass 
densely packed with felty (—0.05 mm.) feldspar 
laths and opaque ferromagnesian dust. The 
anorthoclase phenocrysts are stubby in habit, 
averaging about 4 mm. in length. The mafic 
phenocrysts, as well as the mafics in the ground- 
mass, are altogether altered, but the pheno- 
crysts have the general habit of a pyroxene and 
average 1 mm. in greatest dimension. One 
slide has in addition microlites of a green ferro- 
magnesian, probably aegirite-augite, and a trace 
of biotite bordering the discrete ore grains. 
Hematite and limonite are present as altera- 
tion products, and zircon, apatite, and fluorite 
were found in trace quantities. Opal, chalce- 
dony, calcite, and secondary quartz are com- 
mon amygdules with cristobalite spherulites 
prominent in the vesicular rocks. 


Soda Rhyolitic Vitrophyre 


This glass occurs at the base of lava mem- 
ber No. 1 of the Seven Springs formation; it 
is a black and gray banded glass containing 
shining feldspar phenocrysts, 1-5 mm. in length, 
and fragments of basalt and rhyolite. Some 
samples are easily crumbled due to perlitic 
structure in the glass. The resinous black glass 
is pale brown in section; n = 1.517, +0.002, 
which indicates an SiO, content of about 65 
per cent. The glass has a lower silica content 
than the phenocrysts indicate, these latter 
being anorthoclase, quartz, graphically inter- 
grown quartz and feldspar, augite, and aegirite- 
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* Includes anorthoclase phenocrysts. 
Lenses in the upper Huelster formation: 

Pass. 


1. Porphyritic olivine trachydolerite in tuffs of the Huelster formation, road cut in Wild Rose 

2. Porphyritic olivine trachydolerite in tuffs of the Huelster formation, Wild Rose Pass. 

3. Olivine trachydolerite in tuffs of the Huelster formation, western Barrilla Mountains. 

4. Porphyritic olivine trachydolerite in tuffs of the Huelster formation, western Barrilla Mountains. 
Star Mountain rhyolite: 
5. Riebeckite soda rhyolite porphyry from Star Mountain. 

6. Riebeckite soda trachyte porphyry, lower contact from basal lava in road cut between Fort Davis 
and Toyahvale north of Seven Springs. 
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7. Hornblende (?) soda rhyolite porphyry, gray-green flow. 
8. Soda rhyolite porphyry, Hill P-1, base of Kunz Mesa. 
9. Hornblende (?) soda rhyolite porphyry, southwest side of Kunz Mesa. 
10. Soda rhyolite porphyry, base of formation on CCC road south of Toyahvale. 
11. Soda rhyolite porphyry, igneous hill 4 mile west of Kunz ranch. 
12. Hornblende soda rhyolite porphyry, outlier west of U ranch. 
13. Granophyric soda rhyolite porphyry, basal lava of Star Mountain rhyolite at contact with lower 
ash southwest of Jeff ranch. 
14. Olivine trachydolerite, basal lava directly overlying “flinty” lava above Huelster tuffs south of 
Seven Springs. 
Seven Springs formation: 
15. Spherulitic riebeckite soda granophyre, Kunz Mesa. 
16. Spherulitic riebeckite soda granophyre, Kunz Mesa. 
17. Spherulitic riebeckite soda granophyre, Kunz Mesa. 
Lava member No. 2: 
18. Basalt porphyry, north slope of Star Mountain. 
19. Basalt “conglomerate” lava, north slope of Star Mountain. 
20. Basalt ‘‘conglomerate” lava, Beard Mountain. 
21. Porphyritic olivine trachybasalt, Kunz Meas. 
22. Olivine trachybasalt porphyry, Kunz Mesa. 
23. Olivine trachybasalt porphyry, north of Fort Davis-Toyahvale highway at entrance to McCutcheon 
ranch, 
24. Olivine trachydolerite, north of Fort Davis-Toyahvale highway at entrance to McCutcheon ranch, 
Lava member No. 3: 
25. Granophyric soda rhyolite porphyry, east end of Barrilla Mountains. 
26. Granophyric soda rhyolite porphyry, top of Kunz Mesa. 
27. Soda rhyolite porphyry, road cut at Seven Springs. 
28. Soda rhyolite porphyry, road cut at Seven Springs. 
Lava member No. 4: 
29. Granophyric soda rhyolite porphyry, north slope of Star Mountain. 
30. Granophyric soda rhyolite porphyry along Limpia Creek, isolated hill just east of Jeff ranch turnoff 
from Fort Davis-Toyahvale highway. 
31. Granophyric soda rhyolite porphyry, near McCutcheon ranch entrance. 
32. Granophyric soda rhyolite porphyry, highway cut near McCutcheon ranch entrance. 
Intrusives: 
33. Soda trachyte porphyry, #5 mile toward Toyahvale from fork of old and new Fort Davis-Toyah- 
vale highways. 
34. Soda trachyte porphyry; same as No. 33. 
35. Quartz soda trachyte porphyry, dike (?) in Cretaceous on Fort Davis-Toyahvale highway, south of 
locality for No. 33. 
36. Aegirite-augite soda trachyte porphyry; same as No. 33. 
37. Aegirite-augite soda trachyte porphyry; same as No. 33. 


augite. Fragments of basalt, granophyric rie- Soda Trachyte Porphyry 
beckite soda rhyolite, and plagioclase are also 
present. The glass is devitrified in places border- 
ing amygdules and cracks, and some otherwise 
euhedral crystals are deeply embayed. One 


There are two intrusives of trachyte. The 
rock is a light-gray quartz soda trachyte 
porphyry, with small (1-3 mm.) feldspar pheno- 
specimen contains large amygdules of chalce- Tysts and narrow plates of reddish biotite and 
dony 1-5 cm. long. Under the microscope this hematite. Brown concentric rings 4-% inch 
type shows a groundmass of spherulitic glass across result from weathering. The phenocrysts 
and numerous tiny amygdules (averaging 0.1 are anorthoclase. The groundmass is trachytic 
mm. greatest dimension) of chalcedony. No with alkali feldspar laths 0.04-0.2 mm. in 
perlitic cracks are present. length and a small amount of interstitial quartz 
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varying from a trace to 2 or 3 per cent. A sodic 
amphibole, riebeckite (?), occurs in the ground- 
mass in sponges and shreds and is largely 
altered to an opaque black mineral. Some sam- 
ples contain minute grains of a sodic pyroxene, 
probably aegirite-augite. Magnetite is present 
in discrete grains. Iddingsite (?), hematite, 
and reddish-brown biotite (?) are alteration 
products after mafics, and apatite is an ac- 
cessory mineral. A facies of aegirite-augite 
soda trachyte porphyry occurs in the intrusives. 
This type is distinctive under the microscope 
but megascopically is similar to the weathered 
specimens of the main rock. This peculiarity 
makes the relations of the two kinds somewhat 
equivocal. In hand specimen it is a dark gray- 
green to black fragmental-appearing rock with 
darker-green fragments up to an inch in di- 
ameter seemingly imbedded in the lighter- 
colored matrix. In some specimens the darker 
material is apparently in bands which are 
separated from the light matrix by sharp bound- 
aries. The groundmasses are glass filled with 
alkali feldspar laths and ferromagnesian micro- 
lites, and a greater concentration of the latter 
imparts the darker green color. The feldspar 
laths of the matrix are 0.01-0.2 mm. in length. 
The anorthoclase phenocrysts average 1 mm. 
in length. Veinlets of chalcedony 0.04-0.2 mm. 
wide cut through the rock, and augite and 
aegirite-augite occur in fine grains in the ground- 
mass and in phenocrysts up to 1 mm. in great- 
est dimension. A yellow mineral has completely 
replaced some of the smaller grains and has 
formed an alteration rim around some of the 
larger pieces of the pyroxene. This mineral is 
pleochroic in yellow, 2V about 70°-80° (—), 
ng = 1.68, ny — ne = 0.025, witha prominent 
cleavage, length slow, extinction about 7°-9° 
from the cleavage. It is tentatively identified 
as an amphibole. Small ore grains and apatite 
are accessories. 


Porphyritic Olivine Trachydolerite 


This is the principal rock occurring in lenses 
of lava in the upper Huelster formation; it 
ranges in some places to a labradorite micro- 
syenite. Trachydolerite also makes one lens at 
the base of the Star Mountain rhyolite south 


of Seven Springs (Table 1, No. 14). It is a fine- 
grained black rock with a finely crystalline 
groundmass and small phenocrysts which are 
rare in one occurrence (Table 1, Nos. 3 and 4). 
A few grains of mafic mineral are distinguishable 
megascopically. The phenocrysts are narrow, 
dark, shiny laths of plagioclase (Ang2) commonly 
having an armor of anorthoclase. Labradorite is, 
for the most part, the dominant feldspar, but 
in some cases the alkali feldspar forms the bulk 
of the matrix as well as heavily armoring the 
phenocrysts. The groundmass feldspar grains 
are 0.02-0.3 mm. across, and the texture is 
hypautomorphic granular to trachytic with 
augite ophitic to the feldspar. The augite is 
pale gray to lavender with optically continuous 
patches measuring up to 2 mm. in greatest 
dimension. Olivine, largely serpentinized, is 
present between the feldspar laths in grains 
averaging 0.2 mm. in size. Magnetite is prom- 
inent, and some slides show brown biotite 
formed on the peripheries of a few grains. 
Apatite is a common accessory, and traces of 
calcite and analcime are present interstitially. 


Basalt and Trachybasalt Porphyry 


This is the type found in the flows of lava 
member No. 2 of the Seven Springs formation. 
It is dark gray to dull purple brown; some hand 
specimens are amygdaloidal, and a few show 
narrow plagioclase phenocrysts. Flow structure 
is prominent, and vesicles and amygdules are 
numerous. The groundmass ranges from glass 
containing barely visible plagioclase laths to 
plagioclase laths 0.01-0.3 mm. long, generally 
prominently armored by alkali feldspar. One 
specimen (Table 1, No. 24) is sufficiently coarse- 
grained to be considered a trachydolerite. The 
groundmass plagioclase is slightly less calcic 
than the labradorite (Anss_¢9) phenocrysts. 
Alkali feldspar also occurs in the interstices 
between groundmass laths and along cracks in 
phenocrysts. The texture is ophitic. Abundant 
microlites of augite and ore make the glassy 
matrix opaque to pasty. Iddingsite grains after 
olivine, abundant in the matrix, are rare as 
phenocrysts. Apatite is present in fine needles, 
and chalcedony, calcite, and serpentine (?) fill 
opal-lined amygdules. 
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Vitric Rhyolite Porphyry or Vitrolithic Tuff 


Several lavas are partially vitric or vitrolithic 
tuff. These include a lava in the top of tuff 
member No. 1, the red lava at the base of lava 
member No. 2 (or top of tuff member No. 2), 
and a lava in tuff member No. 4, all of the 
Seven Springs formation. They are brick-red to 
purple-brown dense flinty rocks with notable 
glassy feldspar phenocrysts (usually anortho- 
clase) or fragments visible in hand specimen. 
The groundmasses are usually glassy, n <cana- 
da balsam, with perlitic structure, or finely 
crystalline with prominent flow structure. The 
matrices enclose euhedral and corroded pheno- 
crysts of alkali feldspar and quartz, and a few 
pieces of aegirite-augite, biotite, and magnetite 
along with inclusions of rhyolite and basalt. 
Broken fragments of alkali and plagioclase 
feldspar and quartz are common as are amyg- 
dules and stringers of calcite, chalcedony, or 
hematite. Hematite is also prominent in the 
groundmasses. The tuffaceous parts are similar 
in composition, except the major constituent 
is glass in the form of shards which occasionally 
exhibit a somewhat welded character. 
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MECHANICAL BASIS FOR CERTAIN FAMILIAR 
GEOLOGIC STRUCTURES 


By M. Kinc HvuBBERT 


ABSTRACT 


A simple experiment with loose sand shows that this material exhibits faulting under deformational 
stresses in a manner remarkably similar to rocks. Moreover, the sand experiment is amenable to 
theoretical analysis with good agreement between predicted and observed behavior. The same theo- 
retical treatment, with slight modification, is alsoa pplicable to the behavior of rocks, and appears to 
afford a basis of understanding for a variety of empirically well-known geologic structures. 
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INTRODUCTION 


During the last century and a half, a large 
volume of empirical data on the mechanical 
behavior of rocks in response to tectonic stresses 
has been accumulated in the course of geologic 
mapping over all parts of the earth. In addition 
geologists have given considerable attention 
to the relation between various well-defined 


classes of geologic structures and the associated 
forces or stresses. This has included both small- 
scale structures such as slaty and fracture 
cleavage and joints, and large-scale structures 
such as faulting and folding. During the same 
period in the fields of physics and engineering, 
a great deal of exact knowledge has been 
acquired concerning the properties of stress in 
extended media and the response of matter to 


355 








356 


various stress situations, but only a beginning 
has been made toward an amalgamation of this 
knowledge with geologically well-known situ- 
ations. 

The present paper will examine some of the 
more common types of large-scale geologic 
structures in the light of physical relations 
which have to be satisfied, to see whether a 
better understanding of these structures may 
thus be acquired. Historically, the paper is 
based upon a series of lectures in structural 
geology delivered at Columbia University dur- 
ing the period 1931-1940. For his original 
interest in the geological aspects of the problem, 
the author is indebted to Professors J Harlen 
Bretz and Rollin T. Chamberlin; for the experi- 
ment to be described and the applicable theory, 
he is principally indebted to the writings of, and 
personal conferences with, Dr. A. Nadai (1928, 
1931a, 1931b). 

The essential theory appears to have first 
been established by W. J. M. Rankine in 1857 
in a study of the conditions of stability of 
loose earth. This was elaborated in a series of 
papers by Otto Mohr (1871; 1872; 1882; 1900) 
who, in addition to developing simple graphical 
methods of stress analysis, also enquired into 
the conditions of failure. Within recent dec- 
ades, these principles have been extensively 
employed in the newly developed subject of 
soil mechanics. In geological literature it was 
noted by Chamberlin and Shepard (1923, p. 
511-512), in some model experiments on geo- 
logic structures, that loose sand exhibited re- 
verse faulting, with dip angles ranging from 
20° to 35°. W. J. Mead (1925) reported further 
upon this property as related to dilatancy, but 
without analysis otherwise. Recently E. M. 
Anderson (1942) has reverted to the original 
work and methods of Rankine and has reached 
by a somewhat different route conclusions in 
essential agreement with those presented here. 


EXPERIMENTAL FAULTING IN SAND 


Experimental Apparatus 


Plate 1, Figure 1 shows a box with a plate- 
glass front which is divided into two compart- 
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ments by a rigid partition. The partition is 
attached to a longitudinal screw which can be 
turned by a crank. When the screw is rotated, 
the partition is moved by translation to the 
right or the left, parallel to the long axis of the 
box. 

The box is filled to a depth of about 23 inches 
with a fine-grained, loose, dry sand obtained 
from the beach at Galveston. Embedded in the 
sand next to the plate glass are three horizontal 
white lines composed of dry powdered plaster 
of Paris placed on the sand when the box was 
being filled; except as visible markers, they are 
without significance. 

The sand in the box will stand at a maximum 
angle of repose of about 30° on a free slope. 

When the partition is moved from left to 
right the left-hand compartment is lengthened, 
the right-hand compartment is shortened by 
the same amount and since the behavior of the 
sand in one is independent of that in the other, 
two separate experiments are performed si- 
multaneously. 

The first thing that happens (PI. 1, fig. 2) is 
the development in the left-hand compartment 
of a distinct normal fault whose measured 
angle of dip is 63°. In the meantime, nothing 
other than a slight bulging is observable in the 
right-hand compartment. As the partition is 
cranked farther to the right, reverse faults 
(Pl. 2 figs. 1, 2) with a measured angle of dip 
of about 28° develop. These extend from the 
bottom of the box near the foot of the partition 
to the surface of the sand where clearly visible 
escarpments are produced. This experiment 
has been repeated a number of times and the 
results are presented in Tables 1 and 2. 


Average Dip of Normal Faults 


In the case of the normal faults, the fric- 
tional drag of the sand on the glass plate 
invariably produced a steeper angle at the 
plate than in the interior of the box. Since the 
fault plane extends from the scarp on the 
surface to the foot of the partition its angle 
in the interior can be determined approx- 
imately by the equation, 
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FicurE 1. EXPERIMENTAL APPARATUS 
Box is filled with loose, dry sand. Markers are powdered plaster of Paris. 
Ficure 2. OBLIQUE ViEW oF NorMAL FAULT SHOWING ESCARPMENT 
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a = tan 2/z (1) 


where x is the distance of the escarpment from 
the partition in the interior of the box and z 
the depth of the sand. 


TaBLE 1.—ANGLES OF Dip oF NorMAL FAvLts 






































| Distance z Angle of Dip 
Experi- Depth : of fault 
ment | of sand | scrap from 
No. (inches) | partition Interior At glass 
| | (inches) of sand surface 
1 2.56 1.4 61.3° 64° 
2 2.7 1.4 62.6 63 
3 2.65 1.3 63.9 67 
4 2.8 is 61.8 63 
5 2.65 1.6 58.9 60 
6 2.65 1.6 58.9 60 
7 2.6 1.4 61.7 62.5 
8 2.35 1.4 63.1 64.5 
9 2.65 1.25 63.0 62.5 
10 2.7 1.45 61.8 62 
il 2.6 3.5 60.0 60 
12 ye 4.7 57.8 59 
13 2.5 1.4 60.8 ot 
ee eee 61.2° 63.2° 








Four of the measurements for normal faults 
were made with loosely packed sand; the 
remainder were made after the sand had been 
shaken to a closer packing. The angles for 
loosely packed sand were somewhat less steep 
than those for the closely packed sand, and 
the fault planes were less sharply defined. 


Average Dip of Reverse Faults 


It was not possible to compute the angle of 
dip of the reverse faults in the interior, because 
the fault planes did not originate at the bottom 
edge of the partition, but at a variable distance 
infront. Angles of dip of the reverse faults 
as observed against the glass plate appear to be 
somewhat steeper than the angles in the interior 
of the box. This is indicated by the arcuate 
shape of the fault scarp which is convex in the 
direction of the thrusting. The difference, how- 
ever, is probably not more than a degree or 
two. 


COMPARISON OF RESULTS WITH 
GEOLOGICAL FAULTS 


No field geologist can observe this experiment 
without being impressed by the remarkable 


TABLE 2.—ANGLES OF Dip OF REVERSE FAULTS 

















IN SAND 
Experiment No. Angle of Dip 
i a 
r 22 
3 28 
4 23 
5 25 
6 28 
7 27 
8 26.5 
9 26.5 
10 29.5 
11 24.5 
12 22 
13 26 
21.5 
” 27.5 
24.5 
15 20.5 
16 27 
PO asc ct sane sianis. 25 .2° 








similarity between the faults produced in the 
sand and those which occur in rocks in the 
field. Although statistical evidence on the ob- 
served angles of dip of faults is sparse, it is a 
matter of common geological observation that 
in areas of medium topographical relief and 
mild tectonic deformation (where the deforma- 
tion is insufficient to seriously change the 
initial angle of dip of faults) normal faults occur 
with dips consistently greater than 45°—rarely 
less—and reverse faults have dips consistently 
less than 45°. 

This generalization is confirmed strikingly 
by a recent study by Sax (1946, p. 42-46) of 
faulting in the coal measures in the Nether- 
lands. In this study, measurements were made 
underground of the angles of dip of 2102 
separate faults. Of these, 1651 (79 per cent) were 





PiaTE 2. EXPERIMENTAL FAULTING IN SAND—REVERSE FAULTS 


FicurE 1. DEVELOPMENT OF REVERSE FAULTS 
Ficure 2. OBLIQUE VIEW OF REVERSE FAULTS 
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normal and 451 (21 per cent) were reverse 
faults (thrusts and overthrusts). Frequency 
curves in 5°—interval groups were also plotted. 
The curve for each class of fault showed a 
well-defined peak, that for normal faults at 
63° and that for reverse faults at 22°. 

In 1913, C. K. Leith (p. 55) reported on the 
same object as follows: 


“The dip of a thrust fault is usually low, that of a 
normal fault, high. An average from the United 
States Geological Survey Folios gives a dip of 36° 
for reverse fault planes and 78° for normal fault 
planes.” 


These figures are of only qualitative 
significance since without doubt many faults 
where placed on the cross sections for which 
accurate dip measurements were unobtainable. 

Many normal faults are known from sub- 
surface work in the Texas and Louisiana Gulf 
Coast region. In most of these, dips steeper than 
45° are reported, though some have dips as 
low as 35°. In this region, the faulting in most 
cases has occurred in incompletely compacted 
sediments which subsequently have become 
compacted by amounts up to as much as 50 
per cent of their initial thickness. Hence, an 
initial dip of 60° could be reduced to as little 
as 41°, or an initial dip of 50° to 31°. 

Billings (1946, p. 206) shows a line drawing 
of a graben structure made by Hans Cloos with 
soft mud as his experimental material. The 
principal normal fault even in this “incompe- 
tent” material has a dip of 50°. 

Thus, not only is there a qualitative re- 
semblance between the faults in the sand box 
and those in rocks, but the angles of dip in the 
two cases are in reasonable agreement. This 
strongly suggests that the mechanical basis for 
the phenomena in both cases is essentially the 
same. 


MECHANICAL BASIS FOR THE FAULTING 
Nature of Stresses Produced 


Consider the state of stress at a point in the 
interior of the sand in the right-hand compart- 
ment of the box (Fig. 1). In the undisturbed 
state, the vertical and horizontal normal stresses 
will be approximately equal, and the stress 
state will be hydrostatic with the normal stress 
o on any surface equal approximately to the 


pressure of the overburden, or 
o= Pp = pgz (2) 
where p is the bulk density of the sand, g the 
acceleration of gravity, and z the depth beneath 
the surface of the sand. 
As the partition is advanced to the right 
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FicurE 2.—Srates or Stress WHEN ParTITION 
Is MOVED TO THE RIGHT 








(Fig. 2), the vertical stress o, remains constant 
but the horizontal stress o, is gradually in- 
creased. From the geometry of the deformation, 
these two stresses will also be, respectively, 
approximately the minimum and the maximum 
stresses acting, so we may regard them as 
the minimum and maximum principal stresses 
o3 (minimum) and oa; (maximum). An inter- 
mediate principal stress, o2,' will exist at right 
angles to these two, but in a two-dimensional 
problem it does not require consideration. 
Starting with the maximum and minimum 
principal stresses o; and os, let us first enquire 
as to the nature and the magnitude of the 
stresses acting across a plane which is obliquely 
oriented with respect to the axes of greatest 
and least principal stresses. Let such a plane 
be parallel to the axis of the intermediate 
principal stress, and let it make an angle @ with 
the axis of the least principal stress 3. 


Stress on Arbitrary Plane 


Figure 3 is a section in the plane of the 
greatest and least principal stresses. Let this 
section be of unit thickness and let the line 
MN be the trace of the plane across which the 
stresses are to be determined. Let ABC be a 
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small right-triangular prism whose sides are 
respectively perpendicular to the axes of the 
stresses 0, and a3. Let dS be the area of its 
hypotenuse, and let o and r be respectively 
the normal and shear stresses which act upon 
this surface. 

We take the prism small enough that its 
weight is negligible as compared with the 
surface forces which act upon it. It will be 
accordingly in equilibrium, under the influence 
of the surface forces due to the stresses oi, 
g3,0, and 7. Resolving these forces into hori- 
zontal and vertical components, the sum of 
the vertical components must be zero, and also 
the sum of the horizontal components. 

The forces acting upon the prism are the 
products of the stresses and the areas acted 
upon. The areas of three sides are: 


side AB = dS, 
side BC = dS sin a, (3) 
side CA = dS cosa. 


The sums of the horizontal and vertical 


components are respectively: 


a SS (4) 


o;dS sina —asinadS +rcosadS = O. 


If we eliminate dS and solve equations (4) 
for the two unknowns 7 and a, we obtain 


o = 01 cos?a + o; sin’? a, (5) 
Tt = (co; — a3) sin a cosa. 


From trigonometry: 


1 + cos2a 
cot a = as 
: 1 — cos2a@ 
ava=/|————, 

2 
; sin 2a@ 
sina cosa = . 


When these values are substituted into equa- 
tions (5), those reduce to the simpler forms: 
= oi + o3 a — G3 


2 + 2 cos 2a, 





(6) 
r/ AF tin 2e, 


wherein 7 and o are expressed as simple func- 
tions of ¢; and a3 and the angle 2 a. 


Mohr’s Stress Circle 


These equations are capable of a very simple 
geometrical representation if we choose co- 
ordinates with 7 as the axis of ordinates and ¢ 
as the axis of abscissas. In such a system, o1 
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FIGURE 3.—StTRESSES Across A SURFACE MAKING 
ANGLE a WITH AXIS OF LEAST PRINCIPAL STRESS 


and o3 (Fig. 4) will be points on the o-axis. 
The quantity (0: + o3)/2 will be a point 
midway between oa; and o;, and (a, — o;)/2 
will be a length equal to half the distance be- 
tween o; and o3. 

Now if we let the point (6: + o3)/2 be an 
origin for a radius vector which makes an 
angle 2 @ with the positive c-axis, and if we 
let (0:1 — o3)/2 be the length of this radius 
vector, we shall find that for every value of 
the angle a the coordinates of the terminus of 
this vector satisfy equations (6) and hence 
give the values of o and 7 for the plane which 
makes an angle a with the axis of least stress. 
Also, since (6:1 — o3)/2 is constant, as @ is 
changed the radius vector describes a circle 
which is the locus of all of the values o and r 
for all possible orientations of the reference 
plane. 

This construction (Fig. 4), originally de- 
scribed by the German engineer Otto Mohr 
(1882) and known as Mohr’s circle, enables 
one to determine by inspection the magnitudes 
of both the normal and the shear stresses upon 
a plane of arbitrary orientation a through 
a given point when the principal stresses a; and 
o3 at that point are known. For example, when 
ais zero, ¢ = 01, and t = OQ; upon a surface 
perpendicular to the greatest principal stress, 
the normal stress is equal to the greatest prin- 
cipal stress and the shear stress is zero. 
Likewise, when the surface is perpendicular to 
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the least principal stress (a = 90°, 2a = 180°) 
o@ = 03,t = 0. We can also see that the shear 
stress is a maximum at 2a@ = +90°, ora = 
+45°, and that magnitude of the maximum 
shear stress is 


Tmaz = a (7) 


Another interesting observation pertains to 
the shear stresses upon any two mutually 
perpendicular surfaces a and a + 90°. The 
points on Mohr’s circle for these two planes 
will be at 2 a and 2 a + 180° and so will be 
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a given surface slippage will occur. We can 
get this information by a separate investigation. 
If we place the sand in a box (Fig. 5) whose 
upper half can be skidded with respect to the 
lower half, thereby producing slippage in the 
sand along the plane of shear, we can vary the 
normal stress ao by increasing the load N 
pressing upon the sand. For a given value of 
a, we can increase the shear stress 7 by increas- 
ing the force T until slippage occurs. 

This experiment upon sand has been per- 
formed many times by the students of soil 
mechanics. It is found, to a first approximation, 
that failure occurs when 
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FicurRE 4.—Mour’s Stress CIRCLE 
; : , E r 
diametrically opposite one another Conse a, Pn (9) 
quently N o¢@ 


Ta = — Tota; (8) 


which indicates that the two shear stresses will 
be of equal magnitudes but of opposite sense- 


Properties of the Material 


We now require some additional knowledge 
of the physical properties of the material. We 
have observed that the material (sand) fails 
by slippage along discrete surfaces; it will be 
of interest, therefore, to consider under what 
conditions of tangential and normal stress upon 


This experiment is closely analogous to that 
for determining the friction between two solids. 
In the latter case a block of one solid is rested 
upon the horizontal surface of another and 
loaded by various amounts. For each load, the 
force required to slide the block is determined. 
In this case, too, it has been found that for 
given materials 
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The forces T and N or the stresses 7 and ¢ 
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can at failure can be thought of as being the two about 30°, it follows that this is the steepest 
on. components of a force or a stress, respectively, angle at which the slope is stable (Fig. 6). 
ose acting across the surface of contact at an 
the angle ¢ to the surfacenormal which isa constant Fracture Lines on Mohr Diagram 
the of the material. Hence 
the T + Let us now combine our knowledge of the 

N 3e-* tan ¢. (10) property of the material as expressed by equa- 
of : tion (10) with the general properties of stress 
as- In ordinary friction, @ is said tobe the angle of as expressed by Mohr’s diagram. Upon the 

friction between the given materials; in con- 7o-plane of Mohr’s diagram (Fig. 7), equation 
er- =f tinuous materials, such as the sand here dis- (10) representsa straight line of slope ¢ passing 
oil cussed, @ is said to be the angle of internal through the origin. From symmetry there will 
on, friction of the material. Extensive tests (Kry- be two such lines, one with a slope of + ¢ and 
nine, 1941, p. 179) have shown that for sand the other with a slope of — @. These lines 
the angle @ is not strictly constant but ranges mark the boundaries between regions of stable 
N 
x 
’ 
if bs 4 . ‘ 2 < 7h oe ——_——_—_———— T 
~ % ai ace : 
7 e : ~~ om v | ° . ° 5 . 
EE ee is: ices niente 
. . . n ITZ ; ‘ ‘ e 
5 ae . / c. a : . 
" , *. an - *« ° SAND 
Ficure 5.—TeEst Box For MEASURING r/o RATIO AT WHICH SLIPPAGE OccuRS 
from about 30° to 35°. The angle of 30° is and unstable states of stress for the sand. The 
characteristic of loose sand, whereas the higher state of stress corresponding to any point in 
} angles occur when the sand is in a state of the region between these lines and the o-axis 
y close packing. has a value 
The lower angle of 30° is also the approximate 
9) angle of repose of loose sand. That this should * < tan¢, 
be so can be seen by noting that upon a slope of ” 
at inclination 8, the normal and tangential stresses 4nd hence is a condition of stability; for any 
1s, across the bottom of a surface lamina of sand point outside the lines 
od will be respectively 
nd o = wos 8, (11) > tan ¢, 
he Tt = wsin B, “4 
: or so that fracture will occur before this state can 
be reached. These lines may also be regarded 
= = tan B. (12) as fracture lines, since when the stresses upon 
. any surface reach the values represented by 





a 


Hence, when the slope is tilted until B = 
, slippage will occur. Since in this case ¢ is 


points on either of the two lines, fracture will 
occur. 


















Analysis of Reverse Faulting 


Combining now the fracture lines of Figure 
7 with Mohr’s diagram, we can anticipate what 
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FiGURE 7.—REGIONS OF STABLE AND OF UNSTABLE 
STRESSES ON THE or PLANE 


should happen in the sandbox experiment (Fig. 
8). In the initial stage 


o, = o3 = pgh. 


As the partition is pushed forward, the vertical 
stress o3 remains constant while the horizontal 
stress o; is gradually increased. For each value 
of o; a particular Mohr’s circle will correspond, 
of which o3 and a; will lie at opposite ends of 
the horizontal diameter. As o; is increased, the 
Mohr’s circle will get larger until finally it 
will become tangential to the two lines of 
fracture. At this point, fracture will occur 
upon one or more of those planes in the body 
whose trangential and normal stresses have 
reached the critical ratio. These will be the 
planes whose angles a correspond to the points 
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of tangency of the Mohr’s circle with lines of 
fracture. From the diagram the points of tan- 
gency occur at 


2a = + (90° + ¢) ) 


eazet («s+$). 


Since, as we have noted, @ for sand ranges 
from 30° to 35°, fracture should occur along 
one or both planes, making an angle of about 
60° to 62.5° to the least principal stress, which 
is vertical, or 27.5° to 30° to the horizontal. 
This is in fair agreement with the dip of the 
reverse faults shown in the photographs. It is a 
few degrees higher than the average dips ob- 
served for the reverse faults—a discrepancy to 
which we shall return later. 


(13) 


Analysis of Normal Faulting 


Let us now consider what happens in the 
compartment where the normal faults occur. 
Here, as in the reverse-fault compartment, the 
initial state of stress at a point inside the sand 
is given approximately by 


Oo} = 03 = pg. 


In this case, however, as the partition begins to 
recede, the horizontal stress, which now be- 
comes the least stress o3, begins to diminish, 
while the vertical stress, which is now the 
greatest stress o1, remains stationary, as be- 
fore. 

On Mohr’s diagram (Fig. 9) we have, in this 
case, a fixed maximum stress and a diminishing 
minimum stress, with again a Mohr’s circle of 
increasing radius. Ultimately, as o3 is further 
decreased, the circle will again become tan- 
gential to the fracture lines r = o tan ¢, and 
fracture will occur. As before, the angles a of 
the planes along which fracture will occur will 
be given by 


2a = +(90° + ¢), | 


or 
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MECHANICAL BASIS FOR THE FAULTING 
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FicurE 8.—REVERSE-FAULT CONDITIONS 


Stress o; increases until] Mohr’s circle reaches fracture lines. Fracture occurs at angle a = 45° + <=, 


so that in this case also the fracture should be 
at an angle of about 60° to 62.5° to the least 
principal stress. In this instance, the least 
principal stress is horizontal, so the fracture 
plane should dip about 60° to 62.5°. 

The average dip for normal faults in Table 
1 was 61.2°, which is a very satisfactory agree- 
ment. 

The normal-fault case is instructive in still 
another respect. Geometrically, normal faults 
correspond to a horizontal elongation of the 
section in which the faulting occurs. These 
faults accordingly have commonly been at- 
tributed to tensile stress and are often referred 
to as “tension faults.” Loose, dry sand has 
zero tensile strength and is incapable of sus- 
taining a tensile stress. We observe, however, 
that normal faults occur in it as readily as in 
solid rocks, and entirely as the result of compre- 
Sive stresses. In fact, in both cases the fracture 
occurs along planes making an angle of 45°— 
$/2 to the greatest principal stress. In the 
reverse-fault case this stress is horizontal; in 
the normal-fault case it is vertical. 





¢ 
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APPLICATION TO COHESIVE MATERIALS AND TO 
GEOLOGIC STRUCTURES 


Coulomb Equation 


So far our theory has been restricted to a 
noncohesive material; rocks, however, are co- 
hesive materials, so it remains to be seen how 
we may modify the theory to include these 
also. Since Mohr’s circle represents a general 
property of stress and needs no modification, 
we need only modify equation (10) which 
describes the behavior of the material. This 
modification consists only in adding a cohesive 
shear stress 7) to the stress 7 of equation (10) 
to overcome the cohesion or initial shear 
strength of the material. The equation then 
reads 


-—* = tan¢, 


(15) 
or 


rT =%+ 0c tan d. 
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This equation appears to have been first 
proposed by the French physicist Coulomb 
(1776), and, while it may not be rigorously 
correct, it serves as an excellent first approxi- 
mation to the behavior of cohesive solids, 
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FiGurRE 9.—NorRMAL-Fautt ConpITIONS 
Minimum stress os decreases until Mohr’s circle 
reaches fracture lines. 


including rocks. When plotted upon the or- 
diagram (Fig. 10), equation (15) represents a 
pair of straight lines with slopes + @ and — @ 
intersecting on the o-axis at a pointo = 
To/tan @. 

The geometrical relations (Fig. 11) of the 
stresses at the time of fracture in this case are 
the same as in the noncohesive case. Fracture 
occurs when the Mohr’s circle becomes tan- 
gential to the lines of fracture on the or-dia- 
gram, and, as before, the planes of fracture 
make an angle of 45° + ¢/2 to the axis of 
least stress, or of 45°— @/2 to the axis of 
greatest stress. 


Data upon Angle d for Rocks 


Two methods are available for determining 
the angle of internal friction, ¢, for rocks. One 
consists in computing @ from the observed 
angle between the surfaces of slip and the 
axes of principal stress when the latter are 
known; the second is by means of tests of rock 
specimens under known conditions of triaxial 
stress. 

Prandtl and Rinne (Nadai, 1931la, p. 109 
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BASIS FOR SOME STRUCTURES 


110) carried out a series of tests on polished 
cylindrical specimens of rock by compressing 
them axially to incipient failure. On the cy- 
lindrical surfaces, two families of helical lines 
of fracture were observed intersecting at an 
angle 2 @ bisected by the axis of greatest stress, 
The angle @ is then complementary to the 
angle a, and from the Mohr diagram (Fig. 11) 


?$=2a—9 = 9-28. 


On two photographs of marble specimens 
(shown by Nadai) the angle 2 6 was 41° and 
45°, respectively. The corresponding values of 
@ would therefore be 49° and 45°. 

The same principal may be used in con- 
nection with the dip angles of simple normal 
and reverse faults. When faulting occurs at or 
near the surface of the ground in nearly level 
topography, the principal stresses must be 
respectively nearly vertical and horizontal. Here 
also 


= 9 —2 6 


where @ is the angle of dip of reverse faults, or 
the angle of hade of normal faults. 

Thus from fault data, with the angles of dip 
of simple reverse faults, or the hades of normal 
faults lying between the limits 20° and 40°, the 
angle of internal friction would lie between 
10° and 50°. For the peak frequency in the 
angle of hade for normal faults of 27°, as 
reported by Sax, the peak frequency for the 
angle @ would be 36°. 

The second method, that of testing specimens 
under known conditions of triaxial stress, con- 
sists in testing a cylinder under simultaneous 
radial and axial stress. The radial stress is 
supplied by a fluid pressure and the axial 
stress by a piston. The fluid is set at a chosen 
pressure p = o2 = 03 (or p = a1 = 2) and the 
axial stress, 01, is increased (or decreased if the 
axial stress is o3) until failure occurs. The 
Mohr’s circle for the values of o; and a3 at which 
failure occurs is constructed. This procedure is 
repeated with additional specimens of the same 
material until a family of over-lapping circles 
of increasing radii is obtained. The envelopes 
of these circles obviously represent the em- 
pirically determined loci of the fracture curves 
of the diagram. The angle of slope of the curves 
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with respect to the a-axis is the angle of internal 
friction for the given material. 

Tests of this kind appear first to have been 
made on rocks by von K4rmén (1912); he 


T 


few results on rock and concrete have now been 
published (Jones, 1946; Balmer, 1946; Mc- 
Henry, 1948). On rock cores from Boulder 
Dam (rock type not specified), the results of 
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Ficure 10.—PosiTIoN OF THE FRACTURE LINES ON or DIAGRAM FOR COHESIVE MATERIALS 


ran a series of tests on marble which was found 
to behave as a brittle material at lower radial 
pressures but to flow plastically at higher 
pressures. The envelopes of the Mohr’s circles 
thus obtained were convex outward with 
varying from an initial value of about 40° 
(corresponding to @ = 25°) to a final value of 
near zero. Also conjugate slip surfaces on the 
slightly deformed specimens show an acute 
angle 2 8 = 50°. 

Triaxial tests on concrete were reported by 
Ros and Eichinger (1928). The angle @ as 
observed directly varied from 22° to 33°, which 
would correspond to values from 24° to 46° 
for the angle ¢. 

The U. S. Bureau of Reclamation at Denver 
has recently installed a large triaxial testing 
machine for testing cores of rock and concrete 
in connection with the construction of large 
dams (Blanks and McHenry, 1945), and a 


16 tests gave for @ the value of 35°, which 
would correspond to @ = 27.5°. Tests on 
concrete by both Jones and McHenry have 
given values of ¢ & 45°, or of 0 = 22.5°. 

It is interesting to compare these results 
with tests on asphalt aggregates for road build- 
ing recently reported by Nijboer (1942). Tri- 
axial tests were made upon dozens of different 
mixes under a wide range of conditions. In most 
cases the envelopes of the Mohr circles were 
straight lines though occasionally one with a 
small outward convexity was observed. 

The angles @ found by these tests ranged 
between the extremes of 19° and 38°. The most 
frequently observed value, however, was 28°, 
and, for more than 95 per cent of all the 
measurements made, the observed values lay 
within the range 28° + 4°. For rocks, this 
would correspond to normal faults with hades, 
or reverse faults with dips, of 31° + 2°. Kry- 
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FicurE 11.—Srress CONDITIONS FOR FAULTING IN COHESIVE MATERIALS 


Faults occur at angle a = 45° + : to the axis of least stress, 
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FIGURE 12.—ConpDITIONS FOR ZERO HorizONTAL STRESS IN NORMAL FAULTING 
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APPLICATION TO COHESIVE MATERIALS 


nine (1941, p. 172) reports that, for wet clay, 
@ ranges from a few degrees to 15° or 20° 
This would correspond to normal-fault dips 
between 45° and 55°. 

Thus, although data upon the angle of in- 
ternal friction for rocks are as yet sparse, 
the computed angles of faulting are in good 
agreement with field observations. 


Tension in Normal Faulting 


We have already seen that, in the case of 
cohesionless materials, normal faults are pro- 
duced only under compression. For cohesive 
materials (Fig. 12), it will be seen that, for any 
given material, there must exist at the time of 
faulting a certain critical value of o1, say o:* 
for which og will be zero. Then, for all values 
of 0; > o:*, the horizontal stress ¢3 in normal 
faulting will be compressive, and for all values 
of «1 < o:* the horizontal stress a3 will be one 
of tension. 

Since, in the case of normal faulting, the 
vertical stress at a given depth is equal to the 
pressure of the overburden, there will be some 
critical depth for a given material below which 
the horizontal stress producing normal faults 
will be compressive and the existence of a 
tensile stress will be impossible. At less than 
this depth, the horizontal stress will be one of 
tension. 

The approximate depth at which this transi- 
tion will occur can be obtained from crushing- 
strength data of rocks. The conditions 


a= 0, 
“a = o*, 


are those which are satisfied approximately in 
the crushing-strength tests for rocks at at- 
mospheric pressure. Consequently, for a given 
rock the critical stress o,* is equal approxi- 
mately to its crushing strength. For dry rocks 
the depth at which this vertical stress occurs 
would be 


* 


2* = o,*/pg 


where p is the bulk density of the material. 
In the earth, however, the rocks ordinarily are 
Saturated with ground water, which exerts an 
upward force of buoyancy. In this case, the 
effective stress is that transmitted by the 


367 


solid system. Hence this effective stress ¢, at 
a given depth is 


a = (ps — pw) g (1 — €) 2, 
or for the critical depth 
2* = 3,*/[(o. — pw) g (1 — €)], 


where p, is the grain or mineral density p, the 
density of water, and € the porosity of the rock. 

The crushing strengths of rocks range from 
values near zero, for recent sediments, to as 
high as 3 X 10° dynes/cm? for granites. Such 
tests, however, are made upon nearly flawless 
specimens, so that, for large volumes of rock 
broken by numerous joints, the effective values 
will no doubt be considerably smaller. The 
Beaumont clay underlying Houston for ex- 
ample has a bearing strength of only about 
10® dynes/cm?. Recorded crushing strengths of 
sedimentary rocks are commonly of the order 
of 5—15 X 108 dynes/cm?. 

Introducing these values of o:* into equation 
(16) and solving for z*, we find the maximum 
depths at which tension can exist in rocks of 
various types. In the Gulf Coast sediments, for 
example, this depth is found to be of the order 
of 10° cm., or 10 meters (about 30 feet). In 
the more consolidated sediments, it may in- 
crease to a few hundred meters; and, in the 
strongest crystalline rocks, it may reach values 
of several kilometers. 

Yet the Gulf Coast sediments are cut by 
numerous normal faults extending from the 
surface to depths of thousands of meters. It is 
clear, therefore, that aside from a surface 
veneer, these sediments have been broken by 
a series of normal faults in response exclusively 
to compressive stresses. 


(16) 


PROBLEM OF ASYMMETRICALLY FOLDED AND 
FAULTED BELTS 


One-sided Thrust Hypothesis 


With this background, let us consider an 
equally familiar but more complex problem. 
One of the commoner types of orogenic struc- 
ture is the belt of asymmetrically folded and 
faulted sediments as exemplified in North 
America by the Appalachian Mountains, the 
Ouachitas, the foothills of the Rocky Moun- 
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tains, and some of the West Coast ranges, and in 
Europe by the Alps, the Juras, and the Scottish 
Highlands. In these cases, the folds are uni- 
directionally overturned and accompanied by 
reverse faulting, with the direction of thrusting 
—that is, the direction of the relative dis- 
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A pplication of the Newtonian Laws of Motion 


That this picture is somehow deficient be- 
comes evident when considered in the light of 
the Newtonian laws of motion, according to 
which the sum of all the forces acting upon the 




















FiGurE 13.—CoMpRESSIVE STRESSES ON ENDS OF SECTION OF ASYMETRICALLY 
FOLDED AND FAULTED SEDIMENTS 


placement of the upper block with respect to 
the lower—in the same sense as the overturning 
of the folds. Moreover, where data are avail- 
able, as in the foothills region of Alberta where 
seismic surveys have been made and numerous 
oil wells drilled, the fault surfaces have been 
found to be concave upward, the dip decreasing 
with depth. 

Considerable attention has been given to the 
causes of these structures, and the hypothesis 
most often favored has been that formulated 
by J. D. Dana (1847a; 1847b)—that such 
structures are formed by a one-sided “active 
thrust” from the direction opposite to that of 
the overturning. Thus, the Appalachian Moun- 
tains would have been formed by an active 
thrust from the southeast, which presumably 
died out with distance toward the northwest. 
An auxiliary question often considered is how 
far rocks of a given strength are capable of 
transmitting such a thrust. 

The state of stress thus postulated and its 
associated orogeny is illustrated in Figure 13, 
where a large “active thrust” ¢., acts on the 
left-hand end of the section and dies out to the 
much smaller stress o,, on the right-hand end. 
Not only does this appeal to the intuition, but 
it is supported by the deformational evidence 
as well. Since the rocks on the left are intensely 
deformed whereas those on the right are un- 
distorted, the inference is justified that o., is 
very much greater than ¢;,. 


block in a given direction must equal the 
product of the mass of the block and the 
component of its acceleration in that direction. 
Taking the x-axis as horizontal, this requires 
that in the x-direction 


21 
=F, = [ (oz; —oz2)dz = maz (17) 
0 


As an order of magnitude, (¢., — o:,) may 
be taken to be about 10*® dynes/cm?. Then 
if the block of Figure 13 is assumed to be 100 
kilometers wide, 10 kilometers deep, and 1 
centimeter thick, =F, will be about 10% dynes. 

For the ma-term, the mass of the block will 
be about 3 X 10" grams. Velocities in known 
orogenies are rarely more than a few centi- 
meters per year, so that secular accelerations 
much greater than a centimeter per year per 
year or about 10-® cm/sec? are unlikely. Hence, 
the ma-term would have a magnitude only of 
the order of 10-? gm cm/sec?, which is in- 
finitesimal as compared with the forces acting, 
and the sum of all forces must effectively be 
zero throughout the orogenic process. 

Since, with the stress distribution postulated 
in Figure 13, the sum of the horizontal forces 
acting upon the block is manifestly not zero, 
some essential element must be lacking. Since 
body forces are ineffective horizontally, the 
only remaining possibility is shear stresses 
upon the vertical faces parallel to the plane of 
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the drawing and along the bottom of the block. 
In arcuate structures, shear stresses on the 
vertical faces normal to the strike exist, but 
they vanish as the curvature approaches zero 
and the structure becomes rectilinear. 
Assuming the structure to be rectilinear, 


of opposite rotational senses. Hence there must 
also be a downward-directed shear stress on 
the left-hand vertical face of the block equal 
in magnitude at the bottom corner to the 
stress on the bottom, and approaching zero 
toward the top. A similar shear stress, but of 














FicureE 14.—CoMPLeTe SysTEM OF Two-DIMENSIONAL STRESSES ACTING UPON BOUNDARIES OF BLOCK 


we are then left only with the bottom shear 
stress Tzz, Whose magnitude at any given 
distance x along the bottom must just be 
sufficient to balance the force due to the stress 
difference do, in the horizontal distance dx. 
Thus, at distance x 


21 
T2dx = + [ [ (« + Ove is) _ o | 
0 Ox 


or 





21 a is 
To. = + [ ss dz; (18) 
0 x 


that is the shear stress on the bottom of the 
block at a given point is proportional to the 
gradient of the normal stress in that direction. 
The gradient of the normal stress would be 
approximately the same as the gradient of the 
intensity of the deformation. Hence, in Figure 
13 the shear stress would have a maximum 
value near the lower left-hand corner of the 
block and would decline to near zero at the 
lower right-hand corner. The complete equa- 
tion of forces in the x-direction must therefore 
be 


2) z2 
zF, | (ez, = O22) - i T12dX = 0. (19) 
0 21 
As noted, the components of shear stress 


parallel to a plane normal to both are the same 
on any two mutually perpendicular planes, but 


magnitude near zero and directed upward, 
must exist on the right-hand end. 

In addition, all surface forces and body 
forces must be so related that all turning 
moments acting on the block are zero. When 
these several conditions are satisfied the ap- 
proximate set of boundary stresses shown in 
Figure 14 is obtained as being necessary to 
produce the generalized type of orogeny shown. 


Interior Stress Distribution 


From the boundary stresses shown in Figure 
14, the approximate pattern of the stress dis- 
tribution in the interior of the block can be 
inferred. Trajectories of principal stress are 
curves which at every point are tangent to a 
given principal stress. Since the principal 
stresses are mutually perpendicular, then three 
perpendicular stress trajectories, one each for 
01, 02, and o3, must pass through each point. 
In a two-dimensional stress system on a plane 
parallel to a; and a3, two families of orthogonal 
stress trajectories, one everywhere parallel to 
o, and the other to 2, will exist. 

Recalling that, when the principal stresses 
are unequal, shear stresses exist upon all sur- 


1 The reasoning here is intentionally only quali- 
tative, the objective being to convey an intuitive 
sense of the principles involved. A more formal 
quantitative analysis is given in a companion 
paper by the author’s colleague, W. Hafner (1951). 
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faces except those perpendicular to the stress 
trajectories, the approximate pattern of the 
stress trajectories in the interior of the block 
may be determined from the combination of 
normal and shear stresses on its boundaries. 


% 





Ficure 15.—DrrecTions OF SHEAR STRESSES IN 
VaRI0OUS QUADRANTS DEFINED BY 
PRINCIPAL STRESSES 
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can be definitely known. Across any plane on 
which a compressive normal stress o and shear 
stress 7 exists, the axis of greatest principal 
stress, 01, will lie in the quadrant between the 
plane normal and the direction from which the 
shear stress occurs (Fig. 15). 

The boundary stresses of Figure 14 will 
produce the approximate pattern of principal- 
stress trajectories in the interior shown by the 
solid curves of Figure 16. The family of greatest- 
stress trajectories, 01, is tangential to the upper 
surface of the block, plunging downward and 
divergent to the right, and convex upward. 
The family of least-stress trajectories, 03, is 
orthogonal with this. 


Surfaces of Potential Faulting 


If we now assume that faulting is most 
likely to occur along surfaces tangent to the 
intermediate stress, o2, and at an angle of 


about 45° $ to the greatest stress, where @ 





FIGURE 16.—TRAJECTORIES OF PRINCIPAL STRESSES (SOLID LINES), AND OF POTENTIAL REVERSE- 
Fautt SurFACES (BROKEN LINES) COMPATIBLE WITH THE BOUNDARY STRESSES OF 
FIGURE 14 : 


Arrows indicate directions of potential slip 


On the upper surface of the block, the normal 
stress is the atmospheric pressure and the shear 
stresses are zero. Therefore, at any point upon 
this surface, one of the stress trajectories must 
terminate perpendicularly, and the other must 
be tangential to the surface. 

Across the left-hand end of the block, and the 
bottom, the stress trajectories must cross ob- 
liquely, with the angle of obliquity approaching 
zero near the upper surface, and increasing 
downward. Without knowing the actual magni- 
tudes of both the normal and shear-stress 
components, the angles of incidence of the 
stress trajectories can be determined only quali- 
tatively, but the quadrant in which each falls 


may be taken to be about 30°, we are able to 
sketch in on Figure 16 the traces of the sur- 
faces of potential faults, and indicate their 
directions of relative displacement. These com- 
prise two families of curves (shown by dashed 
lines) which everywhere intersect the o:-tra- 
jectories at an angle of about 30°. In virtue of 
the divergence and curvature of the stress 
trajectories, these surfaces are also curved, one 
set convex and the other concave upward. 

At any given point, the stresses on each of 
these conjugate surfaces are the same; so that 
if they are at the critical value for fracturing 
on one surface, the same is true for the other. 
However, slippage along a finite surface in- 
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volves an integral of conditions along the 
entire surface. Whatever this integral should be, 
in homogeneous materials with symmetrical 
stress distributions, its value over each of a 
pair of conjugate surfaces of potential slippage 
will be the same, so that equal slip should 
occur on each family of surfaces. 

In asymmetrical cases, the integral of any 
stress quantity over one surface or family of 
surfaces will in general be different from that 
over the conjugate surface or family. Because 
of this inequality, slippage should occur on one 
set before on the other, or on only one of the 
two sets. In the asymmetrical geological system 
here considered, faulting develops almost ex- 
clusively on the concave-upward set of sur- 
faces. The same was true for the reverse 
faulting in the sandbox experiment. 

In experiments with symmetrical systems, 
slippage on both sets of surfaces is observed. 

Returning to the sandbox experiment, it will 
be recalled that the average angle of hade for 
the normal faults was about 29°, and of dip for 
the reverse faults was 25°, whereas our theory 
indicated that the two should be the same. 
Perhaps we can now explain this discrepancy. 
We assumed the principal-stress trajectories 
to be strictly horizontal and vertical. However, 
in the reverse-fault case, the sand was pushed 
along the bottom which exerted a frictional 
reaction in the form of a bottom shear stress. 
The effect of this would be to deflect the 
horizontal stress trajectories downward and 
thus to lower the dip angle of the reverse 
faults. 

Other shear stresses along the vertical faces 
of the partition must also have existed which 
would likewise tend to produce some deviation 
in observed angles from the idealized con- 
ditions assumed. 


CONCLUSION 


It thus appears that some of the more 
common large-scale geological structures are 
in satisfactory agreement with theoretical de- 
ductions based upon the stress patterns in- 
ferred from the observed deformation and New- 
ton’s laws of motion, and the empirically 
determined properties of rocks under known 
stress conditions. 
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STRESS DISTRIBUTIONS AND FAULTING 
By W. HAFNER 


ABSTRACT 


Tectonic deformations result from a condition of internal stress caused, in turn, by primary and second- 
ary forces. In the geological literature, a great deal of discussion is based on a direct connection between 
forces and deformation, completely by-passing the concept of stress. This paper is a contribution in the 
intermediate field of stress relations. It presents the complete solutions of certain stress systems caused by 
various forms of boundary forces. Furthermore, the location and attitude of the fault surfaces likely to be 
associated with them is determined. 

The basic concept of stress is briefly reviewed and some of the fundamental differences between the 
force-vector and the stress-tensor are pointed out. The fallacy of applying the familiar methods of vector- 
addition of forces to problems in stress is demonstrated. 

For certain systems of external boundary forces acting on a portion of the earth’s crust, the internal 
stress distribution can be calculated by means of the familiar equations of elasticity. Appropriate calcula- 
tion methods for two-dimensional cases are shown and the basic equations applicable to a series of important 
boundary conditions are derived. The examples here presented include: (1) superposed horizontal compres- 
sion with constant lateral and vertical gradients; (2) horizontal compression with exponential attenuation; 
and (3) sinusoidal vertical and shearing forces acting on the bottom of a block. The latter equations provide 
solutions for differential vertical uplift and for the important case of drag exerted on the bottom of the 
crust by convection currents in the substratum. Diagrams show configuration of the stress trajectories 
and distribution of the maximum shearing stress for the resulting stress systems. 

A parallel series of diagrams shows the disposition between the relatively stable and unstable segments 
of the blocks and the probable attitude of the fault surfaces likely to be associated with the individual 
stress systems. The construction of the fault surfaces is based on the original stress distributions alone, 
the influence of local stress alterations due to the occurrence of fracture being disregarded. The full effect 
of this inter-action is not known, due to the extreme complexity of the problem. The fault patterns shown 
are strictly applicable only to the initial stages of fracture, but may also represent fair approximations dur- 
ing the more advanced stages, since the original stress remains the dominating influence and stress-altera- 
tions due to faulting diminish rapidly with distance. ; 
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INTRODUCTION 


The earth’s crust is unceasingly undergoing 
deformations because of the interplay of forces 
acting upon it. The primary forces are suffi- 
ciently powerful to produce disturbances of 
enormous magnitude, such as mountain chains 
or the crustal down-buckling deduced from 
gravity observations; the secondary ones tend 
to bring about partial or complete restoration 
of equilibrium. 

Many hypotheses have been advanced re- 
garding the origin and nature of the deforma- 
tive forces but, as yet, all these theories are 
speculative in varying degree and a wide di- 
vergence of viewpoint exists. Knowledge of 
the results of these processes is much more 
advanced, at least with respect to the upper- 
most portion of the crust. A vast and ever- 
growing body of observational data on the 
nature of the deformations is becoming avail- 
able. In fact, the validity of the various hy- 
potheses on the origin of the deformations 
rests largely upon the degree of agreement with 
this accumulated store of data, aside, of course, 
from the inherent soundness of the principles 
employed. 

In all quantitative studies on the relationship 
between the original forces and the resulting 
deformations, an intermediate field of investi- 
gation enters, the condition of stress in the 
earth’s crust. The original forces—whether pri- 
mary or secondary, whether internal body 
or external boundary forces—set up a state 
of stress in the earth’s crust which, in turn, 
produces the observed distortions. In the geo- 
logical literature, we quite often find arguments 
attempting to explain the deformations directly 
in terms of the original forces. While this 
approach may be justified in certain simple 
qualitative considerations, it is an inherently 
fallacious procedure and more often than not 
leads to confusion, rather than elucidation, of 
the problem discussed. It is especially inap- 
propriate in quantitative analysis. 


This paper offers the results of some studies 
in the intermediate field of stress relations. It 
presents the complete solutions of the stress 
distributions produced by certain systems of 
external boundary stresses. The nature of the 
faulting most likely to be associated with them 
is analyzed, although certain complications for 
which no complete solutions are available have 
to be disregarded in this phase of the investi- 
gation. The stress distributions here analyzed 
were selected with the view of covering some of 
the geologically more important groups of 
boundary stresses acting on portions of the 
earth’s crust, such as lateral compression, dif- 
ferential vertical uplift, and drag exerted on 
the bottom due to viscous currents in the 
substratum. No attempt is made to enter the 
more speculative field of the origin of these 
forces. However, the author believes that any 
advance made in our understanding of the 
relation between the observed tectonic phe- 
nomena and the causal stress systems will 
indirectly promote progress in the more im- 
portant, but much more illusive, field of the 
ultimate causes of the earth’s deformations. 
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STRESS AND FORCE 


In the geologic literature we encounter not 
infrequently a misunderstanding of the concept 
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of stress. Many authors use the terms “force’’ 
and “stress” synonymously, in the sense of 
dynamical forces, and are not aware of the 
fundamental differences between the two con- 
cepts. In view of this fact, it appears advisable 
to review briefly the general properties of the 





FicureE 1.—LiwitepD DEFINITION OF STRESS 


state of stress and its relation to forces. Such 
considerations, though of elementary character, 
may be helpful to readers not thoroughly 
familiar with these concepts. 

Stress is force per unit area. A more precise 
definition (Timoshenko, 1934, p. 3) states that 
it represents the limiting value of the ratio 
of a surface force divided by the area upon 
which it acts (Fig. 1); in formula: 


5F 
S = lim — 
im 3A 


This force may be resolved into a normal com- 
ponent o,, and a tangential component 7, . 
Depending upon whether ao, is directed outward 
or inward, it represents a tension or a pres- 
sure, respectively. The tangential component 
t,, Which is acting parallel to the surface 
element, is a shearing stress. 

This definition brings out one of the differ- 
ences between stress and force: that the two 
quantities are dimensionally at variance. The 
dimensions of force are (m, I, t-*); those of 
stress are (m, 1-1, t-*). Hence, the statement 
that stress is force, as found in a well-known 
textbook on structural geology, is incorrect. 
However, the difference between stress and 
force goes much deeper. If we are dealing 
with dynamical forces only, it would always 
be possible by the process of resolution of 
forces to determine the normal and tangential 
components with respect to amy direction of 
the surface element 8A passing through the 
point. When dealing with stress, on the other 
hand, we know nothing whatever of the mag- 


nitude of stress after rotation of the element 
into a different direction. To solve this simple 
problem, we need an expression for the state 
of stress which defines this quantity for every 
surface element leading through a given point, 
irrespective of direction. Mathematically, this 
more comprehensive expression is a tensor 
quantity and, more specifically, a symmetric 
tensor of second degree. This is in marked 
contrast to force, which is a vector quantity. 
A tensor requires for its description more than 
three components and is therefore a more 
complicated concept than a vector. The rela- 
tion is thus similar to that between a vector 
and a scalar. It is as misleading to think of 
stress in terms of a vector as it would be to 
treat a problem involving force in terms of a 
scalar function. The specific properties of the 
stress tensor can be adequately described with- 
out recourse to the general theory of tensors. 
The most important fact is that a total of six 
quantities is needed for a complete determina- 
tion of the state of stress at any one point in 
a solid. In contrast, force, being a vector, 
requires only three quantities for specification. 
This explains why the above simple problem 
is indeterminate if only the three parameters 
of a force-vector are given. 

There are two convenient ways to describe 
the six parameters of the stress tensor (Fig. 2): 
(1) State the stress components acting on the 
surface of a cubic element whose faces are 
parallel to the coordinate axes or (2) Use the 
well-known stress ellipsoid. The stress ellipsoid 
is defined by three mutually perpendicular 
directions, called the principal directions, and 
the intensities of the stresses in these directions, 
called the principal stresses. Each method re- 
quires six independent quantities; if known at 
each point, the condition of stress throughout 
the body is completely determined. The stress 
systems analyzed in the present paper will be 
restricted to two-dimensional cases of plane 
deformations. These are characterized by the 
fact that variable displacements occur only in 
parallel planes, namely those planes which 
are oriented perpendicular to the axes of the 
structures. Elongation in the third, or longi- 
tudinal, direction may or may not take place; 
if it does, it is constant over the entire cross- 
section. In such systems, two of the independent 
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shear stress components are automatically zero, 
and the normal stresses parallel to the longi- 
tudinal axes do not vary in this direction. The 
solution of stress problems is then reduced to 
the determination of three components instead 
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a) Stress components on cubic 


element. 
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simultaneously at a point, the resultant force 
may be obtained by the well-known method 
of vector addition. Not infrequently we find 
this method applied to stresses too, sometimes 
quite specifically, more often in rather vaguely 
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b) Stress ellipsoid. 


FIGURE 2. COMPLETE EXPRESSIONS OF STRESS TENSOR 


of six, as in the general three-dimensional case. 
Taking the x-axis as the horizontal direction 
in the cross-section, the y-axis vertically down- 
ward (the z-coordinate being in the longitudinal 
direction), the three components then are o. , 
o,, and 7,z,, and the stress ellipsoid is reduced 
to a stress ellipse in the xy-plane (Fig. 3). 
The stress ellipse provides a simple and 
convenient means of illustrating some of the 
fundamental differences between the stress- 
tensor and the force-vector. Two of them con- 
cern the methods of resolution and composition 
of forces and stresses. To illustrate the first 
case, assume a surface force acting in the 
direction of o, (Fig. 3). Its vector component 
at right angle, in the y-direction, is zero, 
whereas the corresponding stress component 
o, may be of any magnitude. This again shows 
that the process of resolution of forces for 
the purpose of determining the stress com- 
ponents acting upon a surface in any direction 
other than that specified in the limited defini- 
tion of stress is inapplicable. A similar situa- 
tion exists with respect to the composition of 
forces and stresses. When several forces act 
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FicurE 3.—Two-DIMENSIONAL STRESS Com- 
PONENTS AND STRESS ELLIPSE 


phrased statements. Such deductions are bas- 
ically wrong and lead often to complete cen- 
fusion in the process of reasoning. When a 
solid is subjected to several simultaneous 
stresses, they can, of course, also be combined 
into a single stress system, but the procedure 
is quite different from the method of vector 
addition. The same is true for the respective 
results. The combination of stresses is based 
on the principle of superposition. It states that 
addition—in the form of linear functions—of 
the stress components from two, or several, 
valid stress systems, produces a new one which 
is likewise valid, provided the resulting defor- 
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mation is still within the elastic limit. An 
important example of this principle will be 
given later. Again the stress ellipse is useful 
to illustrate the difference between the meth- 
ods of combining stresses and forces. Let us 
assume a uniform and uniformly directed pres- 
sure o, acting in the direction of the x-axis 
and a second similar pressure o, acting in the 
y-direction. Each of these systems satisfies 
the required equations and, hence, represents 
a permissible stress system. The method of 
superposition yields a two-dimensional stress 
system characterized by the stress ellipse shown 
in Figure 4. Its principal directions coincide 
with the co-ordinate axes, and the intensity 
of the stress acting on any surface element 
at right angles to the xy-plane is represented 
by one of the radius vectors of the stress 
ellipse. It is to be noted, however, that the 
direction of the normal to the surface element 
and that of the corresponding radius vector 
are not identical (except for the principal 
planes); the exact relation between these direc- 
tions is somewhat complicated but need not 
be considered here. If, on the other hand, we 
should apply the method of vector-addition 
to forces of equal intensity as o, and o,, the 
result would be totally different. Compare 
Figures 4a and 4b. The vector sum of f, and 
f, yields a force f, which is substantially in 
excess of any radius vector in Figure 4a, re- 
gardless of direction. Hence this procedure 
cannot be applied to stress relations. It is, 
nevertheless, frequently employed in discus- 
sions on structural problems, and leads to 
serious errors, both in the quantitative results 
as well as the more general qualitative deduc- 
tions. 

A further difference between stress and force 
concerns statements of direction. One may 
speak of a force as acting in a certain direction, 
say towards the north. When applied to internal 
stresses, such statements have no sense. A 
stress component acting upon one side of a 
surface element exists only in co-ordination 
with a component of equal intensity but op- 
posite direction, acting on the other side. This 
is true for normal, as well as shear stresses. 
Hence a pressure or a tension may exist in a 
north-south direction, but not towards the 
north or towards the south. 
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The co-existence of stress and counter-stress 
on opposite sides of a surface element is a 
general property existing on any imaginary 
internal surface of a body, as well as on all 
its boundary surfaces. One very important 
implication of this condition concerns its ap- 
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(c) Superposition of stresses (b) Vector addition of forces 


FicurE 4.—CoMPOSITION OF STRESSES AND 
Forces 


plication to the surface of the earth. Since the 
air is incapable of sustaining a shearing stress, 
Tzy must be zero everywhere along this surface. 
Furthermore, the vertical pressure is reduced 
to the normal atmospheric pressure and cannot 
exceed it. This pressure is a very small quantity 
when compared to the internal stresses in the 
crust in which we are here interested; hence, it 
can safely be neglected. These considerations 
impose two specific boundary conditions which 
must be satisfied by every stress distribution 
applicable to geological problems: that o, and 
tzy vanish for every point along the free sur- 
face of the earth. 


METHOD OF CALCULATING STRESS SYSTEMS 


The methods applied to derive complete 
solutions of geologically significant stress dis- 
tributions will be briefly outlined. Only two- 
dimensional cases will be treated. This simplifi- 
cation is fully justified for numerous structural 
problems but is not applicable, of course, to 
those of essentially three-dimensional charac- 
ter. A further assumption is that the stressed 
medium is homogeneous. This second restric- 
tion is geologically much more undesirable; 
however, the mathematical treatment of non- 
homogeneous cases would be extremely compli- 
cated. 

The internal state of stress in a block of the 
earth’s crust is adequately portrayed by the 
following sets of curves: 
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(1) The stress trajectories, defined by the 
property that the directions of the principal 
stresses are the tangents to these curves at all 
points. One set determines the direction of the 
maximum principal stress and a second one 
that of the minimum principal stress. The two 
sets are orthogonal everywhere. 

(2) Lines of equal maximum shearing stress, 
consisting of a series of lines each of which 
connects all points of equal magnitude of this 
quantity. 

The intensities of the two principal stresses 
can be illustrated by similar curves; however, 
they are less important for our purposes and 
are omitted from the illustrations to avoid 
overcrowding. Conversely, the external forces, 
or boundary stresses, form an integral part 
of each analysis and are, therefore, plainly 
marked on the diagrams. 

The method here employed for the calcula- 
tion of the internal stress distributions is based 
on the use of the Airy stress function, which 
offers an especially convenient approach for 
two-dimensional cases. It can be shown that 
for plane stress distributions—as well as for 
plane deformations—all three stress compo- 
nents can be expressed as second partial deriva- 
tives of a single function ® which is a scalar 
function of x and y. Thus its role is somewhat 
similar to the potential function as applied to 
gravity, magnetism or other force fields. As is 
well known, the three force-components of such 
fields can likewise be expressed as partial deriva- 
tives of a single scalar function—the potential 
function—which greatly facilitates the mathe- 
matical treatment of the respective fields. How- 
ever, the vector-components of the force fields 
are the first partial derivatives of the poten- 
tial function, whereas the three components 
of the stress-tensor are the second partial de- 
rivatives of the stress function. 

Now the theory of elasticity shows that 
the equations of static equilibrium are auto- 
matically satisfied if we define the stress com- 
ponents by the expressions 
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when body forces are absent, and by the ex- 
pressions 
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when gravity represents the only body force, 
Here p is the density of the material, assumed 
to be constant, and g is the constant of gravi- 
tational acceleration. (When applied to en- 
gineering problems, expressions (1)’ are usually 
stated in a different form, the term —pgy here 
added to o, being replaced by a term —pgz, 
attached to 72, . (See Timoshenko, 1934, p. 25), 
The variation here adopted likewise satisfies 
the equilibrium equations and is preferable 
for geological investigations because the physi- 
cal interpretation of the gravity term is more 
directly apparent. The resulting loss of sym- 
metry is of no consequence.) 

In order to yield a valid stress system, an- 
other requirement must be fulfilled. This arises 
from the fact that the stress components are 
connected by Hooke’s law with the strain 
components (elongations), which, in turn, are 
defined as functions of the displacements. In 
two-dimensional problems, there are three stress 
and three corresponding strain components but 
only two displacements; hence, the former are 
not independent but must be connected by an 
equation establishing a condition of internal 
consistency of the displacements. This equation 
is called the compatibility equation. When 
applied to the stress function, it can be shown 
that a necessary—but also sufficient—require- 
ment is that this function satisfy the equation; 

O*p Oh Che 
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If this requirement is fulfilled, the stress func- 
tion provides a method for the direct determi- 
nation of the three stress components, making 
it unnecessary to first determine the displace- 
ments. The method then consists in finding 
suitable functions ® which satisfy equation (2) 
and simultaneously yield a system of boundary 
stresses which are of interest to geological 
investigations. 

Having found a suitable function of this 
type, the solution of the complete stress sys- 
tem is simple and straightforward. The three 
stress components are given directly by differ- 
entiation of &, according to formulae (1). The 
remaining significant stress quantities are then 
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obtained by the well-known equations of elastic- 
ity. The maximum shearing stress at any 


— ¢y\* 4 
me f(252) ol 


The intensities of the two principal stresses 
are: 


point is: 


0; OF dg = 4 (G2 + oy)  Tmas; (4) 


and the directions of the stress trajectories 
are given by: 


Bes 


tan 28 = " (5) 


Cz — Oy 
where 6 represents the angle between the 
outward-directed normal of the surface element 
and the positive x-axis, in the direction from 
the latter towards the positive y-axis. In the 
following calculations, the two normal stress 
components oz and a, are taken positive if 
they represent tensile stresses, negative if they 
are compressive stresses. 

The relation between the internal stress 
distribution of a body of limited extent and 
the forces causing it need some further explana- 
tion. The stress components defined by equa- 
tions (1)’ include the effect of gravity, whereas 
those of equations (1) are due only to the 
surface forces acting upon the boundaries of a 
block. In contrast to many engineering prob- 
lems where gravity is a negligible factor, it 
constitutes one of the primary forces in geo- 
logical systems and hence requires special con- 
sideration. This will be discussed in more 
detail in the next section. First we will focus 
our attention on the surface forces. The stress 
systems to be analyzed are defined by analytical 
functions which, in most cases, extend to in- 
finity. In geological problems, on the other 
hand, we wish to investigate finite portions 
of the earth’s crust which are exposed to the 
action of surface forces along their boundaries. 
Since such blocks are essentially stationary, 
the surface forces cannot be chosen arbitrarily 
but only such combinations are admissible 
which, on the whole, satisfy the equations of 
static equilibrium. Now the relation between 
the internal stress distribution defined by equa- 
tions (1) and the boundary forces responsible 
for it is directly given by a fundamental prin- 


ciple. It states that any imaginary surface 
inside a stressed medium can be made a bound- 
ary surface without causing a disturbance in 
the internal stress distribution, provided a 
system of surface forces equivalent to the 
former internal stress components is applied 
to this boundary. By this process we can 
arbitrarily select any desired location and lim- 
its of a block within the confines of an analyti- 
cally defined stress distribution and can also 
immediately state the system of surface forces 
which produces it. This principle will be used 
extensively in the later treatment of specific 
examples. 

The method of determining stress systems 
here outlined is, in some respects, an indirect 
approach. The stress distribution is defined 
first and the causal boundary forces are then 
fixed automatically. This is obviously not an 
altogether ideal procedure. Yet it has decided 
advantages and in reality is much more flexible 
than one might expect it to be. First the system 
of boundary forces thus derived automatically 
satisfies the requirement of static equilibrium. 
This imposition might constitute a tedious 
complication in all but the simplest cases if 
the direct approach of starting with arbitrarily 
selected boundary conditions were used. Sec- 
ond, numerous boundary systems represent- 
ing some of the geologically most interesting 
cases can easily be obtained by suitable choices 
in the selection of the Airy stress function. 
Another important factor greatly enhances the 
flexibility of this method. As seen from equa- 
tion (2), the stress function ® is obtained by a 
solution of a fourth order differential equation. 
Hence it contains four independent integration 
constants which will also enter the expressions 
for the stress components (1). Now, as pre- 
viously pointed out, any stress distribution 
applicable to our problems must satisfy the 
boundary conditions at the surface of the 
earth. Due to this restriction, two of the four 
integration constants must be expended for 
this purpose. The remaining two, however, 
are unassigned by general theoretical require- 
ments and, therefore, can be used for the 
modification of the boundary stresses. It will 
be seen later that this freedom provides the 
means for a considerable degree of variability 
in the selection of geologically interesting cases. 
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Tue STANDARD STATE 


E. M. Anderson (1942) has introduced the 
term “standard state,” defined (p. 137) as 


“‘a condition of pressure which is the same in all 
directions at any point, and equal to that which 
would be caused by the weight of the superincum- 
bent material, across a horizontal plane, at the 
particular level in the rock. It is assumed in this 
definition that the surface is flat, and that the strata 
are of uniform specific gravity.” 


It will be shown that the standard state—ac- 
cording to this definition—represents a stress 
system composed of two parts: (1) the effect of 
gravity, and (2) a superposed horizontal stress 
which is constant in any horizontal plane but 
increasing uniformly with depth. The manner 
in which gravity affects our investigations must 
be clearly understood, since it represents one 
of the major causes of stress. Being a body 
force, it plays a somewhat different role in the 
analytical treatment than do the surface forces. 
For this reason—and also to illustrate the 
discussed calculation methods—the standard 
state is briefly derived as follows: 

We select the Airy stress function in the 
form of a polynomial of third degree 


© = his? + bea? + karyy + be’, 


which satisfies equation (2). The stress com- 
ponents (1)’ then are: 


a 
o; = oy = 2ksx + Okay; 


Coy = — pgy = 6kix + 2key — pgy 
° Ox ; < 
Tey axa —Lkyt — ay. 


In view of the boundary conditions at the 
surface, k,; and ke must be set zero; kg and ky 
can be chosen arbitrarily. We select ks = 0; 
kg = —42pg and obtain: 


Oy = Oy = —psy; (6) 


Tzy = 


These are the stress components for the stand- 
ard state. The normal stresses are equal in all 
directions and for all points on a horizontal 
plane (y = constant) and the shear stress is 
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zero throughout the body. This stress system 
is, therefore, equivalent to the hydrostatic 
pressure in a liquid. Its physical significance 
for solid bodies is apparent from the following 
consideration. In a homogeneous solid body 
of infinite horizontal extent, uniform lateral 
extension is prevented by lack of space. There- 
fore, if the weight of the body is the only source 
of stress, the strain components (elongations) 
in the horizontal directions are zero. From this 
the relation 





Thor —* vert 
is derived, where oir and Or: represent the 
horizontal and vertical stress components, re- 
spectively, and y is Poisson’s ratio. It is seen 
that Ghor = Geers Only if vy = 0.5. Hence, the 
hydrostatic state is produced by gravity alone 
only in liquids. In solid bodies the horizontal 
stress components due to weight are substan- 
tially smaller than the vertical component; 
they are only a third of the latter if Poisson’s 
ratio is taken to be 0.25, which is a good ap- 
proximation in most cases. The term 6k in 
the above derivation of the standard state 
thus includes two effects: roughly a third of it 
is due to gravity, while the remaining portion 
represents a superposed horizontal stress, con- 
stant in any horizontal plane, but directly 
proportional to depth. 

The important result derived from these 
considerations is that the éffect of gravity can 
be incorporated into a stress system which is 
essentially hydrostatic in character. This is a 
consequence of the fact that a stress system of 
the form: 


Cy = Ty = 0 


oz = ky; 


is likewise a valid one—as shown by the above 
derivation—and thus can be combined with 
the weight of the body to produce a hydrostatic 
system. It does not matter whether the com- 
plementary horizontal stress ky actually exists 
in any particular case; the standard state 
represents an idealization which is probably 
only seldom realized in nature. Its principal 
usefulness—aside from offering a convenient 
“standard of reference—” arises from an appli- 
cation of the already discussed principle of 
superposition of stresses, which can now be 
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formulated more precisely. Let ©; and 2 be 
two stress functions, each of which represents 
a correct solution. Then the sum of the two is 


G=4+4 (7) 


also a valid solution. Now let ®, be the standard 
state which includes the effect of gravity. 
Since no other body force is involved, equation 
(7) states that we may superimpose upon it 
any other valid stress system #2 caused by 
surface forces alone. Anderson has named such 
systems “supplementary stresses.” Their math- 
ematical treatment is now divorced from the 
effects of gravity, and the stress components 
are given by the simpler expressions of equa- 
tions (1). The physical implication of this 
procedure is, of course, merely a statement of 
the fact that the standard state contributes 
nothing to the shearing stress and has no 
influence on the configuration of the stress 
trajectories. Conversely, it forms a major com- 
ponent of the total stress and, therefore, must 
be taken into account in phenomena depending 
upon the magnitude of the confining pressure. 
One of these phenomena is faulting. 


FAULTING ASSOCIATED WITH 
Stress DISTRIBUTION 


While the determination of the stress proper- 
ties is amenable to exact mathematical analy- 
sis, the occurrence of faulting presents an 
intricate problem for which no complete theory 
has yet been worked out. It is, therefore, im- 
portant to discuss its critical aspects and to 
point out the shortcomings of the present 
treatment. 

The principal object of the fault analysis is 
two-fold: (1) to know the attitude of the fault 
surfaces and their variations in the calculated 
stress systems, and (2) todetermine the location 
of faulting in the stressed medium, i.e., the 
distribution of the stable and unstable portions. 

The first problem is readily solved if we 
accept a widely held theory on the relation 
between the fault surfaces and the directions 
of the principal stresses. It states that fracture 
occurs along two planes which have a specific 
angle @ with the direction of the maximum 
principal pressure. This angle is a material 
constant and, therefore, variable. It is, how- 


ever, always less than 45°, and for most types 
of rock falls between 20° and 40°. A value of 
30° appears to form a good overall approxima- 
tion. In three-dimensional stress systems the 
direction of the two fault surfaces is parallel 
to that of the intermediate principal stress, 
the angles 6 being measured in the plane con- 
taining the maximum and least principal 
stresses. In the two-dimensional cases here 
analyzed, we shall presume that this is the 
plane containing the variable stress compo- 
nents—that the third dimension is parallel 
to the intermediate principal stress. Then the 
attitude of the fault surfaces is obtained by 
drawing two sets of curves having everywhere 
a selected constant angle with the direction of 
maximum principal pressure. This is evidently 
a simple procedure, once the configuration of 
the stress trajectories has been completely de- 
termined. While the choice of @ is arbitrary 
within certain limits, a value in the neighbor- 
hood of 30° appears most appropriate. Cor- 
responding to the plus and minus sign of 
6, the method always yields two sets of po- 
tential fault surfaces, which intersect each 
other at an oblique angle. While there is no 
theoretical basis for a discrimination between 
the two systems in a homogeneous medium, 
geological observations indicate that in most 
areas only one of them has been utilized, or 
that it strongly predominates. Nevertheless, 
examples of the co-existence of both sets within 
the same structural belt are known for all 
classes of faults. 

The validity of the above theory on the 
attitude of the fault surfaces has been ques- 
tioned by some authors (Jeffreys, 1936), but 
it finds much support both from the results 
of experiments and from geological observa- 
tions. It is based on the assumption that frac- 
ture in a stressed medium takes place along 
that plane which offers the least resistance to 
it. To produce faulting, the available energy 
represented by the condition of strain has to 
overcome not only the strength of the material 
but also the internal friction opposing the 
motion. The latter is proportional to the normal 
pressure across the fault surface. The shear 
stress is always greatest in a direction at 45° 
to those of the maximum and minimum prin- 
cipal stresses. Now if we rotate a plane from 
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this position towards the axis of maximum 
principal pressure, both the shear stress—which 
tends to produce faulting—and the normal 
pressure—which opposes motion—become 
smaller. Depending upon the specific proper- 
ties of the material, there will be one par- 
ticular angle for which the combined effect 
of the two opposing factors becomes most 
favorable for the occurrence of faulting; this 
is the angle @ above referred to. A detailed 
mathematical formulation of this principle is 

found in Anderson (1942, p. 9-10). 

The second problem of faulting accorded 
attention here deals with the determination of 
the boundaries between the stable and un- 
stable portions of the stress systems. We will 
first examine the problem solely upon the 
basis of the statical theory as applied to the 
original stress distribution. 

The basic principle is easily stated: faulting 
occurs where the shearing stress exceeds the 
strength of the material. The magnitude of 
strength depends, first, upon the material, 
falling in the range of approximately 500- 
3000 kg/cm? under atmospheric pressure for 
the most common rock types. Second, it varies 
widely for a given material as a function of 
other factors, such as the confining pressure, 
temperature, presence or absence of solutions, 
and even the time of application of the shear- 
ing stress. Only the first of these has been 
adequately investigated in laboratory experi- 
ments. It was generally found that the breaking 
pressure increases with the confining pressure 
and that in many cases the rate of increase 
is nearly linear over a wide pressure range. 
Since little information is available on the 
influence of the other factors, the boundary 
equation for the separation of the stable and 
unstable portions adopted in the present paper 
is: 

min S Nomaz + 0 (8) 
where: 

Omin = Maximum compressive stress (compressive 
stresses being given the negative sign, the 
maximum compressive stress is, algebraic- 
ally, the smallest), 

maz = minimum compressive stress, 

o) = breaking strength under atmospheric con- 
fining pressure, 

nm = material constant. 
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As pointed cut, while the standard state 
contributes nothing to the shearing stress in a 
homogeneous medium, it adds a vital com- 
ponent to the total confining pressure; hence 
min @Nd gomez in (8) represent total principal 
stress, including that due to the standard 
state. They are obtained by adding to the 
components from equation (4) the hydrostatic 
pressure defined by (6): 


onin = 4(oz + oy) — Tmaz — _ (9) 


Omaz = 4(¢2 + oy) + Tmar — psy 
For the calculation of numerical examples it 
will be necessary to assign the constants pg, 
oo and m specific values. The selections for the 
present paper are: 


pg = 250 kg/cm? per kilometer or 400 kg/cm? per 
mile, 

oa = 1000 kg/cm?, 

n= 4, 


While in nature each of these quantities is 
variable, the figures chosen represent the right 
order of magnitude. The value for m is based 
mostly on geological evidence, but it is also 
substantially in agreement with the results of 
experimental data. 

In the practical application, a simple method 
was used for the solution of equation (8) which 
is illustrated in Figure 5. The curve o,; repre- 
sents the standard state, to which the prin- 
cipal supplementary stresses derived from 
equation (4) are added algebraically (compres- 
sive components being negative). Then the 
curve: (ndmez + 0) is drawn (dashed line). 
Faulting occurs where the latter is to the left 
of omin. The intersections between the two 
curves mark the boundaries between the areas 
of stability and those in which the shear stress 
exceeds the strength. In the specific examples, 
the boundary lines were determined for not 
only one, but several, values of the significant 
constants appearing in the formulae for the 
external boundary stresses. 

In the preceding discussion, the distribution 
and mode of faulting were examined solely in 
relation to the original stress system. With 
the first occurrence of faulting, an entirely new 
problem arises. The process of fracturing causes 
an instantaneous and drastic alteration in the 
pre-existing stress system, accompanied by the 
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introduction of a new internal boundary sur- 
face. It is thus evident that any subsequent 
faulting must take place in accordance with 
the changed conditions, and the value of the 
original stress system as a guide to further 
deformation may rightfully be questioned. This 
interaction between repeated changes in the 
stress distribution and the successive occur- 
rences of faulting introduces a complication 
which presents great difficulties to an ade- 
quate mathematical analysis. No attempt is 
made here to deal with this problem. However, 
some qualitative arguments are advanced which 
may shed some light on the degree of control 
which the original stress distribution still exer- 
cises during the more advanced stages of fault- 
ing. 

A mathematical solution of the alteration of 
stress due to the occurrence of a fault is given 
by Anderson for a simple case (1942, p. 144- 
156); this is the only example known to the 
writer. The results demonstrate that pro- 
nounced changes in both direction and in- 
tensity of the stresses occur in the immediate 
neighborhood of the fault surface, but that 
they diminish rapidly with distance. This is 
especially true in a direction at right angle to 
the fault where the alterations soon become 
insignificant. Near the central portion of the 
fracture the shear stress is greatly reduced, 
which has the effect of preventing further 
faulting in this vicinity. Conversely, there is a 
substantial increase in stress at both ends of 
the fracture, indicating a tendency for lateral 
extension, once the process has been initiated. 
While the solution applies to a specific case 
only, it is probable that the basic results are 
similar in other, more complicated situations. 

In the light of these data we may visualize 
the following sequence of events. Under the 
action of external forces, a condition of stress 
is produced in a block of the earth’s crust. 
As long as the shearing stress is everywhere be- 
low the strength of the material, the stress 
distribution is static and can be determined 
by means of the discussed methods. The vari- 
ous examples given below are rigorously cor- 
tect during this stage only. With further in- 
tensification of the stresses, a point is reached 
where the strength is exceeded at some par- 
ticular place, and faulting will occur there. 


The general location of this place is determi- 
nable from the initia] conditions. How far the 
fracture extends itself laterally into regions 
where the original shear stress is still well 
below the strength is an unsolved problem. 
Theoretical results suggest the possibility of 
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pronounced extension, whereas geological ob- 
servations demonstrate that the majority of 
faults are of rather limited extent. Regardless 
of this apparent conflict, the principal result 
will be a drastic diminution of stress in the 
neighborhood of the origin of the fracture, 
while farther away the initial stress distribu- 
tion remains substantially unaltered. The first 
fault thus provides local relief in the most 
intensively stressed portion of the block. After 
its occurrence, the further building-up of the 
stress system is presumed to go on, and addi- 
tional portions of the block will gradually reach 
and exceed the strength limit. Renewed fault- 
ing then occurs some distance away from the 
first fracture, at a place where the stress-relief 
is of negligible amount. Hence, while no longer 
exactly the same, the conditions producing the 
new fracture should still conform closely to 
the original stresses. A further implication of 
this process of reasoning is the probability that 
the later faulting belongs to the same set as 
the initial one, since a secondary fracture utiliz- 
ing the alternative complementary direction 
would lead straight into the protective zone 
of the preceding one. This, perhaps, is the 
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explanation of the relative scarcity of com- 
plementary faults. 

In this manner the process may continue 
until the maximum development of the system 
has been reached. At this stage the extent of 
the unstable areas has likewise reached its 
maximum expansion which, depending upon 
the ultimate intensities of the causal forces, 
will include a smaller or greater portion of the 
entire block. Due to the limited extent of the 
influence of previous fractures in adjacent areas, 
it appears probable that the initial faulting 
in each portion of the unstable segments is 
primarily controlled by the original stresses. 
The local relief of stress provided by the exist- 
ing fractures will then gradually disappear 
again due to the tendency of the original stress 
system to re-establish itself. This can lead 
to renewed faulting in already fractured seg- 
ments. The process thus becomes more and 
more complicated, on account of the intro- 
duction of a new internal boundary surface 
along each additional fault and the adjacent 
zone of pronounced stress alterations. In ex- 
treme cases, the stress system may be sub- 
stantially changed during the final stages of 
deformation and the further course of events 
is then beyond the scope of theoretical analysis. 

From this line of reasoning it appears justified 
to assume that the original stress system consti- 
tutes the controlling influence on the location 
and type of faulting at least during the initial 
stages of deformation in each faulted segment. 
The author, therefore, believes that illustra- 
tions showing the attitude of fault surfaces 
and the boundaries between stable and unstable 
segments based purely on the original stress 
system serve a useful purpose in studies of this 
nature. 


EXAMPLES OF Two-DIMENSIONAL 
Stress SYSTEMS 
Introduction 
Using the methods outlined in the preceding 
sections, we will now present a series of ex- 
amples of two-dimensional stress systems. One 
or several numerical examples for each group 
are illustrated graphically in Figures 6-8 and 
on Plate 1. In each figure, the upper diagram 
shows the complete solution of the internal 
stress distribution in the form of stress tra- 
jectories and lines of equal maximum shearing 
stress; the lower diagram shows the attitude 
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of the fault surfaces and the distribution of the 
stable and unstable portions of the block. All 
graphs show the respective systems of boundary 
stresses. 


Supplementary Horizontal Stress without 
Superposed Vertical Stress 


In the first group to be analyzed, we assume 
the presence of a supplementary horizontal 
stress but absence of an associated vertical 
stress component, i.e., there is no pressure or 
tension across any horizontal plane in the body 
in addition to the normal hydrostatic com- 
ponent. Mathematically this is expressed by: 


ap 
oy = eg 0 for all values of y. 


Integrating we obtain the stress function: 

© = of(y)x + ax + bf2(y) + d. 
To satisfy equation (2) 
ap a ab “ - 
— +2——_ + — = y+ = 
bat asta t Oy cxfy" (y) + bf2""(y) = 0 
the fourth order derivatives of f, and fe must 
be zero. Hence the second order derivatives 
may be either linear functions of y, constants, 
or zero. The stress components then are (equa- 
tions (1)) 
arp 
dy? 
ah 
 Oxt 


fi™(y)x + bf2"(y); 


o; = 
0; 


oes 
= ——— = I 
Try axdy fi (y). } 
The boundary conditions at the surface re- 
quire that fj(y) = 0 for y = 0. Keeping within 
the limits of the above restrictions we can set 
up the following three subgroups: 


a) fio) =0; fi) =y9+¢4 
Ox = by + d; 


oy = 0; 





Try = 0. 
filo) =9; fry) =0 
Oz = CX; 
oy = 0; 
Try = —CY. 
fig) = 4%; fz) =0 
oc = CXY; 
oy = 0; 
c 
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EXAMPLES OF TWO-DIMENSIONAL STRESS SYSTEMS 
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INTERNAL AND BOUNDARY STRESSES: LEGENO 
o, *ce* {Cf,(y) - Df,(y)} ——— = Trojectories of moximum principal pressure = Omin 
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¢ from 1000 kg/cm? to 3000 kg /cm* 


Ficure 8. SUPERPOSED HorizoNTAL STRESS DECREASING EXPONENTIALLY IN HorIzONTAL DIRECTION 


In these equations, b, c, and d are arbitrary 
constants and may be assigned any value, in- 
cluding zero. Furthermore, any linear com- 


bination of equations (10a) to (10c) represents 
another allowable system. 
The simplest case of a supplementary stress 
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system is given by (10a) if we let 6 equal zero. 
Then the superposed stress is restricted to a 
constant horizontal component o, = d, the 
other two components go, and 7,, being zero 
everywhere. From (3) and (5) it follows that 
the maximum shearing stress is also constant— 
equal to half the horizontal stress—and that 
the stress trajectories are everywhere hori- 
zontal and vertical, respectively. The associated 
fault planes will be inclined at angles of about 
+30° to the horizontal (thrust faults). 

We will next examine the stress systems of 
the first two sub-groups. The constant d is 
now selected to be zero since its contribution 
has already been discussed. The combination 
of (10a), (105), and the standard state gives 
the total stress components: 


Ts = cx + by —ay, 
ge — 
Try = —Cy; 


where a@ = pg. The term -ay represents the 
hydrostatic pressure and is the only one ap- 
pearing in the vertical stress component. The 
horizontal component oe: contains two addi- 
tional terms, the first a linear function of x only 
and the second a linear function of y only. Hence 
o; has constant gradients in both the horizontal 
and vertical directions. The shear stress 1:, is 
seen to be constant in any horizontal plane; 
further, it increases vertically at a constant 
rate which is equal to the lateral gradient of 
oz. The term ky in o, appeared previously as 
a component of the standard state; in the 
present formulae, it denotes the difference (ex- 
cess or defect) between the actual magnitude 
of this stress and that portion of it which is 
absorbed in the standard state. 

Figure 6! illustrates the case where the supple- 
mentary horizontal pressure is constant with 
depth (b = 0) and Figure 7 where it has a 
vertical gradient equal to the horizontal one 
(6 = c). The external boundary stresses are 
shown on the end-surfaces and along the bottom 
of a rectangular block. On the two sides, the 
horizontal stress o, is broken down into its 
individual components. It is seen that addition 
of the term by does not change the equilibrium 
of the block in the lateral direction, a fact which 


1 Similar to Figure 25 in the companion paper by 
M. King Hubbert (1951). 
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explains that the shearing stress along the 
bottom is the same in both cases. 

The trajectories of maximum principal pres- 
sure are curved lines, dipping downward away 
from the area of maximum compression. The 
curvature is stronger if the vertical gradient 
of a, is small. Since the trajectories are curved, 
so are the potential fault surfaces. The latter 
obviously belong to the class of thrust faults. 
The set dipping towards the area of maximum 
pressure is slightly concave upwards, the com- 
plementary set concave downwards. Thrust 
faults of the former type are very common in 
nature and the theoretically deduced curvature 
is frequently observed. The latter type appears 
to occur only rarely and little is known regard- 
ing the curvature. 

The lines of equal maximum shearing stress 
are expressed in terms of multiples of the con- 
stant c. The shearing stress naturally increases 
towards the area of greater horizontal pressure 
but it also increases with depth, the latter effect 
being more pronounced if o: has a vertical 
gradient. To determine the boundary of the 
stable portion of the block, the constant c has 
to be fixed numerically. This is done best by 
selecting certain values of the ratio c/a, thus 
expressing the lateral gradient of the super- 
posed horizontal stress in terms of a fraction 
of the vertical pressure gradient due to weight. 
In nature this ratio is evidently widely variable 
with place and time and, in fact, may assume 
any value from zero up to a magnitude large 
enough to cause thrusting. The most useful 
procedure, therefore, consists of assuming 
several ratios and calculating their respective 
boundary lines. In both Figures 6 and 7, two 
such lines are shown corresponding to the ratios 
c =aandc = a/2. 

These boundary lines are seen to be dipping 
towards the area of greater horizontal pressure. 
The rate of dip is primarily a function of the 
ratio c/a, being greater the larger this fraction 
is. It is steep only if the latter is near unity, 
when the horizontal pressure gradient ap- 
proaches the magnitude of the vertical pressure 
gradient due to weight. For smaller values of 
c, the boundary lines are nearly straight and 
their inclination diminishes rapidly. Numerical 
examples on the dip 8 are: for ¢c = a/2, 8 is 
about 10° in Figure 6 and 12}° in Figure 7, 
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fo. . = @/10, 6 is less than 2° in the former 
case. 

These results suggest the following con- 
clusions. If the supplementary horizontal pres- 
sure has only a small lateral gradient, say of 
the order of half the vertical pressure gradient 
or less, the potential area of thrusting is con- 
fined to a shallow, gently dipping wedge. The 
resulting deformation will consist of a series of 
slice-thrusts covering a broad belt but extend- 
ing only to shallow depth. The presence of a 
vertical gradient slightly steepens the wedge. 
However, the boundary becomes very steep or 
nearly vertical only if the horizontal gradient 
approaches the magnitude of the vertical in- 
crease of the hydrostatic pressure. In that 
event, thrusting can take place throughout a 
thick, but probably only narrow, zone of the 
crust. Such a condition appears to have oc- 
curred, for instance, in the marginal belts of 
several ranges in the Rocky Mountain province. 

The boundary stresses of the third sub-group 
(formulae 10c) are also characterized by the 
fact that the shearing stress is constant in all 
horizontal planes. Its vertical gradient, how- 
ever, is now a function of the second order of y, 
instead of the first order. The shearing stress 
is balanced by a horizontal pressure, again 
having constant lateral and vertical gradients 
but increasing more rapidly in, the diagonal 
direction. This case has not been further 
analyzed. 

The most general expression for the stress 
systems satisfying the assumption of absence 
of a vertical stress component is given by the 
superposition of equations (10a), (105) and 
(10c), as follows: 


o, = ial 
yey | (11) 


It is seen that the stipulation ¢, = 0 is asso- 
ciated with two additional general properties 
of the internal stress system: (1) that the 
shearing stress is a function of y only, i.e., 
constant in all horizontal planes, and (2) that 
cz has linear gradients in both the horizontal 
and vertical directions. That the first two 
properties are reversible can readily be demon- 


c 
Try = —g? — 9 
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strated by deriving the stress system based on 
the assumption: 
t= Ee = — hf). 

Using the same direct integration method as 
before and satisfying equation (2) and the 
boundary conditions at the surface, this leads 
again to equations (11). Hence, the reversed 
statement: “constancy of the shearing stress in 
all horizontal planes is associated with absence 
of a vertical supplementary stress throughout 
the body,” is also true. Thus the heading of 
this Section could equally well read “Supple- 
mentary shearing stress constant in all hori- 
zontal planes.” 


Supplementary Horizontal Stress Decreasing 
Exponentially in the Horizontal Direction 


We have investigated the condition of the 
presence of a supplementary horizontal stress 
coupled with absence of a superposed vertical 
component. In this section we again take a 
supplementary horizontal stress, but drop the 
assumption that o, be zero throughout the 
block. Instead, we now impose the condition 
that o, decrease exponentially in the horizontal 
direction. 

For this purpose we take the Airy stress 
function in the form: 


& = ce*f(y). (12) 


Applying equation (2) to this function yields 
the differential equation 


fy) + 2f") + f°") = 0 
the general solution of which is 
fly) = Asiny + Boos y + Cysin y + Dy cos y. (13) 


Substituting (13) in (12) and differentiating we 
obtain the stress components (1). The com- 
pulsory boundary conditions cy = rzy = 0 for 
y = 0 can be taken care of by a proper dis- 
position of the two constants A and B. After 
eliminating these constants, the stress com- 
ponents become: 


o, = c@{C(2 cosy — ysiny) — D(sin» + ycos y)}; 
oy = ce*{Cysiny — D(siny — ycosy)}; 
Tey = ce*{—C(siny + ycos y) + Dysiny}. 
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Let: 


fey) = sin y + y cos 9; 
f(y) = sin y — ycos 9; 


fly) = 2cosy — ysin y; 
Sly) = ysin y; 


then the stress components can be written in 
the form: 


oz = ce*{Cfi(y) — Dfx(y)} = ceF,(y);) 
oy = ce*{Chily) — Dfaly)} = oo (14) 
Tzy = ce*{—Cfo(y) + Dfs(y)} = ce*Fs(y). 


The direction of the stress trajectories is 
given' by: 


2 Fay) 


F,(y) - F,(y) : as) 


tan 28 = 
The functions F; to F3 are constant in the 
horizontal direction since they are functions 
of y only. Therefore, equation (15) shows that 
all stress trajectories are parallel curves and 
the same applies to the fault surfaces under 
the assumption of constancy of angle between 
the latter and the directions of maximum prin- 
cipal pressure. 
When writing equations (14) in the form: 


o = ¢ ex, 4), 


it is evident that the constant c has the dimen- 
sion of stress. Its value can be selected arbitrar- 
ily and determines the absolute magnitude of 
the stress components, whereas the function 
€(x, y) determines their areal variations. The 
two remaining integration constants C and D 
are dimensionless and, being likewise fully in- 
dependent, permit the selection of numerous 
variations. The degree of variability is further 
enhanced by the fact that the position of the 
bottom of the block can be chosen at any de- 
sired lower limit of the y functions, the unit of 
length also being arbitrary. The range of oppor- 
tunity thus provided can easily be exploited 
by drawing graphs of the four functions f; to f 
and making visual estimates of the combina- 
tions obtainable after multiplication with vari- 
ous constants. 

Figure 8 illustrates the case where C = —1 
and D = +2. In this and all succeeding illustra- 
tions, only the supplementary boundary stresses 
are shown; the normal hydrostatic components 
are omitted to simplify the drawings. The 
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shearing stress acting on the bottom of the block 
is no longer constant but, like the superposed 
horizontal stress, increases exponentially 
towards the pressure area. The same is true 
for the vertical component which now forms a 
part of the stress system. All three boundary 
stresses diminish rapidly away from the area of 
compression. The zone of potential faulting 
consists of a narrow, nearly vertical, belt and 
a shallow outwardly protruding wedge near 
the surface. Figure 8) indicates three boundary 
lines limiting the area of stability; they cor- 
respond to values of the constant ¢ of 1000, 
2000, and 3000 kg/cm, respectively. 

The general fault system in the shallow por- 
tion of the unstable segment is similar to that 
obtained in the preceding stress systems, except 
for a slight change in attitude. The thrust 
faults dipping towards the pressure area are 
less steeply inclined. This divergence increases 
with depth where, because of a gradual clock- 
wise rotation, the fault surfaces of this set 
become practically flat and finally even over- 
turned. It appears doubtful, however, that the 
extreme stage is ever realized, because the 
sudden decrease of oz near the lower left-hand 
corner of the block does not represent a 
boundary condition which one would expect 
in nature. As already stated, we are at liberty 
to match the bottom of the block with any 
desired point of the (c: , y) curve, for instance, 
that for which o, reaches a maximum. In that 
case the lateral pressure increases throughout 
the entire thickness of the block, thus providing 
a stress system which fully satisfies our intuitive 
concept of geological probability. Figure 8 
demonstrates that even within these limitations 
the thrust planes reach a practically horizontal 
attitude (in the vicinity of approximately half 
of the total depth to which calculations were 
carried out). 

It is thus evident that stress systems of this 
type produce thrust faults which are only 
gently inclined at shallow depth and become 
nearly horizontal at greater depth. Such atti- 
tudes are frequently observed in nature and are 
usually interpreted to indicate subsequent tilt- 
ing of the fault surfaces. It has been shown here 
that the gently inclined attitudes may equally 
well be explained as original features of the 
stress system. 
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Variable Vertical and Shearing Stress 
Acting on the Bottom of a Block 


The preceding two sections were principally 
concerned with various forms of horizontal 
compression acting on the two sides of a block. 
We now focus on stresses acting on the bottom 
of a block. These may be superposed vertical 
stresses, shearing stresses, or a combination of 
both. It is evident that a vertical component 
which is constant laterally is of no interest for 
geological purposes. Such a condition may 
exist within a limited width of the crust, but it 
then forms only a part of a broader system in 
which these stresses eventually die out sideways. 
Hence, only a laterally variable vertical com- 
ponent will be investigated. The most suitable 
form of lateral variation is that of a sinusoidal 
curve. This is the most convenient assumption 
from the mathematical standpoint and in all 
probability also forms the most satisfactory 
approximation to many actual cases of differen- 
tial vertical uplift. 

Examples of both constant and variable 
shearing stresses along the bottom of a block 
were already encountered in association with 
superimposed horizontal pressures; however, if 
variable, they had a strictly unidirectional gra- 
dient. This property renders them inapplicable 
for many important phenomena, such as the 
drag produced along the bottom of the crust 
by the action of convection currents in the 
underlying substratum. The treatment of this 
problem likewise leads to a sinusoidal form of 
the primary boundary stresses. 

To obtain stress systems of the desired type 
we select the Airy stress function in the form 


® = sin ax-f(y), (16) 
where f is a function of y only and 


2 
a= — 


l 


Application of equation (2) to this stress func 
tion yields the following equation for f(y): 


af(y) — 2orf(y) + f(y) = 0. (17) 


The general solution of this differential equa- 
tion is (Timoshenko, 1934, p. 45) 





S(y) = Ci cosh ay + C; sinh ay (18) 
+ Czy cosh ay + Cyy sinh ay. 


The stress components (1) are obtained by 
inserting this function of y into the stress 
function (16) and subsequent differentiation. 
They are: 


oz = sin ax{Cjo? cosh ay + Co? sinh ay 
+ Cya(2 sinh ay + ay cosh ay) 
+ Cya(2 cosh ay + ay sinh ay) } 
oy = —a* sin ax{C; cosh ay + Cz sinh ay 
+ Cry cosh ay + Cy sinh ay} 
Tzy = —a cosax{Ca sinh ay + Ca cosh ay 
+ C,(cosh ay + ay sinh ay) 
+ C,(sinh ay + ay cosh ay) } 


The mandatory boundary conditions at the 
surface are satisfied if we let: 


C, = 0; C.- — 
a 


This imposition reduces the stress components 
to the form: 


oz = asinax{C;(sinh ay + ay cosh ay) 
+ C,(2 cosh ay + ay sinh ay) }; 
oy = —a sin ax{—C;(sinh ay — ay cosh ay) 
+ Ciay sinh ay}; 
Try = —a cos ax{Cyxy sinh ay 
+ C,(sinh ay + ay cosh ay)}. 


(19) 


The remaining two integration constants C3 
and Cy, are available for suitable selections of 
the boundary stresses. From the form of equa- 
tions (19) the proper choice evidently consists 
of functions of the type 


oy = —A sin az; 
Tzy = —B cos ax; (20) 


where A and B are arbitrary constants. These 
boundary stresses are placed along the bottom 
of a block, the position of which is determined 
by y = c, ¢ being likewise a constant. Sub- 
stituting functions (20) for the left-hand sides 
of the corresponding equations (19) and re- 
placing y by c, we obtain two equations from 
which the as yet unassigned integration con- 
stants Cz and C;, can be calculated in terms of 
the remaining constants a, c, A and B, They 
can then be eliminated from equations (19). 
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Performing these steps we arrive at the follow- 
ing final expressions for the stress components: 


o, = sinax|—hifi(y) + hofely)} = sinaxF,(y); ) 
oy = sinax|—hifsly) — hefa(y)} = sin arxFa(9)5t (21) 
Try = cosax{hifs(y) — kefi(y)} = cosaxF;(y);) 


where: 


f(y) = sinh ay + ay cosh ay; 
fly) = 2 cosh av + ay sinh ay; 
fily) = sinh ay - wy cosh ay; 
f(y) = ay sinh ay; 


(22) 


and: 


Aac cosh ac — Bacsinh ac + A sinh ac ; 
? 





: sinh? ac — atc? 
(23) 
k Aac sinh ac — Bac cosh ac + B sinh ac 
2 => 





sinh? ac — a* 


Equations (21) represent the general solution 
for the internal stress components produced 
by sinusoidal boundary stresses acting on the 
bottom of a block. 

It can easily be verified that equations (21) 
satisfy the imposed boundary conditions, both 
the mandatory ones at the surface, as well as 
the selective ones at the bottom. All three stress 
components are seen to vary in the horizontal 
direction according to sinusoidal functions, and 
in the vertical direction according to the more 
complicated hyperbolic functions /; to f4 (equa- 
tions 22). The latter are multiplied by two 
constant factors, k; and k; which, in turn, are 
expressions of four other constants, a, c, A, and 


B (equation 23). Since a is defined by om, the 


terms ac, ay, and functions thereof, are dimen- 
sionless, whereas the constants A and B and, 
therefore, k; and ke, have the dimensions of 
stress (m, [-!, t-*). The latter result is, of course, 
also directly obtainable from equations (20). 
The significance of the four constants is: 


is the depth of the bottom of the block 


along which the boundary stresses are act- 
ing; 


GC @ Ve 


Qn. . , ; rao 
| = W— is the distance, in the horizontal direction, 
Qa 


of a full wave length of the sinusoidal vari- 
tions of the stress components; 

A and B represent the maximum values of the sup- 
plementary vertical and shear components 
at the bottom. 


All four constants are mutually independent 
and may be assigned any desired value. For 
instance, any combination of thickness and 
lateral extent (wave length) of the block is 
permissible. The same applies to the ratio of the 
amplitudes of the two supplementary boundary 
stresses. Letting A equal zero yields a case of 
pure drag along the bottom, whereas B equal 
zero corresponds to differential vertical uplift 
alone. Though it is questionable that these 
latter conditions exist in nature, or are even 
closely approached, it is of interest to examine 
the general behavior of these two extreme cases, 
On the whole, it is seen that equations (21) 
provide the solution for a substantial variety 
of situations. 

The form of the stress trajectories is given by: 


1 2F;(y) 
tan ax F,(y) — F2(y)* 





tan 28 = (24) 
This equation shows that, for values of x de- 
noting uneven multiples of a quarter wave 
t 3 er , 
74 etc.), the two principal direc- 
tions are vertical and horizontal, respectively, 
and for values of x denoting even multiples of a 


length (x = 


21 
quarter wave length (x = 0, 7 etc.), they are 


both inclined at 45°. For all other values of z, 
the directions of the trajectories are complicated 
functions.of y. In the horizontal direction they 
are constant for those specific depths for which 
either F; or F, — Fo are zero. The former case 
always includes the surface of the block, where 
the trajectories are horizontal and _ vertical; 
the latter, where the trajectories are at 45°, 
is only occasionally encountered. Of special 
importance in the construction of the stress 
diagrams are those points for which the right- 
hand side of equation (24) is indeterminate 
(for instance, because of the condition of zero 
over zero). These are singular points for which 
the directions of the principal stresses are un- 
defined. They occur at places where the shear- 
ing stress is zero. If located internally they are 
associated with a 90-degree turn of the stress 
trajectories when passing the point in a vertical 
or horizontal direction. Owing to the complexity 
of the expressions for the stress components, 
the complete solution of equation (24) and the 
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general equations (3) and (4), is best carried out 
by numerical and graphicel methods. 

“A series of specific examples is shown on 
Plate 1. The boundary stresses along the bottom 
of the block are clearly illustrated in each case. 
In all instances, a full wave length of the 
sinusoidal stress variations is shown, but the 
lateral extent of the block is left undetermined. 
It should be remembered that we are at liberty 
to terminate the blocks at any desired place 
by the simple expedient of applying the appro- 
priate boundary stresses. We might be in- 
terested, for instance, in a full wave length 
of the diagram, cr in any one of several signifi- 
cant portions, such as the positive or the nega- 
tive zone of the supplementary vertical stress. 
It will be seen later, however, that in some 
systems the associated horizontal stress com- 
ponent ¢, varies in the vertical direction in a 
manner which probably excludes it from con- 
sideration as a primary boundary stress. Its 
role should then be restricted to that of a purely 
internal stress. In such cases, the logical choice 
of the lateral limits of the block is restricted to 
the 0- and 180-degree positions (Pl. 1, B 
and D). To facilitate this aspect of the analy- 
sis, the horizontal and lateral shearing stress 
components have been added to the diagrams 
at the most important positions. From equa- 
tions (21), it is seen that ¢; is zero for all points 
* = 0 or 1/2 but reaches maximum values 
(positive or negative) for all points x = 1/4, 
31/4 etc. The relations are similar for the shear 
stress components, except for a phase change 
of one-quarter of a wave length. 

The boundary lines segregating areas of 
stability from those subject to faulting were 
again calculated for several values of the maxi- 
mum supplementary boundary stresses, the 
usual selections being 1000, 2000, and 3000 
kg/cm*. In this connection, a remark concerning 
the symmetry of the diagrams is in order. Equa- 
tions (21) show that all quantities which de- 
pend solely upon the supplementary stress 
components have the full degree of symmetry 
characteristic of the trigonometric functions. 
These include the stress trajectories, the lines 
of equal maximum shearing stress and the 
attitudes of the potential fault surfaces. Con- 
versely, the boundary curves of the stable areas 
do not display this symmetry between the 


positive and negative halves of a full wave 
length. The reason can be seen from a study 
of Figure 5. For any two points half a wave 
length apart, omin and omaz Change sign, whereas 
the line representing the standard state always 
remains negative. Hence, there is a switch in 
position between the two principal stresses 
with respect to the line representing the hydro- 
static pressure. This causes a completely differ- 
ent distribution between the stable and un- 
stable segments in the positive and negative 
half wave length. The only exception is for 
the points x = 0 and //2, for which the inter- 
change is of no consequence due to the fact 
that omez = —Omin. This lack of symmetry of 
the boundary lines represents a striking feature 
of the diagrams. It is another consequence of the 
fact that the occurrence of faulting is a function 
not only of the maximum shearing stress but 
also of the total confining pressure. Gravity 
once more exerts a decided influence upon the 
results. 

Plate 1, A represents an example of a rela- 
tively narrow wave length. The selected con- 
stants are: = 31.4 miles, c = 10 miles, and 
A = B. The wave length here chosen roughly 
equals the distance between adjacent mountain 
ranges in portions of the Basin-and-Range 
province, notably Nevada. This statement is 
made to facilitate a mental association of the 
selected width with certain well-known struc- 
tures. Although the diagram covers the two 
adjacent negative quadrants, we shall limit 
our attention to the positive half wave length. 
The supplementary stress system is composed 
of a vertical pressure gradually decreasing from 
a maximum value at the center (a = 90°) to 
zero at both ends (a = 0 and 180°), and of a 
bottom shear stress increasing from the center 
towards both ends. The two primary stresses 
set up a horizontal tensile component ¢, which 
is nearly constant with depth and has the same 
lateral distribution as o,. Only a downward- 
directed shearing stress acts on the two sides 
of the block. The necessity of its presence is 
immediately realized from a consideration of 
the equilibrium requirement in the vertical 
direction; it counteracts the upward push at 
the bottom due to the superposed vertical 
pressure. 

The stress trajectories are vertical and hori- 
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zontal in the central portion of the block, but 
gradually turn into a 45-degree position towards 
both ends, except within a thin zone immedi- 
ately below the surface. This rotation is asso- 
ciated with corresponding changes in the 
attitude of the fault surfaces. In the central 
upper portion of the block, normal faults occur, 
having the characteristic dips for this class of 
faults. Away from the center, the dip of one 
of the two sets gradually increases, finally be- 
coming vertical, or even slightly overturned, 
near the ends. The complementary set de- 
creases in dip, changing from a normal 60° at 
the center to as low as 15° in the deeper end 
portions. These changes in attitude do not 
apply to the surface, where the directions of 
the stress trajectories remain unaltered, be- 
cause of the required boundary conditions. 
There is, however, a rapid transition with 
depth from this normal position along the 
surface into the significant changes in orienta- 
tion below. 

During the gradual build-up of a stress 
system of this type (lower diagram of PI. 1, A), 
instability with resultant faulting is first reached 
in a broad, but relatively thin, layer below the 
surface. With continued intensification of the 
primary stresses, this zone expands with depth 
and, to a lesser degree, sideways (at the surface 
it never reaches the two endpoints, since they 
represent singular points of zero shear stress). 
When their maximum value has attained a 
magnitude of the order of 3,000 kg/cm’, it has 
spread out to embrace nearly the full depth of 
the outer portions of the block, but it does not 
yet cover the deeper central portion where the 
maximum shear stress is still below the strength. 
The process of this gradual expansion of in- 
stability is illustrated graphically by the three 
boundary lines corresponding to values of 1000, 
2000 and 3000 kg/cm, respectively, of the 
maximum value of the uplifting force at the 
bottom. 

The changes in attitude of the normal faults 
associated with this stress system are of con- 
siderable interest. Faults which are steeply 
dipping, vertical, or even slightly overturned, 
are frequently observed in areas where the 
dominating tectonic force appears to have been 
a differential vertical uplift. Numerous ex- 
amples of this type have been described, for 
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instance, from the Basin-and-Range province 
and from those portions of the Rocky Moun. 
tains in which vertical deformation predom. 
inates. If the dip is less than 90°, such faults 
are still assigned to the class of normal faults; 
if overturned, they are usually designated as 
reverse faults. Anderson has emphasized that 
faults of this nature belong in many instances 
to an entirely different class, the transcurrent 
or lateral class. This is often, but by no means 
always, true. An alternative explanation, ap- 
plicable to those cases where lateral displace- 
ment is lacking, is provided by the present 
analysis. It shows that steeply dipping normal 
faults and reverse faults can be the direct prod- 
uct of an original stress system such as that 
here analyzed. They occur on the flanks of the 
uplift, beginning at a shallow depth. To be ob- 
servable at the present surface, we merely have 
to assume that the uppermost layer has since 
been removed by erosion. The rather wide 
occurrence of such faults is not unexpected in 
view of the probability that the causal stress 
system is frequently realized in nature. How- 
ever, this system is essentially of local character, 
which explains the fact that near-vertical normal 
faults and reverse faults form no part of 
Anderson’s classification, based exclusively 
upon regional stress systems. 

The complementary set of low-angle normal 
faults also deserves some attention. Faults of 
this type appear to be rare and are seldom men- 
tioned in the literature. An outstanding case, 
however, has been described by C. R. Longwell 
from the Desert and Sheep ranges in southern 
Nevada (Longwell, 1945). Here Longwell found 
the presence of an extensive system of low- 
angle normal faults, all located on the flanks 
of a major anticline and dipping towards its 
axis. The range of the dips is mainly from 10° 
to 25°, and the form of the fault surfaces is 
further characterized by a slight upward con- 
cavity. This combination of facts agrees re- 
markably well with the general behavior of the 
low-angle normal faults emerging at the surface 
inside the positive quadrants. The correspond- 
ing dips in Plate 1, A are not quite as low as 
those described by Longwell from southern 
Nevada, but it can logically be explained by 
assuming a further flattening due to a subse- 
quent rise of the central portion of the uplift. 
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Another characteristic feature shown on the 
diagram is the reversal in curvature of the low- 
dipping normal faults with depth. Its verifica- 
tion in nature would be contingent upon finding 
an example similar to the Desert-Sheep ranges, 
but exposed at a lower level. No clear-cut 
example of a combination of the near-vertical 
and low-angle, normal faults is known to the 
author, probably because of the protective in- 
fluence of preceding fractures. The question 
of why very steep normal faults occur much 
more frequently than flat ones is less easily 
answered. 

Plate 1, B—D presents a series of examples 
with a broader wave length. While the thick- 
ness of the block is the same as in the preceding 
case (c = 10 miles), the wave length is now 
twice as large, 62.8 miles. This distance is more 
comparable to that between adjacent ranges 
in certain portions of the Rocky Mountains. 
The difference between the three cases is in 
the ratios of the superposed vertical and shear- 
stress components. They are: A = 0 in Plate 1, 
B, A = B/2 in Plate 1, C and A = 2B in 
Plate 1, D. 

The results in the second case (Pl. 1, C) are 
quite similar to those in the preceding example. 
Evidently, the doubling of the wave length is 
roughly compensated by a reduction of the 
ratio A/B by half. The only notable difference 
lies in the relative proportions of the stable and 
the unstable segments. It is seen that for equal 
values of Byoz, the depth penetration of the 
faulted portions is much greater in Plate 1, C 
than in Plate 1, A. For Braz = 3,000 kg/cm’, 
it now embraces the entire positive half of the 
wave length, excepting two small areas in the 
immediate vicinity of the corners of zero shear- 
ing stress. 

The negative half of the wave length illus- 
trates the conditions due to the presence of a 
variable vertical tensile stress acting on the 
bottom of the block, in combination with a 
shear stress component of equal size but oppo- 
site directions as compared with the positive 
quadrants. There is, of course, no absolute ten- 
sion at this depth, the supplementary tensile 
stress being always smaller than the standard 
State. It represents a partial relief of the pres- 
sure due to the weight of the overburden. The 
boundary stresses at the bottom set up a hori- 


zontal compression which may produce thrust 
faults within a thin layer below the surface. 

Major changes in the character of the stress 
system accompany alterations of the above 
ratio A/B in either sense. Although it is doubt- 
ful whether these cases still fall within the 
domain of geological reality, they demonstrate 
the nature of the changes in both directions. 
Plate 1, B represents the extreme case where 
the supplementary stress applied to the bottom 
of the block consists of a sinusoidal shear stress 
only, A being zero. Now the internal horizontal 
stress component is no longer approximately 
constant with depth, nor even of uniform sign. 
At the 90-degree position o, changes from a 
maximum compressive stress at the surface to 
a maximum tensile stress at the bottom; at 
270°, the signs are reversed. At some particular 
depth along these positions, singular points of 
zero shearing stress occur, accompanied by a 
90-degree turn in the direction of the stress 
trajectories. A second pair of singular points is 
found at the 0 and 180-degree positions. The 
latter points are further characterized by the 
fact that the horizontal plane containing them 
is, like the surface, a principal plane where all 
trajectories are vertical or horizontal. The 90- 
degree switch here takes place in the horizontal 
direction. 

A peculiar and important feature of a stress 
system of this type is that it produces two layers 
of faulting, separated by an intervening stable 
zone. Corresponding to the described variations 
in oz in the first half wave length, there is a 
zone of normal faulting in the lower part of the 
block, underlying a thin zone of thrusting in 
the uppermost portion. In the second half 
wave length the distribution is reversed, but 
thrusting in the lower portion does not occur 
unless the boundary stress at the bottom reaches 
very high values, exceeding 3,000 kg/cm? (al- 
though these areas are essentially stable, the 
attitude of a few potential fault surfaces has 
been indicated in the diagram). Conversely, 
at the 90-degree position, the strength limit 
at the bottom is reached for comparatively 
small values of Bnaz, due to the large tensile 
stress component. 

The effect of a change in ratio in the other 
direction is illustrated in Plate 1, D. Here the 
ratio A/B is equal 2, the vertical stress com- 
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ponent now being the larger one. Although the 
shear stress component is still of appreciable 
magnitude, the resulting effects are equally 
pronounced as in the preceding case. They are 
also quite similar, though of opposite sense. 
The horizontal component oe; again has a large 
vertical gradient, accompanied by a change in 
sign and by the occurrence of singular points. 
(They are absent, however, at the 0 and 180- 
degree positions.) The distribution of faulting 
is the reverse of that encountered in Plate 1, B. 
The upper portion of the first half wave length 
is a zone of normal faulting, whereas the lower 
portion is stable up to very large magnitudes 
of the boundary stresses. The second half wave 
length is roughly divided into three parts, 
namely, an upper layer of thrusting, an inter- 
mediate stable one, and a lower one of normal 
faulting. One important difference between the 
two cases—aside from their opposite char- 
acter—lies in the relative proportions of the 
stable areas. This is especially marked in the 
vicinity of the 0- and 180-degree positions, 
where almost the entire thickness of the block 
is stable in Plate 1, B and unstable in Plate 1, 
D. In the latter case, the faulting here en- 
countered is quite similar to that previously 
found in Plate 1, C, one set consisting of nearly 
vertical normal or reverse faults and the com- 
plementary set of gently dipping normal faults. 

The stress system illustrated in Plate 1, D is 
obviously due to tectonic forces causing differ- 
ential uplifts and depressions in adjacent belts 
of the crust. The former correspond to the 
positive, the latter to the negative half wave 
lengths of the diagrams. The primary vertical 
stresses applied at the bottom—representing 
an uplifting force in the positive belts and a 
partial relief of pressure in the negative zones— 
are probably accompanied by shear stresses of 
smaller but still appreciable magnitude. Pos- 
sibly the ratio selected in Plate 1, D is already 
weighted too strongly in favor of the vertical 
component to be geologically significant. Dia- 
grams B and D are nonetheless of interest due 
to the result that a departure in either direction 
from a certain “‘balanced”’ ratio produces stress 
systems in which a zone of thrusting may occur 
above a zone of normal faulting and vice versa, 
the two faulted layers being separated by a 
stable zone. 


CONVECTION CURRENTS 


Another important usefulness of the equa- 
tions derived in the preceding section (equations 
20 to 23) consists of their application to con- 
vection currents. Several authors (Griggs, 1939, 
Vening Meinesz and others, 1948), have demon- 
strated the likelihood of a periodic existence of 
large-scale convection currents in the weak sub- 
stratum underlying the crust of the earth, 
Griggs shows that of all the primary forces so 
far suggested as the basic cause for mountain 
building, only two are potentially of sufficient 
magnitude to provide acceptable working hy- 
potheses. These are (1) lateral compression due 
to thermal contraction, and (2) the viscous 
drag exerted on the bottom of the crust by 
convection currents in the substratum. He re- 
jects the first one on account of other difficulties, 
but demonstrates the adequacy of the second 
on the basis of theoretical considerations, ex- 
perimental data, and the results obtained from 
an ingenious scale model with dimensionally 
correct physical properties. In view of the 
results of these and related investigations, the 
effect of convection currents on the deformation 
of the crust has become of considerable geo- 
logical interest. 

The shearing stress exerted on the bottom 
of the crust by the viscous flow in the sub- 
stratum is evidently best expressed in the form 
of a sinusoidal curve, the half wave length 
of which represents the width of the convection 
cell. Hence, it conforms to the second of the 
two primary boundary stresses analyzed in the 
preceding section (equations 20). Furthermore, 
it appears logical to assume that the region of 
the crust above the rising current endures 4 
supplementary upward pressure, and the por- 
tion above the sinking column experiences @ 
downward drag causing a reduction in the hy- 
drostatic pressure. Therefore, the first of the 
two boundary stresses of equations (20) is 
probably also effective. This relationship be- 
tween the location of the convection cells and 
the boundary stresses acting on the bottom 
of the crust is schematically illustrated in Figure 
9. On the other hand, the numerical value of the 
ratio of the two stresses A/B is not safely 
determinable from qualitative considerations. 
Whatever it may be, equations (20) to (23) 
provide the basic solution for the determination 
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CONVECTION CURRENTS 


of the stress distribution in the crust due to the 
action of convection currents in the substratum. 
The selection of the critical constants is as 
yet highly speculative. In a recent paper on 
convection currents, Vening Meinesz (1948, 
p. 200) discusses an apparent example from the 
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pressure is small, say of the order of half, or 
less, of the vertical pressure gradient due to 
weight, then the potential zone of faulting is 
limited to a thin wedge below the surface in 
which numerous shallow slice thrusts are ex- 
pected to occur. If the gradient is sufficiently 
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Indonesian Archipelago and assumes the follow- 
ing constants: depth of crust c = 30 kilometers, 
distance between rising and sinking currents 
1/2 = 340 kilometers, and maximum shear 
stress at the bottom of the crust B = 3,320 
ke/cm®. This case, making the further assump- 
tion A = 0, has been calculated and was found 
to yield a stress system almost identical to that 
shown in Plate 1, B. It is not reproduced here 
since the only essential difference between the 
two pictures consists of a reduction in the 
ratio of thickness over wave length by the 
order of about three to one. A further reason 
for this omission lies in the uncertainty regard- 
ing the proper ratio A/B. The presence of an 
associated vertical stress component being sur- 
mised, the correct solution may correspond 
more probably to a flattened-out version of 
Plate 1, C or to another closely related system. 


CONCLUSIONS 


A superposed horizontal compression with 
constant lateral and vertical gradients is always 
accompanied by a shear-stress component being 
constant in all horizontal planes. The asso- 
ciated thrust-fault surfaces are slightly curved, 
concave upwards for the set dipping towards 
the region of increased compression. If the 
lateral gradient of the superposed horizontal 


large (approximating the vertical gradient of 
the hydrostatic pressure), thrust faults may 
penetrate to great depths, but may then be 
confined to a comparatively narrow zone. The 
same is true if the superposed horizontal com- 
pression has an exponential—instead of a linear 
—trate of decrease. It is then associated with 
vertical and shear stress components likewise 
decreasing exponentially in the horizontal direc- 
tion. The inclination of the resulting thrust- 
fault surfaces diminishes rapidly with depth. 
Hence, a gently dipping, or even horizontal, 
attitude may not necessarily imply a subse- 
quent rotation but can also result directly from 
a local stress system of this type. 

A laterally variable vertical boundary stress 
acting on the bottom of a block is most likely 
associated with a corresponding shear-stress 
component, the latter being 90° out of phase 
if the vertical component is assumed in the 
form of a sinusoidal curve. For certain “bal- 
anced” ratios of the two components, this com- 
bination causes a supplementary horizontal 
tensile stress in the belt subjected to the in- 
creased upward pressure and a corresponding 
compressive stress in the adjacent depressed 
belts. The boundary zones between these belts 
are characterized by absence of a superposed 
horizontal pressure and the occurrence of maxi- 
mum shear-stress components. This causes a 
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45-degree rotation of the stress trajectories, 
which, in turn, leads to a pronounced change 
in the attitude of the fault surfaces. The 
standard normal faulting occurs over the cen- 
tral uplifted portion, but towards both flanks 
one of the two sets gradually changes to very 
steep, vertical, and finally overturned positions, 
whereas the complementary, but rarely realized, 
set of normal faults assumes unusually low 
inclinations. Thus the nearly vertical normal 
faults, so frequently encountered in tectonic 
provinces characterized by differential vertical 
uplift, are explained as a primary feature of such 
stress systems. The anomalous attitudes are 
caused by the described rotation of the stress 
trajectories. They cannot be ascribed to vari- 
ations in the ratio of supplementary horizontal 
and vertical pressures, as is occasionally ex- 
pressed in the literature, based on a false appli- 
cation of the methods of vector analysis. On 
the other hand, as Anderson has pointed out, 
vertical and reverse faults are in numerous 
instances not variations of normal faults but 
belong to the dynamically entirely different 
class of transcurrent faults. 

A pronounced departure of the ratio of the 
two primary stress components in either direc- 
tion from the discussed balanced condition pro- 
duces theoretically interesting cases which may, 
or may not, be realized in nature. One of the 


principal features of these systems is the co- 
existence of two diverse zones of faulting at 
different levels in the crust, one situated above 
the other but separated from it by an inter- 
vening zone of stability. 

Finally, the equations derived provide the 
basic solution for the stress distribution in the 
crust caused by convection currents in the 
substratum. A detailed application to this phe- 
nomenon must await a more secure knowledge 
of the pertinent constants. 
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MEASUREMENT OF PALEOTEMPERATURES AND TEMPERATURES 
OF THE UPPER CRETACEOUS OF ENGLAND, DENMARK, 
AND THE SOUTHEASTERN UNITED STATES 


By H. C. Urey, H. A. Lowenstam, S. Epstern, AnD C. R. McKinney 
ABSTRACT 


Since the abundance of the O" isotope in calcium carbonate varies with the temperature at which it is 
deposited from water, the variation in abundance can be used as a thermometer. This paper discusses the 
following problems: (1) the magnitude of the effect expected, (2) the mass spectrometer of high sensitivity, 
(3) the preservation of the record during geological time, (4) the constancy of the isotopic composition of 
the ocean, (5) the impossibility of using skeletons of air breathing animals, and (6) the temperature of marine 
animals relative to their surroundings. 

A Jurassic belemnite is used to show that the record has been retained since Jurassic times, and belemnites 
of the Upper Cretaceous of the United States, England, and Denmark are used to determine the tempera- 
ture of this time at these localities. The temperatures are about 15°-16°C and indicate nearly uniform tem- 
perature over this latitudinal belt. Because of the possible variation in O¥8 content of the oceans and the 





limited number of samples, these temperatures are regarded as preliminary. 
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INTRODUCTION 


Geologists have drawn many conclusions 
from the purely qualitative evidence of geo- 
logical studies in regard to the past climatic 
conditions on the earth. These deductions are 
based upon a great variety of evidence, and the 
ability of the geologists to deduce as much as 
they have in regard to these conditions excites 
the wonder and admiration of all the uninitiated 
who examine their work even casually. Palm 
trees do not grown in Spitzbergen, according 
to all our evidence, unless the weather is mild— 
about the same as that of Florida. Heavy 
shells are deposited in warm water, and thin 
shells in cold water. Tree rings are definite 
evidence of seasons. 

During the last 15 years a number of people 
have studied chemical properties of isotopic 
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substances. These studies include the demon- 
stration of differences in the chemical proper- 
ties of the hydrogen isotopes, the first calcula- 
tions of such differences in the case of other 
of the lighter isotopes, extensive applications 
of such differences to the separation of isotopes 
and considerable evidence of such differences 
from surveys of the abundances of the lighter 
isotopes in nature. From the standpoint of 
the natural abundances the early demonstra- 
tion that deuterium was not constant over the 
surface of the earth, the demonstration by 
Dole (1935, p. 2731; 1936, p. 268) that oxygen 
also varies, particularly in the air, and the 
variation of carbon from inorganic and organic 
sources by Nier (Nier and Gulbranson, 1939, 
p. 697; Murphy and Nier, 1941, p. 771) con- 
tributed to geology. Recently Thode, Mac- 
Namara, and Collins (1949, p. 361) have dem- 
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onstrated substantial variations of the sulfur 
isotopic abundances, and Thode, MacNamara, 
Lossing, and Collins (1948, p. 3008) have 
demonstrated such variations in boron. The 
present paper presents a way in which the 
natural variations in the abundances of the 
lighter isotopes makes possible, in principle 
at least, quantitative measurement of paleo- 
temperatures. 
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FRACTIONATION OF ISOTOPES IN CHEMICAL 
PROCESSES 


Differences in the chemical properties of the 
elements, such as differences in vapor pressures 
and in equilibrium constants of exchange reac- 
tions, are at least partially responsible for these 
variations in the lighter isotopes in nature. 
Quite certainly the variation of the hydrogen 
and oxygen isotopes in the fresh waters of the 
earth as compared with the sea is due to the 
difference in vapor pressures of the varieties of 
water containing hydrogen and deuterium or 
the three oxygen isotopes, oxygen 16, 17, and 
18. The variation in the abundance of the 
carbon isotopes discovered by Nier cannot 
certainly be assigned to this cause. He found 
that carbon 13 exists in organic matter in a 
decreased abundance as compared with the 
carbon in carbonates. Whether the photosyn- 
thetic process concentrates the carbon isotopes 
by some equilibrium or near-equilibrium proc- 
ess is not certain. Since the oxygen isotopes 
of the air are not in equilibrium with the 
oxygen isotopes of the sea or fresh water, this 
effect is probably also produced by the photo- 
synthetic process or by the oxidation of organic 
matter by atmospheric oxygen,! and if a vital 


1 However, Dole has recently proposed other ex- 
planations for this effect (1949, p. 77). 


selection of the isotopes occurs in the case of 
one of these isotopes it may well do so in 
others. 

Variation in isotopic abundances which are 
due to differences in thermodynamic properties 
—ie., processes which are in equilibrium or 


TABLE 1.—CALCULATED FRACTIONATION FAcrtors 
FOR O'§ anp O!8* 


C COs —H:0() SO. —H20(1) PO. = —H,00) 


0 1.0220 1.0204 


1.0104 


25 1.0176 1.0157 1.0087 
Ratio 1.0044 1.0047 1.0017 
Ko/Kes 





* McCrea (Doctoral Dissertation, Chicago, 1949) 
by attempting to include the effects of the potential 
fields in crystals instead of assuming no such fields 
as was done by Urey (1947) secures slightly differ- 
ent values for the CO;~ — HO constants. 


very close to equilibrium—are understood in 
principle and in many cases in detail. On the 
one hand, it is necessary to know the quantum 
states of these molecules, and on the other 
it is necessary to have the exact theories of 
statistical mechanics in order to relate these 
quantum states to the physical chemical prop- 
erties. For Molecules in the neighborhood of 
ordinary temperatures and above, it is un- 
necessary to consider the rotational quantum 
states, since these are very close to classical 
at these temperatures in the case of all mole- 
cules except those involving the hydrogen iso- 
topes. In this case the entire difference in the 
chemical properties resides in the differences 
in the vibrational energy levels. A great many 
calculations covering exchange reactions of the 
lighter elements were made by Urey and Greffi 
(1935, p. 321) and have been recalculated by 
Urey (1947, p. 562), together with the mathe- 
matical formulae and spectroscopic data by 
which the calculations are made. In general 
differences of a few per cent in the ratio of ise- 
topes in two chemical compounds which are 
brought to equilibrium with each other with 
respect to the exchange of isotopes are possible 
for elements even as high in the periodic system 
as chlorine. The heavier elements can be ex- 
pected to show but very slight effects of this 
kind. The calculations indicate that there 


























should be differences in the ratio of the oxy- 
gen isotopes in water and carbonate ion, water 
and sulfate ion, and water and phosphate 
ion, (Table 1). 

In these calculations the results are reported 
for equilibrium between liquid water and the 
carbonate ion rather than between the gaseous 
water and the carbonate ion. A slight difference 
of a few per cent is expected in the isotopic 
composition, and slight temperature coeffi- 
cients are predicted. The calculations show that 
this difference for the carbonate-water equilib- 
rium shoud amount to about 0.44 per cent 
in 25°C, or 0.176%» per degree C. The calcula- 
tions are not exact because experimental data 
do not supply us with the vibrational energy 
levels for those molecules containing oxygen 18, 
and hence the vibrational frequencies of these 
molecules must be calculated from the known 
vibrational frequencies of the more abundant 
oxygen 16 varieties. The calculation depends 
upon a choice of constants. In fact, in the case 
of the carbonate ion five constants are needed 
to describe the potential energy function to the 
precision of quadratic terms in the co-ordinates, 
whereas only four independent frequencies are 
available for the calculation of these constants 
from the observed vibrational frequencies of 
the CO3"* ion. A similar situation exists in the 
phosphate and sulfate ions. 

If a relation exists between the oxygen 
abundances in these compounds, and tempera- 
ture, and the water in which the compounds 
were deposited, it should be possible to set 
up a thermometer for present-day organisms 
living in the sea, if the isotopic composition of 
the water is known not to differ from the mean 
of the present seas, or, in case it does, if both 
the isotopic composition of the carbonate and 
water are determined. Thus it would only be 
necessary to establish an empirical temperature 
scale and to measure the isotopic abundances 
of oxygen in a shell deposited by an animal in 
the sea in order to determine the temperature 
at which deposition took place. This applica- 
tion may be of interest in some cases, as for 
example in determining the temperature history 
of an animal, but the process which we de- 
scribe is difficult, and wherever possible ordinary 
thermometers are likely to be considerably 
more convenient. On the other hand, these 
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effects may be preserved during geological 
time, and hence this method affords a possible 
method for measuring paleotemperatures. 


PROBLEM OF MEASURING PALEOTEMPERATURES 


Since the predicted temperature coefficient 
for the ratio of oxygen isotopes in calcium 
carbonate relative to water is small, mass 
spectrometers of very great precision and sen- 
sitivity are needed for this problem. It would 
be desirable to measure paleotemperatures 
within less than a degree Centigrade, which 
corresponds to a change in the atomic weight 
of oxygen in the calcium carbonate of 7 X 1077 
atomic weight units, or a change in the abun- 
dance of oxygen 18 by 1 part in 2.5 million of 
oxygen. We have set 0.01 per cent or better, 
corresponding to about 0.5°C or 1°F, as our 
objective in these experiments. As measure- 
ments of the ratio of isotopes have been con- 
ducted in the past, reports in the literature of 
0.1 per cent or 1%p are about as high a pre- 
cision as has been attempted. Hence the first 
problem in the application of this method to 
paleotemperatures is the construction and oper- 
ation of very sensitive mass spectrometers. 

Since we have examples of animals and 
plants laying down isotopes out of equilibrium 
with their surroundings, as for example the 
carbon of organic material, we may well ask 
whether the deposition of calcium carbonate 
by an animal or a plant occurs so as to leave 
the oxygen isotopes in equilibrium with the 
water from which deposition takes place; that 
is, we may ask whether there is a vital effect. 
Probably the equilibrium deposition should 
be closely approximated in many cases, for 
the calcium carbonate deposited is constantly 
bathed with water from the surroundings and 
hence may exchange its oxygen with the sur- 
rounding water during the deposition process. 
However, one could quite readily imagine mech- 
anisms which would involve the deposition of 
organic oxygen and thus invalidate these as- 
sumptions. The whole question could be an- 
swered only by experiment, so far as we can 
judge from the present state of our knowledge 
of the physiology of deposition of such sub- 
stances. The question of deposition of phos- 
phate on the other hand would appear to be 
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even more doubtful than the deposition of 
carbonate, since phosphate ion does not ex- 
change its oxygen with water in which it is 
dissolved, and hence perhaps the phosphate 
deposited by animals generally does not have 
its oxygen in equilibrium with its surroundings. 
Since phosphate is so intimately involved in 
biochemical reactions, its oxygen may be 
brought to equilibrium with water by them. 
Only careful experimentation can determine. 
the feasibility of the phosphate thermometer? 

Animals rarely deposit sulfate, and it is not 
an important substance for investigations of 
the type we envisage immediately. Silica is 
often deposited by plants and animals as a 
hydrous silica. Since the relative proportions 
of water and silica might change the relative 
abundance of the oxygen isotopes, this does 
not appear to be a very promising line of 
attack, though again experiments should be 
made to determine whether there is any possi- 
bility of using a silicate thermometer. Of the 
possible thermometers—that is, the carbonate, 
phosphate, sulfate, and silica the carbonate 
thermometer seems the most likely to give 
satisfactory results. A combination of the car- 
bonate and phosphate would eliminate the 
water, and hence the isotopic variation of water 
always present in brackish water and fresh 
water would be eliminated. Dr. Martin Stein- 
berg is investigating the possibilities of this 
thermometer. 

Assuming that the shell is laid down in 
equilibrium with the water, the next important 
question is whether this composition is pre- 
served during geological time. In the first 
place, if radical metamorphosis takes place, 
as for example the conversion of limestone to 
dolomite, no satisfactory record of the isotopic 
abundance can be expected, for the disolved 
carbonate should exchange with the water in 
which it is dissolved and then be redeposited 
with the composition partially dependent on 
its original composition and partially on the 
water which was the means of redeposition. 


2Dr. Mildred Cohen (private communication) 
informs us that she has found that some enzymatic 
hydrolyses of phosphoric esters split the oxygen 
phosphorus bond though others do not. Hence a 
mechanism is available which may establish equi- 
librium between water and phosphate in animal 
bodies. 


Also, if aragonite is converted to calcite through 
solution in water, the record would be de- 
stroyed, but, if recrystallization without the 
possibility of exchange occurs, the record might 
be preserved. 

The second important method for destruc- 
tion of the temperature record is diffusion 
through the solid state. Shells as usually de- 
posited appear to be fairly compact, but micro- 
scopic examination shows that they consist 
mostly of very small crystals, often with definite 
canals through the shell, as for example in the 
case of the brachiopods or the echinoderms, 
Thus in considering the problem of diffusion 
it is not the gross size of the shell that is im- 
portant, but rather the small microscopic par- 
ticles, since in the course of geological time 
water would probably pass through the shell 
by hydrodynamic flow, or diffusion in the liquid 
would bring the isotopic composition of water 
within the shell to that of water in which it is 
immersed. The diffusion problem for consider- 
ing this situation exactly would be very com- 
plicated. Assuming that the initial conditions 
are that (1) at the surface of the crystals the 
composition is that required for equilibrium 
with the surrounding water, and that it remains 
so throughout the diffusion process, and (2) 
that initially the composition of the crystal 
differed from that which would be in equilib- 
rium with the water in which it has been 
immersed during the time since the shell was 
formed, and that it changes from this initial 
condition in accordance with the diffusion laws 
at some constant temperature, a calculation 
can be made which relates the size of the crys- 
tal, the time, and the coefficient of diffusion 
to the composition of the shell remaining. The 
fraction of the excess mean concentration over 
that in equilibrium with the surrounding water, 
assuming that the coefficient of diffusion is the 
same in all directions, is 


—_= = > —____— e7a(tmt+)) .3 
(2m - 1)? 


} 


where a = D is the diffusion coefficient, 


Dx'*t 
ie 
t the time, and / the length of the edge of the 
cube. 
If we require that the crystal shall retain its 
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original composition within some arbitrary per- 
Drt 





centage, the value of can be estimated 


from this equation. (If one uses different values 
of D for different axes, the result is the same 
as though different lengths were used for the 


TABLE 2.—VALUES OF Q/Qo FOR VALUES OF a* 














Dri Q 
ace’. Qo 
0.1 0.461 
0.05 0.590 
0.01 0.795 
0.005 0.854 
0.001 0.915 
0.0001 0.964 
0.00001 0.972 





* We are indebted to Mr. L. G. Montet for the 
numerical results given in this table. 


direction parallel to the axis and for the two 
directions perpendicular to this.) 

Table 2 gives a few values of a for a few 
values of Q/Qo. 

The coefficient of diffusion of one ion can be 
calculated from the electrical conductivity pro- 
viding the current is carried by that ion and 
not by the other. Joffe (1923, p. 461) concluded 
that only the calcium ion carries the current, 
but to what extent the carbonate ion carries 
current was not determined except that its 
contribution was not detected. If we assume 
that the carbonate ion carries some fraction, 3, 
of the current, the coefficient of diffusion of 
the carbonate ion can be calculated from the 
resistivity of calcite from the equation, 


os RTé 
~ mtF*Cp’ 


where m is the number of charges on the ion— 
i. 2—, F is the faraday, 9650 EMU, C is the 
number of gram molecules of carbonate ions 
per cm, p is the specific resistance in absolute 
units, and 6 is the unknown fraction of the 
conductivity due to carbonate and assumed 
to be 0.01 here. If p is in ohms, this formula 
must be multiplied by 10-*. Curie (1889, p. 
385)? measured the specific resistance of calcite, 





3 See Landolt-Bornstein Tabellen II, 1062 (1923). 


and his values and the values of D calculated 
from them are given in Table 3. 

Using the first two values of Doo3 in the 
table, we can estimate the time required for 
diffusion to destroy the record. Thus at 20°C 
a crystal of 1 mm. dimensions will retain 96.4 


TABLE 3.—MEASURED SPECIFIC RESISTANCE OF 
CALCITE AND CALCULATED COEFFICIENTS OF 











DIFFUSION 
eC p Doo; 
|| axis 20 | 5.5 X 10 |4.4 x 107% 
100 | 5 X10" |4.8 x 10-™ 
160 |3 X10" |8.1 x 10- 
Laxis 15 | 9.5 X 10 |2.5 x 107% 
100 | 2.4 X 10" |1.0 x 10-” 
150 | 1.3 X 10" |1.9 x 107» 














per cent of the original concentration relative 
to the equilibrium concentration of a changed 
environment for 2.3 xX 10% seconds (7 « 108 
years), while at 100°C the same crystal will 
retain this record only 2.1 xX 10" seconds 
(6.4 X 104 years). Crystals of this size will not 
lose the temperature record due to diffusion 
in times sufficiently long to be of considerable 
interest, and the coefficients of diffusion may 
be even smaller than those given in Table 3. 
A shell may grow a centimeter per year or 
3 X 10-* cm. per day within an order of magni- 
tude. For simplicity we assume that the con- 
centration of O% with distance due to daily 
variation of temperature varies according to a 
cosine law. Then the diffusion equation gives 


x 
C = e-D#'tll? cos x ; ‘ 


where we take the mean concentration as zero. 
If we can detect a concentration difference 
equal to 0.5 of the original difference between 
maximum and minimum, 7.e., at x = 0 and 
x = 1, we have 


Dxr't 
2= —. 
In p 


Substituting D = 4.4 X 10% and/ = 3 x 
10-*, the time becomes 1.4 X 106 or half a 
billion years. If we ask for 0.75 of the original 
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variation the time becomes fifty milion years. 
Thus even a daily variation of temperature 
might be partially preserved. 

There is the question of whether the isotopic 
composition of the ocean has remained con- 
stant during geological time. Such a change 
would be produced, if igneous rocks of one iso- 
topic composition are weathered and sediment- 
ary rocks of another composition are produced. 
At high temperatures the isotopic compositions 
of all chemical substances should be the same, 
and, if this temperature was that of molten 
lava or higher, differences in isotopic composi- 
tion of water and silicates should not be large 
and at least they should be less than at lower 
temperatures. The calculations referred to 
above show that O" concentrates preferentially 
in the solid phases and not in the water, and 
it can be expected that this will be true with 
respect to the silicate rocks. Thus 0" should 
enter the sedimentary rocks and its concentra- 
tion in the oceans should decrease with time. 
Assuming that the ratio of concentration in 
water and silicates is the same as for sulfate 
and water at 25°C, namely, 1.016, and that 
600 g. of rock have been eroded for each 1000 
g. of water (Goldschmidt, 1933, p. 112), it is 
easily shown that the water should have 
changed its ratio of oxygen isotopes by 3.8%p» 
due to this cause. Dr. P. Baertschi has shown 
that there are differences between ocean water 
on the one hand and the Disco Island Green- 
land basalt containing native nickel-iron and 
other igneous rocks, and, on the other hand, 
sedimentary rocks of just the kind and magni- 
tude predicted, though the total difference 
between ocean water and sedimentary rocks 
is not 16%» as predicted. These results will 
be published elsewhere. 

It is of interest to inquire as to the time at 
which these changes took place. The carbonate 
equilibrium should have taken place early in 
the earth’s history, and largely since then no 
important change in the O* content of carbon- 
ates occurred. The shales which make up much 
of the sedimentary rocks should exchange their 
oxygen with water less readily than the car- 
bonates, but no information is available. Baer- 
tschi finds that sandstones have an increased 
O08 content relative to the Greenland basalt 
and thus may have exchanged their oxygen 
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to some extent, though the difference is small 
enough to be explained by equilibrium differ. 
ences in concentrations between the silica and 
other phases of igneous rocks. Perhaps it is a 
conclusion based on little more than prejudice, 
but we believe that it is most likely that the 
major part of this change occurred before the 
Cambrian. Until more data are available we will 
make no correction for this effect. 

If the water of the oceans has been steadily 
produced during the earth’s history from plu- 
tonic activity, as maintained by some (W. W. 
Rubey, private communication), the entire 
problem of the isotopic composition of the 
oceans must be considered from this point of 
view and must involve measurements on the 
O* and D concentrations in waters which might 
be expected to be juvenile, as from geysers, 
hot springs, and volcanoes. Such measurements 
have not been made, but are desirable from the 
standpoint of the work presented here as well 
as for the evidence they might give relative 
to this suggestion in regard to the origin of the 
earth’s water. 

The oxygen of the air has a higher isotopic 
composition than that of the sea and fresh 
water. The difference in isotopic composition 
amounts to about 2 per cent for air oxygen as 
compared with the sea. All animals use atmos- 
pheric oxygen, and hence produce water that 
contains an increased concentration of oxygen 
18. In the case of marine animals having gills 
which exchange water between the sea and their 
blood stream, it can be expected that in spite 
of the oxidation processes taking place within 
their bodies the water content of their bodies 
will be very close to that of the sea. One can 
reason that this must be true, for the water 
oxygen must pass through the gills approxi- 
mately as readily as does the elementary oxygen, 
and since the amount of water present at the 
surface of the gills is larger than the amount of 
oxygen by a factor of about 105, the rate of 
transfer of the water oxygen must be much 
more rapid than the transfer of elementary 
oxygen, and hence the isotopic composition of 
the blood stream must remain close to that of 
the sea. On the other hand, either marine or 
land animals that breathe air are suspect for 
our purposes, since atmospheric oxygen is high 
in oxygen 18 (i.e., plus ~ 20% ). Land animals 
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drink fresh water with a decreased concentra- 
tion of oxygen 18, and animals of the sea eat 
food containing a very considerable amount of 
water of the composition of water of the sea. 
Without knowing the relative importance of 
these factors it is difficult to have any confi- 
dence in any measurements made on air-breath- 
ing marine animals, and even less on land 
animals. It would be possible to investigate 
the blood of marine reptiles or mammals to 
determine how nearly the oxygen 18 content 
of their blood agrees with that of the sea. 
This has not yet been done. 

The question arises as to whether marine 
animals maintain their body temperature above 
that of the surrounding sea. They do not. 
The source of heat is the combustion of food 
by oxygen received through the membranes of 
gills, usually, or of the general body surface. 
The membrane through which the oxygen dif- 
fuses is the same as that through which heat 
is lost to the surroundings. The transport of 
oxygen determines the rate of heat production, 
and this must equal the rate of loss by thermal 
conduction. Assuming the maximum concen- 
tration of oxygen in sea water—.e., 6 ml. per 
liter—and that the oxygen combines with sugar 
to produce 56,000 calories per gram atom of 
oxygen, and knowing the coefficient of diffusion 
of oxygen in water and the coefficient of con- 
duction of heat in water, we can set up the fol- 
lowing equation: 


AC At 
56,000D — = K — 
,000. i w’ 


where D is the coefficient of diffusion, K the 
coefficient of thermal conduction, AC the con- 
centration difference of oxygen between the 
blood stream and the water and not greater 
than that corresponding to 6ml. per liter, or 
0005 gram atoms per cm’, Af the difference 
in temperature between the water and blood 
stream, and Al the membrane thickness through 
which both diffusion and heat conduction take 
place. The quantity Al is unknown, but is the 
same for both processes and can be cancelled 
from the equation and the equation solved for 
At. Substituting D = 2 X 10-5 cm*sec~ and 
K = 1.42 x 10-* cal. cm.sec~!°°, we find that 
4t = 0.38°C. This is the maximum difference 
in temperature possible under normal aerobic 


PROBLEM OF MEASURING PALEOTEMPERATURES 


405 





conditions. Animals can expend energy for 
limited periods of time under anaerobic con- 
ditions and thus raise their temperature above 
that calculated here, but such periods would 
seem to be negligible. 

The preparation of carbon dioxide gas for 
use in the mass spectrometer has constituted 
one of the most difficult problems in this re- 
search. It is so intimately related to the prob- 
lem of the empirical temperature scale that 
both these problems and related ones will be 
discussed in a separate paper. Impurities in the 
carbon dioxide must not produce peaks at 
masses 44 and 46 such that the ratio of inten- 
sities of these masses is changed by more than 
one part in 10‘, which means an impurity in 
mass 44 of less than this amount and in mass 
46 of less than one in 2.5 X 10°. We are for- 
tunate in using a relatively clean part of the 
mass spectrum. The carbon dioxide is prepared 
from calcium carbonate by the action of ortho- 
phosphoric acid to which phosphorus pentoxide 
is added to bring the concentration to approxi- 
mately 100 per cent. To this enough chromic 
oxide is added to give a slight permanent red 
tint, and the excess oxide is reduced with 
hydrogen peroxide and the excess decomposed 
by heating. The calcium carbonate sample is 
first heated in a stream of carbon dioxide-free 
helium for 15 minutes at approximately 475°C 
in order to destroy organic matter. Approxi- 
mately 9 cc. of the phosphoric acid and 50 mg. 
of the calcium carbonate are put in separate 
arms of a reaction vessel, evacuated, and 
brought to 25°C in a thermostat. The acid is 
then poured onto the carbonate. A vigourous 
evolution of carbon dioxide results, sometimes 
with a little froth which breaks during half an 
hour in the thermostat. The gas is then trans- 
ferred to sample tubes and analyzed on the 
mass spectrometer. All stop cocks and tapered 
joints are greased with Apiezon N cock grease. 
We have never detected any effect of impurities 
due to this material. 

With this procedure reproducible analyses 
are secured for both recent specimens and fos- 
sils. Using the shells of present-day marine 
animals grown at known temperature and cor- 
recting for the variation of isotopic composition 
of the various sea waters, the following equa- 
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tion for temperature is secured (Epstein, 
Buchsbaum, Lowenstam, and Urey, 1950): 


# = 11.88 — 5.916, 


where ¢ is measured in °C, and 4 is the differ- 
ence in % > between the ratio of O'* to O'* in 
a standard working gas, and this ratio for the 


sample, i.¢., 
-1) 10, 


pare ( Reample 
Retandard 

where Rompe and Ryrandara ate the ratios of 
COO" to CO’ for the sample and reference 
gas respectively. The reference gas is carbon 
dioxide prepared from Belemnitella americana 
from the Peedee formation of South Carolina. 
The slope of this curve agrees within the limits 
of error with that secured by McCrea (1949; 
1950, p. 849) from inorganic crystals of calcium 
carbonate grown at constant temperature. This 
equation is believed to be accurate to +1°C 
if the calcium carbonate were deposited in sea 
water of the mean isotopic composition of sea 
water of salinity 34.8% . If the calcium car- 
bonate was deposited in water of other isotopic 
composition, appropriate corrections must be 
made. The maximum variation in temperature 
for fully marine conditions due to the variation 
of isotopic composition of the oceans is +4.5°C. 

The variation of isotopic composition of the 
present oceans, as discussed by Epstein, Buchs- 
baum, Lowenstam, and Urey (1951), presents 
the greatest difficulty so far encountered in the 
measurements of paleotemperatures. These au- 
thors point out that isotopic composition varies 
approximately linearly with salinity. Their data 
on surface waters are not very extensive, but 
this rough linear relation applies also to surface 
waters. The isotopic composition varies in a 
rough linear way with the temperature of sur- 
face waters. This indicates that, if no correction 
for the isotopic composition were made, a tem- 
perature of 30°C would be read as about 25°C, 
while a temperature of 10°C would be read as 
15°C. High temperatures would be read too 
low, and low temperatures too high. However, 
no correction to past temperature determina- 
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tions can be made by use of their curve because 
(1) A general average change in the earth’s 
temperature would not change the average 
isotopic composition of the oceans, and hence 
a measured isotopic composition could not be 
corrected with their curve. (2) The slope of the 
curve of 08 composition vs. temperature is due 
to evaporation and condensation of water, and 
hence can be expected to be greater in times 
of high precipitation and lower in those of 
low precipitation. In general one would expect 
that high rates of evaporation and precipitation 
would accompany higher average temperatures 
of the earth’s surface. However, we have the 
phenomena of glacial periods which somehow 
appear between limits of extremely low and 
extremely high temperatures, and these make 
estimates in regard to the relation of tempera- 
ture to precipitation very difficult. For the 
present we do not attempt to estimate the 
correction for isotopic variations in the sea, 
but this problem is being studied from several 
angles and it is hoped that corrections can be 
applied in the future. 


MEASUREMENTS OF IsoToPIc COMPOSITION, 
AND TEMPERATURES OF THE 
Upper CRETACEOUS 


In considering geological specimens for our 
first attempts at measuring these temperatures, 
we have been guided by our considerations of 
probable diffusion of ‘the carbonate in solid 
calcium carbonate and the improbability that 
recrystallized carbonates would retain the 
record. Since aragonite is unstable and com- 
monly recrystallized, all recrystallized arago- 
nitic shells will probably lose the record. 

A study of fossil carbonate shells indicated 
that the belemnites would be the most likely 
to retain the record. These have a very compact 
structure with calcite crystals megascopically 
visible and so closely fitted together that in- 
ternal reflections do not occur since the speci- 
mens are characteristically translucent. The 
polarizing microscope shows that the calcite 
crystals are arranged with their axes perpendic- 
ular to the axis of the guard. The orientation 





Pirate 1.—SEASONAL GROWTH OF JURASSIC BELEMNITE 
, Figure 1.—Cross SECTION OF A JURASSIC BELEMNITE . 
Showing relation of samples to growth rings. W and S refer to winter and summer regions respectively. 
FIGURE 2,—VARIATION OF TEMPERATURE AND C! CoNTENT WITH RADIUS OF JURASSIC BELEMNITE 
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respect to carbon and oxygen. A flat disk 
about 3 mm. thick was cut from the specimen, 
polished, and photographed in transmitted light 
(Pl. 1, fig. 1). Samples of carbonate were 
ground from the specimen following the growth 


of the crystal axes and the shape and size of the 
crystals are evident. 

Several paleontologists advised investigation 
of the possibility that various representative 
species of the genus Imoceramus might retain 


TABLE 4.—ANALYSES OF A JURASSIC BELEMNITE 























6(C18) 5(O'8) 
Ring No. eC 
R L Av. R L Av. 

1 3.34 2.50 —0. 65 —1.8& —0.65 (?) 15.7 
E 3.57 3.35 3.4; —0.4, —0.9; —0.7; 16.; 
3 3.8 3.% 3.8 —0.4, —0.3o —0.3s 14.; 
4 2.65 2.88 2.76 —1.0, —0.8 —0.9, ° 
5 2.41 2.37 2.39 —1.3, —1.3, —1.3; 19., 
6 3.27 3.6 3.45 —1.3, —1.1; —1.2. 19.; 
7 4.2. 4.3, 4.3 —1.2; —1.1, —1.1, 18.5 
8 3.52 _ 3.5 —1.1, _ —1.1, 18.¢ 
9 2.65 _ 2.65 —1.04 _ —1.0, 18.¢ 
10 —_ 2.75 i _— —1.0; —1.0; 18., 
11 3.46 3.62 3.54 —0. 62 —0.4, —0.5; 15., 
12 3.25 3.49 3.34 —1.0; —0.8; —0.% 17. 
13 3.37 3.37 —0.8 —0.8 17.; 
14 _ _ — — _— 
15 3.06 | 3.08 3.07 —1.5, | —1.3; —1.4, 20.3 
16 3.15 3.15 —0.7; —0.7; 16.; 
17 3.0, | Ske 3.12 —0.8 | —0.64 —0.7; 16.2 
18 2.75 yy A —0.75 —0.7; 16.3 
19, 20 2.7% ay —1.2; —1.2; 19.9 
21 2.59 2.59 —1.2, —1.2, 19.2 
22, 23 2.66 2.66 —0.8 —0.8, 17., 
24 2.65 2.65 —1.3o —1.3o 19.4 
Av. 17.6 




















the temperature record. These fossils show 
evidence of massive calcite crystals of substan- 
tial size, though they are opaque. The calcitic 
layer of theshell is also present in many present- 
day mollusks, though usually not developed to 
the remarkable extent that it is in this extinct 
group. Experiments in progress confirm our 
expectation that this group retains the temper- 
ature record. 

In the belemnites, a particularly convincing 
test for the durability of the record has been 
made. A Jurassic belemnite* approximately 2.5 
cm in diameter, brownish and translucent, 
with very well-defined growth rings, has been 
investigated for a seasonal effect and bilateral 
symmetry of isotopic constitution both with 





‘Unfortunately, we do not know the species of 
this specimen. It was collected along with others on 
the Island of Skye by Dr. Cyril S. Smith, who kindly 
made them available to us. The specimens are of 
various shapes and sizes and do not appear to be a 
single species. As nearly as we can learn they are 
Jurassic, probably from the Oxford clay. 


lines from each side of the disk, as shown by 
the numbered rings shown in the photograph. 
In some cases right and left samples were not 
taken and analyzed, either because of accident 
in the cutting or preparation or because the 
samples were combined to secure sufficient 
carbonate for our samples. In one case both 
the right and left samples were lost. 

The samples were analyzed for the C* and 
O" isotopes (Table 4). 

The right and left samples of the outermost 
ring do not agree either in C" or 0%, and the 
check for the two samples of the second ring 
is rather bad for the O*. Otherwise the agree- 
ment on the two sides lies within 0.2%. 
Considering probable errors of measurement, 
and the difficulty of cutting samples of the 
same thickness from the two sides, the agree- 
ment is very satisfactory. A radial variation of 
isotopic composition occurs both in the carbon 
and oxygen isotopes. The variation of the 
carbon isotopes does not correlate with any 
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other variable found either in fossil or present- 
day specimens, and the variation in this speci- 
men is of the order found in present-day 
specimens. The erratic variation is probably 
due to fractionation of the carbon isotopes by 
plants and the large effect that this has on the 
very small amount of carbon in the sea (ap- 
proximately 28 mg. per liter). 

Figure 2 of Plate 1 plots temperatures deduced 
from our O¥ abundances against a radius meas» 
ured along a line from the center of growth 
outward. The temperatures are plotted as hori- 
zontal lines over the range of the radius for 
which they are the mean temperature. The 
true temperature curve should intercept equal 
areas above and below these lines. Such a 
curve is drawn trying to take account of this 
requirement and the uncertainty of the data. 
Where analyses for both the right and left 
sides are available the line is the average, and 
points above and below represent the two data. 
The plot of & (C¥) is also given and shows no 
correlation with the seasons. The lack of struc- 
ture at the smaller radii is probably due to the 
fact that the samples are too thick to show the 
variations that may exist. 

It seems difficult to believe that any “acci- 
dental” process of diffusion or exchange could 
have produced such a distribution of carbon 
and oxygen isotopes, and we conclude that the 
record has been substantially preserved for over 
a hundred million years by this belemnite. In 
our work on other belemnites, we have never 
had any reason to doubt the isotopic record 
providing the specimen was translucent. 

This Jurassic belemnite records three sum- 
mers and four winters after its youth, which 
was recorded by too small amounts of carbonate 
for investigation by our present methods; 
warmer water in its youth than in its old age, 
death in the spring, and an age of about four 
years. The maximum seasonal variation in 
temperature is about 6°C. The mean tempera- 
ture was 17.6°C. 

Dr. C. J. Stubblefield of the British Museum 
has very kindly co-operated in securing speci- 
mens of belemnites, oysters, and brachiopods 
from the Upper Cretaceous Chalk of England. 
Thus we have well-preserved specimens im- 
bedded in a matrix which has been modified 
but slightly with time. Table 5 lists the speci- 
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mens secured from him and his record of the 
places and zones from which they originate. 

Dr. J. Troelsen of the Mineralogical and 
Geological Museum of Copenhagen also sup- 
plied us with three belemnites from the lower 
and upper Maestrichian of Denmark. The data 
regarding their collection are as follows: 

1. Belemnitella mucronata, Upper Maestri- 
chian, 17m. below the top of the white 
chalk, Stevns, Denmark. 

2. Belemnitella mucronata, Upper Maestri- 
chian, 30-40 mm. below the top of the 
white chalk, Stevns, Denmark. 

3. Belemnitella mucronata, Lower Maestri- 
chian (white chalk), Mgen, Denmark. 

In studying the isotopic composition of these 
specimens we first removed carefully the chalk 
from the specimen. In the case of the belemnites 
the surface layers and any chalky infiltrations 
in cracks of the specimen were carefully re- 
moved. When possible a complete cross section 
of the guard was used to secure a true average 
temperature. In any case a piece of the solid 
clear calcite was gound to a fine powder and 
carbon dioxide prepared from these samples. 
In the brachiopods and oysters the chalk again 
was carefully removed, and parts of the clean 
shells were gound to fine powder and used 
for the preparation of the carbon dioxide sam- 
ples. Whenever the chalk was available samples 
of the chalk in the immediate neighborhood 
of each specimen were analyzed. The results 
of the carbon and oxygen analyses for the 
English samples are recorded in Tables 6 and 
7. In these tables the zones have been numbered 
as indicated in Table 5. 

The oysters, brachipods, and their chalks 
agree in O* composition only very roughly, and 
their C" concentrations do not agree within 
the limits of our measurements. The belemnites 
do not agree in O¥ concentrations with their 
chalks nor with the oysters and brachiopods. 
We conclude that only the belemnites retain 
the temperature record. Probably the crystals 
of the other specimens are so small and the 
structure so open that diffusion has destroyed 
the record. In other studies now in progress we 
have found that a chalky belemnite—i.¢., one 
that has lost its translucent character—agrees 
in isotopic composition with its surrounding 
chalk. 
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TABLE 6.—TuHE Qvantity 4(C) ror UPPER 
CRETACEOUS OF ENGLAND 




















Brachiopod Oyster | Belemnite 

od amy RE Oe RR Ny le 
—_ | Chalk Speci- | Chalk | Spon | Chalk 

1 | 1.3) | 0.9% | 0.7) | 0.7, | 0.45 | 1.6, 
F 2.83 1.8; 2.60 | 1 5s 2.4: 1.3, 
3 2.53 | 1.61 | 1.8 | otek Bade | Bude 
412.9% | — |2.2|1.%/0%| — 
6124! thi thelial — | — 
6 | 2.3: | 1.67 | 1.8 | 1.6 | 1.8 | — 
/ —_ — | — —_ — 
8 | — | 1.5. ] 1.6 | 1.52} — — 
9 1.53 | 0.75 | 2.0, | 1.0, —_ —_ 
10 1.97 | 1.17 | 1.30 0.95 | — _ 
mu | — | — }1.%/)/1.2}); — | — 
2)/—-i;—-|;—-{]—-/;-|]-— 
13 —_— —_— 2.8; 2.8: | 3.20 _ 




















Some specimens of Table 5 were not used because 
of bad preservation. 


constant composition had been reached, all 
the analyses should be identical, and thus we 
conclude that this is not the case. 

The oxygen analyses for the Danish speci- 
mens together with derived temperatures are 
given in Table 8. The specimens are numbered 
as given previously. 

Upper Cretaceous paleotemperature data for 
North America were secured from the Exogyra 
costata zone of a number of widely scattered 
outcrops of the Gulf Coast and Eastern Sea- 
board of the United States between western 
Tennessee and northeastern Mississippi north- 
eastward to New Jersey.* The specimens from 
Tennessee were derived from the sandy marls 
of the basal third of the exposure on Dave 
Weeks’ place, 34 miles south of Enville, Mc- 
Nair County, the type locality of the Coon 
Creek tongue of the Ripley formation (Wade, 





TABLE 7.—THE Quantity 6(0!*) ror Upper CRETACEOUS OF ENGLAND 














. Brachiopod | Oyster | Belemnites | Belemnite 
one 
: 3 
Specimen | Chalk | Specimen | Chalk Specimen | Chalk | ‘emperature 
| | | | i 
1 —1.46 —2.35 —1.85 | —2.2. —0.5;5 —1.8s 15.0 
2 | 1.9 —2.9 | -2.6 | 2.9 | O42 | -2.52 | 14a 
3 —2.40 —2.4, | —2.9, } —2.66 —1.2, —3.8; | 19.2 
4 | —1.95 —_ —1.9 | —2.6¢ | —1.2¢ _ 19.3 
5 —2.4; —2.2 | —2.73 —2.7%4 —_ _— 
: —2.1o —2.4¢ —2.4, —2.56 —2.05 -- 23.8 
si — | =-2@ | -—2m 2% -- oe: ty 
9 —2.83 | —2.9% —2.02 —2.5, — — 
10 | —2.60 —2.4, —3.6s —2.7; | _— _— 
11 | = | ~ —3.62 —3.16 | —_ — 
12 | -” | - = — ate | om | 
3); -—- | =— —1.8 —3.54 —0.40 | — | 14 
14 | mane | sca oa —_ — = 





Present-day specimens show marked varia- 
tion in C** concentrations even when they grew 
in the same region at the same time. Thus the 
samples considered here should have had iden- 
tical or nearly identical O* concentrations, but 
not identical C' concentrations when first 
buried. If both the oxygen and carbon con- 
centrations changed with time we would expect 
that the oxygen concentrations would tend to 
change by comparable amounts, and hence 
remain somewhat the same as observed. On the 
other hand the carbon concentrations should 
show no such approximate agreement. If com- 
plete equilibrium with water and carbon of 


1926). Though belemnites do not occur in these 
deposits, depriving us of our reliable standard 


5 We have analyzed a few samples of ooze from 
the Pacific and Atlantic Oceans in co-operation with 
W. F. Libby and J. R. Arnold, who are working on 
C* concentrations in such materials. We are unable 
to interpret the results in any consistent way as yet, 
and evidently more data are required. : 

6 The specimens from Western Tennessee, Mis- 
sissippi, North Carolina, and South Carolina were 
collected by H. C. Urey and H. A. Lowenstam. We 
are indebted to Dr. L. W. Stephenson for advising 
us on certain collecting localities in South Carolina, 
to Dr. Horace C. Richards who guided us to the 
New Jersey outcrops, and to Mr. P. A. Bethany of 
Macon, Mississippi, who advised us on and collected 
with us from the outcrops of the Macon area. 








_ 


nn thie ie 








S we 


peci- 
3 are 
ered 


a for 
gyra 
ered 
Sea- 
tern 
rth- 
rom 
aris 
dave 
Mc- 


ade, 








ISOTOPIC COMPOSITION AND TEMPERATURES OF UPPER CRETACEOUS 


for comparison, the excellent shell preservation 
of all the molluscan elements with originally 
aragonitic shells, elsewhere in the outcrop areas 
preserved only as casts, implied exceptionally 
favorable burial conditions suggestive of pres- 
ervation of original isotopic abundances. The 
specimens analyzed from northeastern Missi- 
ssippi were secured from the Prairie Bluff 
Chalk, at two localities, the glades on the Vanny 
Jackson farm west of the Rocky Hill Church 
road, 3.6 miles north-northeast of Starkville, 
Octibbeha County, which is the general area 
of fossil locality 6846 of Stephenson and Monroe 
(1940), and from glades on the A. Lindley farm 
4 miles southwest of Macon, Noxubee County 
(Stephenson and Monroe’s 1940 collecting lo- 
cality 17242, 17484). At both localities, the 
specimens were found loose in the residual 
fossil concentrations which littered the glades, 
since Prairie Bluff specimens satisfactorily pre- 
served and located im situ were unavailable. 
The positions marked 10 feet and 13 feet above 
the base of the Prairie Bluff formation in Table 
9 for the specimens from the Lindley locality 
refer to collecting levels on the glade and not 
to stratigraphic positions. Despite the rich 
assemblages in the glade concentrations, only 
originally pure calcitic skeletal forms on the 
calcitic layer of calcite and aragonitic inter- 
layered shells are preserved as such, reducing 
the number of potential elements for paleo- 
temperature determinations to crinoidal re- 
mains, echinoids, Terebratulina floridana, Os- 
treidae, and Pectinidae among the pelecypods 
and guards of Belemnitella americana, the bal- 
ance consisting of molds, prevalently in phos- 
phatic form. Shark teeth occur with the inver- 
tebrate assemblages composed of carbonate 
skeletons and furnish source material for a 
possible future check of the carbonate ther- 
mometer by the phosphate thermometer. 

The South Carolina specimens were secured 
in situ from the basal section of the Peedee 
formation which is well exposed in the west 
bank of the Peedee River at Burches Ferry, 
154 miles from Florence. This is the locality 
from which Stephenson (1923) described and 
figured a number of well-preserved specimens 
of Belemnitella americana and described a new 
species. At the time of our visit, the river was 
low enough to expose a total of 134 feet of Pee- 
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dee strata, consisting of calcareous, argillaceous, 
glauconitic sand above 34 feet of the shales of 
the Black Creek formation. Belemnitella ameri- 
cana was found in small numbers through the 
entire local Peedee section; it was concentrated 
in conspicuous numbers in a narrow 9-inch 
zone, 57 inches above the base of the formation. 


TaBLE 8.—OxyGEN Isotopic ANALYSES AND 
TEMPERATURES OF THE DANISH 














BELEMNITES 
Specimen | 5(018) | a 
1 | —hie | the 
2 | —0.4; 14.3 
3 —0.05 12.; 





The Belemnitella guards are throughout oriented 
more or less horizontally, with occasional de- 
viations ranging up to about 10° from the hori- 
zontal plane. Within the horizontal plane the 
belemnite rostra display random orientation 
as far as could be ascertained in the vertical 
exposure. Exogyra costata occurs sparingly dis- 
tributed through the basal 7} feet of the 
section, the individuals being generally small 
and thin shelled. Larger, but comparatively 
fragile individuals are abundant in the over- 
lying 22 inches of the section which is con- 
spicuously fossiliferous, consisting of a semi- 
coquina of pelecypod shell debris. The speci- 
mens analyzed were taken in situ. 

The North Carolina specimens were collected 
from a 6-inch calcareous sandstone ledge of the 
Peedee formation, exposed above water level, 
and 6 inches of soft marl below water level at 
the Jackson’s Hole Landing, locality of Ste- 
phenson (1923, p. 25), Northeast Cape Fear 
River. A number of large, heavy-shelled Ex- 
ogyra costata individuals, with strongly resorbed 
surfaces, and a single rostrum of Belemnitella 
americana, were secured in situ from the water- 
covered marly layer, which contrasts strongly 
with the rich mollusk assemblages of overlying 
calcareous sandstone ledge. Pectinidae and 
Belemnitzlla americana, originally calcitic skele- 
tal types alone are well preserved; the originally 
aragonitic shells of the associated pelecypods 
is often represented by casts and molds. As the 
Peedee section lacks top and bottom it is not 
possible to determine the stratigraphic position 
of this section within the formational range. 
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TABLE 9.—OxyceEN Isotopic ANALYSES AND TEMPERATURES FROM Belemnitella americana FROM THE 
SOUTHEASTERN UNITED STATES 





—__—_—___.. 





| Av. 























Formation Locality Level Analysis | Temp. | Temp. 
Prairie Bluff | Starkville, Miss. 3.6 miles | Unknown (1) —0.5s | 14., 
N-NE. on Rocky Hill (2) —0.97 | 17.6 
Church Road (3) —1.2; 19.; 
(4) —1.0, 18.2 
(5) —0.1; 12.5 | 16., 
Macon, Miss. A. Lindley | 13 feet above base of Prairie (1) 0.00 Sie.) Bie 
place 4 miles SW. of Bluff 
Macon 
10 feet above base of Prairie (1) —1.0; 17.5 
Bluff (2) —0.75 16.2 | 
(3) —0.5; 14.5 | 
(4) —1.4, 20.6 | 17., 
Peedee Peedee River, S. C. 105 inches above base of (1) —0.5; 14.5 | 14., 
formation 
57 inches to 68 inches above (1) —0.4, 14., 
base of formation (2) —0.3o 13.6 
(3) —0.26 13.4 
(4) —0.3, 13.7 
(5) —0.27 13.4 | 
(6) —0.67 15.5 
(7) —0.55 15., 
(8) —0.4¢ 14. 
(9) —0.2s 13.5 
(10) —0.6¢ 15.3 
(11) —0.6; 15.5 | 14.4 
Cape Fear River Unknown (1)'—0.54 15.; 
(2) —0.8,5 17.1 | 16.; 
Navesink Mullica Hill, N. J. Base of formation (1) —1.2; 19. 
(2) —1.1s 18.5 
(3) —0.9 17.2 
(4) —0.8; 16.7 
(5) —1.0o 18.5 
(6) —0.9. 17.6 | 18.6 
Hornerstown, N. J. Base of formation (1) —0.8, 17.0 
(2) —0.5s3 15.0 
Crosswicks Creek (3) —0.4, 14. 
(4) —0.7; 16.3 
(5) —0.85 16.6 
(6) —0.8o 16.6 
(7) —0.3, 14.2 
(8) —0.6 15.5 | 16.0 
Cream Ridge Unknown (1) —0.5, 15.4 
(2) —0.7o 16.0 15.7 
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The New Jersey specimens analyzed were 
secured from the Navesink marl at the well- 
known fossil locality in the hillside just south 
of the railroad trestle at Mullica Hill (Weller, 
1907, locality 169), in the cut of a tributary of 


convexa, and Belemnitella americana. The dense 
packing of the pelecypod shells with the ran- 
domly oriented, interwedged Belemnitellas in- 
dicates post mortem concentrations swept to- 
gether in shallow agitated waters. Burial 
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FicurRE 1.—CRETACEOUS SECTIONS FOR LOCALITIES: STUDIED 


Crosswicks Creek near Hornerstown, and from 
the dump piles of the abandoned “marl pit” 
on Edgar Schank’s farm near Cream Ridge 
(Weller, 1907, locality 127). The Mullica Hill 
specimens were derived from a 4-inch shell bed, 
consisting of Exogyra costata, Pectinidae, and 
Belemnitella americana in a marly glauconitic 
sandy matrix about 12 feet above the base of 
the bluff exposure. The preservation of the pele- 
cypods and belemnites is generally poor due 
to leaching and partial limonitic replacement, 
while conspecific forms occurring in shell beds, 
located 10 feet lower and about 3 feet higher 
up in the bluff section are represented only in 
the form of casts and molds. The Hornerstown 
specimens are derived from a shell bed, com- 
posed principally of Exogyra costata, Gryphaea 





grounds foreign to the habitat are thus strongly 
implied for the Belemnitella element. According 
to Spengler and Peterson (1950), the shell beds 
mark the base of the Navesink member of the 
Monmouth formation, previously interpreted 
as a formation and group respectively. A basal 
position for the Mullica Hill and Hornerstown 
specimens in the Navesink is thus indicated, 
while the stratigraphic position of the Cream 
Ridge forms, secured loose from the dump piles, 
cannot be determined. 

All the North American specimens included 
in the present study fall within the range of the 
Exogyra costata zone, as shown by the correla- 
tion of the formations from which they were 
derived (Fig. 1). In comparing the stratigraphic 
positions of the specimens secured from Tennes- 
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see, Mississippi, and South Carolina, and the 
Mullica Hill and Hornerstown localities in New 
Jersey within their respective formations and 
in relation to the formational ranges and their 
correlation, it will be seen that the Coon Creek 
forms, from southwestern Tennessee, occupy a 
lower; those from the Peedee formation of South 
Carolina and from the Navesink localities of 
New Jersey an intermediate, and the Prairie 
Bluff forms from northeastern Mississippi a 
higher position within the zonal range. Their 
stratigraphic relations to the North Carolina 
and Cream Ridge locality of New Jersey cannot 
be determined because of the unknown strati- 
graphic positions of the latter within their re- 
spective formations. To judge from the narrow 
range of the Navesink marl, the New Jersey 
specimens may be anywhere between the top of 
the Coon Creek formation of southwestern Ten- 
nessee and the base of the Prairie Bluff chalk 
in northeastern Mississippi, those from the Pee- 
dee formation of North Carolina a considerably 
wider interval. 

The North American forms serve for an ini- 
tial survey of the special paleotemperature 
range within the time interval occupied by the 
Exogyra costata zone, even though in some cases 
our specimens came from different levels within 
that zone and in other cases the relationships 
within the zone are uncertain. 


DISCUSSION OF THE NoRTH AMERICAN 
TEMPERATURES 


The Exogyra costata zone is correlated with 
the Maestrichian of the European standard sec- 
tion (Stephenson and Reeside, 1938, p. 1637; 
Stephenson et al., 1942, Chart No. 9). Belem- 
nitella mucronatz which occupies in the Upper 
Senonian of northern and central Europe a 
stratigraphic position similar to Belemnitella 
americana in southeastern North America, has 
been considered closely related to the North 
American species (Stephenson, 1923, p. 56). 
Approximate time equivalency of their enclos- 
ing strata is further indicated by other faunal 
and floral similarities. Keeping in mind the as- 
yet existing uncertainties in intercontinental 
correlation, particularly with reference to de- 
tails in time equivalency and length of time 
intervals, it is nevertheless permissible to com- 
pare in general terms the Maestrichian paleo- 
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temperature data from both sides of the Atlan- 
tic that is, of the Gulf Coastal area in North 
America and of southern England and Den- 
mark in northwestern Europe. There is sur- 
prisingly little difference in the temperatures 
recorded by the reliable forms from either side 
of the Atlantic despite the latitudinal range of 
the collecting localities, which extend from 33° 
to about 4i1°N. in North America to between 
about 52° and 56°N. in England and Denmark. 

In surveying the meaning of the temperature 
records only those of the belemnites can be at 
present considered reliable. As the temperatures 
recorded constitute at best mean temperatures 
of skeletal growth periods, they do not neces- 
sarily represent annual mean temperatures of 
the waters in which they lived. This is clearly 
indicated by the temperature range recorded 
by the shell of a recent Chama macerophylla 
from the Bermuda inshore waters, which has 
been analyzed for seasonal effects. Shell ma- 
terial is laid down only after the surface waters 
reach 24°C and extends through the warm part 
of the year, which reaches a maximum at 29°C, 
while the yearly temperature range of the in- 
shore waters covers a range from 18°C to 29°C. 
The mean temperature of 264°C recorded by the 
Chama mazerophylia shell is thus 24°C higher 
than that of the waters. Organisms with a 
narrow temperature range for skeletal growth, 
particularly with a growth period largely off the 
average with reference to the yearly range of 
their waters, would constitute unreliable data 
for climatic records unless forms from identical 
habitat with wider or overlapping skeletal 
growth periods are available for comparison. 
The Inoceramus lineage now under investiga- 
tion may furnish a comparative element for 
correction at a future date. 

Returning to the belemnites, preference for 
cooler waters is strongly implied by their geo- 
graphic distribution. Dacqué (1915, p. 423- 
426), in discussing the paleoclimatology of the 
Cretaceous, pointed out that broad climatic 
zoning is introduced in the beginning of the 
Cretaceous and becomes sharply established 
by Upper Cretaceous time. Differentiation into 
a boreal and Mediterranean-aquaetorial zones 
is fairly well documented, the aquaetorial zone 
being characterized by prominent carbonate 
sediments, reef-forming corals and rudistids, 
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large Foraminifera, and heavy-shelled Nerine- 
idae and actaeonellids. Most convincing is the 
appearance of dwarfed and crippled rudistids 
in the bordering boreal zones both to the north 
and south, as shown by Dacqué’s distribution 
map for the rudistids. Belemnitella was shown 
to be confined to the northern boreal zone of 
the Upper Cretaceous seas. In comparing the 
size ranges of the collections of Belemnitella 
americana secured from the various localities of 
the Gulf Coastal area, an increase in maximum 
sizes attained is clearly noticeable from Missis- 
sippi northeastward to New Jersey. The Missis- 
sippi localities were apparently located in the 
southern border area, as Stephenson and Mon- 
roe record representatives of the rudistid genus 
Durania from the Exogyra cancellata subzone 
of the upper Selma Chalk in which Belemnitella 
americana makes its first appearance. 

The northeastward trend of increasing sizes 
in Belemnitella americana could be correlated 
with cooler water preference, but their mean 
temperature records do not seem to support 
this view. Refinement in the sampling procedure 
contemplated for studies of seasonal variations 
may decide this question. Throughout our dis- 
cussion we have neglected the possibility and 
even probability of some variation of the iso- 
topic composition of the waters, and this may 
explain the nearly identical temperatures of the 
Carolina and Mississippi region as contrasted 
to New Jersey. 

The depth range is another factor yet to be 
considered in the evaluation of the tempera- 
tures recorded by nectonic elements such as 
belemnites. The sediment character and the 
burial assemblages of the American representa- 
tives in the locations examined indicate that, 
though the bottoms were largely below effective 
wave base (except for New Jersey examples, 
which occur in wave concentrated shell beds), 
the abundant and varied mollusk forms point 
to comparatively shallow bottoms, well within 
the euphotic zone. This agrees with Naef’s 
(1922) conclusions ¢f prevalent shallow-water 
habitat for the belemnites in general as implied 
by functional morphology. Thus the tempera- 
ture variation should be characteristic of a 
limited depth range. That Belemnitella lived in 
shallow cooler waters is further indicated by the 
absence of stray individuals in the southern 
Tethys waters. 
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The narrow gradient indicated by the mean 
temperatures of the Belemnitellas across the 
latitudinal width from Mississippi to Denmark 
is surprising. Unless the temperature tolerance 
of skeletal growth of the European and Ameri- 
can species happens to coincide, rather uniform 
climatic conditions are suggested over most of 
the boreal zone during the Maestrichian. Grad- 
ual lowering of the temperatures through Upper 
Cretaceous time toward the Maestrichian are 
indicated in England and agree with the estab- 
lishing of the climatic zones noted by Dacqué. 
Increased mixing of the oceans due to stronger 
current development should have resulted. 
North America and South America were proba- 
bly separated by the circumequatorial Tethys, 
and the waters from the pole area previously 
connected with it via the western interior 
strait became separated at about the critical 
time that the lowering of the late Cretaceous 
temperatures is indicated to have taken place. 
Free mixing of the northern waters due to lack 
of the present existing north-south land barriers 
should have tended to equalize the tempera- 
tures toward the pole. Subtropical climatic 
conditions could have thus been maintained on 
the northern land masses where the dinosaurs 
continued to flourish. 

Cowles (1939; 1940; 1945) has suggested ex- 
tinction of the dinosaurs at the closing interval 
of the Cretaceous due to increased tempera- 
tures. On the other hand continued lowering of 
the temperatures as suggested by the English 
downward trend would have equally caused 
their extinction. We do not wish to emphasize 
this apparent downward trend since our tem- 
peratures depend on single samples in each zone 
and our experience indicates that a variation 
between a number of specimens of the same 
locality of as much as 5°C is possible. The in- 
vestigation of more specimens from these zones 
and seasonal variations of such specimens may 
enable us to reach a definite conclusion. Longi- 
tudinal paleotemperature traverses for the en- 
tire Upper Cretaceous extending into early 
Tertiary, based on a large number of associated 
faunal elements and including benthonic forms 
should furnish a clearer trend of paleoclimato- 
logic conditions as well as aid in quantitatively 
evaluating the temperature niches of the marine 
organisms. 
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ADDENDUM 


It is necessary to consider the possible effect of the changed temperature scale mentioned in 
the Addendum of the preceding paper on the temperatures recorded here. Repetition of the analyses 
using the modified method of preparation and the corrected temperature scale on 10 specimens 
reported in this paper changed their measured temperatures by an average of 1.1°C with a spread 
of —0.6° to +2.2°C. The temperatures remain remarkably uniform though systematically higher 


by about 1°C. 
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CARBONATE-WATER ISOTOPIC TEMPERATURE SCALE 


By SamuEet Epstein, Ratpu BucusBAuM, He1nz LOWENSTAM, AND HAROLD C. UREY 


ABSTRACT 


A relationship between temperature and relative O'* abundance in calcium carbonate in marine shells 
has been determined. If the relative O08 abundance of the water in which the shell grew is known, the tem- 
perature can be determined to an accuracy of +1.0°C. The relative O contents of marine waters increase 
with salinity and, in the case of surface marine waters, with salinity and temperature. 
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INTRODUCTION 


Urey e al. (1951) discuss a method for 
measuring paleotemperatures that depends ona 
variation in the oxygen —18 abundances in cal- 
cium carbonate deposited in the sea. McCrea 
(1949; 1950, p. 849) investigated the isotopic 
composition of calcium carbonate precipitated 
inorganically from water and determined the 
temperature coefficient for the abundance of 
O in calcium carbonate relative to water as a 
function of temperature; he found a value of 
-5.4°C permil increase in the ratio of O* to 
O' in the calcium carbonate. McCrea gave a 
calculation for this quantity including estimates 
of the effect of the translational and hindered 
rotational degrees of freedom in the solid state 
on the distribution functions, and arrived at a 
value of 5.1. Assuming free rotation in the solid 
State Urey (1947, p. 562) calculated a value of 
~—5.7 for the same quantity. This research aims 
to determine this quantity and to set up an 
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equation relating temperature to the difference 
in isotopic composition of samples of carbon 
dioxide prepared from organically precipitated 
calcium carbonate and an arbitrarily selected 
reference sample of carbon dioxide. 

The relationships between temperature and 
O" abundance in calcium carbonate from ma- 
rine shells and that from inorganic calcium car- 
bonate precipitates are similar. But a number 
of special considerations must be taken into 
account in the case of marine shells to find the 
temperature of the sea water in which the ma- 
rine animals lived: (1) the isotopic composition 
of the oxygen of the oceans; (2) the collection 
of marine shell-bearing animals which grow in 
water of constant and known temperatures; and 
(3) the problem of impurities in the carbon 
dioxide prepared from marine shells. 
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RELATIVE O88 CoNTENTS OF MARINE WATERS 


Since the ratios of vapor pressures of H,O"* 
to H,O* are 1.008 and 1.011 at 25°C and 0°C 
respectively (Riesenfeld and Chang, 1936, p. 
127), evaporation and condensation continu- 
ously change the O" contents of bodies of water. 
The magnitude of the variation of O”* in the 
oceans is important to the determination of a 
valid temperature scale and to the general prob- 
lem of the determination of temperature of un- 
known samples, because this thermometer de- 
pends essentially on the magnitude of the 
difference between the O'* contents of the cal- 
cium carbonate and the water in which it was 
precipitated. Consequently a series of O" deter- 
minations were made on samples of sea water, 
most of which were made available to us by H. 
Pettersson, who collected them on the ALBA- 
TROSS expedition (1949, p. 182). The method of 
determining the relative O content of water 
samples involves the equilibration of carbon 
dioxide with the water at 25°C and analysis of 
an aliquot of the gas. The carbon dioxide so 
equilibrated is compared with our standard 
working sample of carbon dioxide prepared from 
the calcium carbonate of belemnites of the Pee- 
dee formation using a high-precision mass spec- 
trometer (McKinney, McCrea, Epstein, Allen, 
and Urey, 1950). The analyses are relative 
only, but are correct within a probable error of 
about 0.2%p of the isotope ratio. The results of 
these analyses are given in Table 1. The unit 
used in this table is defined by the following 


equation: 


3(0%8) = (fe =~ ) 1000, 


standard 


where Reampte aNd Ryandara are the ratios of 
CO"*0" to CO!" for the sample and standard 
reference gas, respectively. 

Apparently, the isotopic composition of oxy- 
gen is not homogeneous throughout the oceans. 
Deviations from the mean isotopic composition 
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TABLE 1.—RELATIVE O CONTENT OF Marine 


























WATERS 
Pepth T | | 
of poin em- 
toctin | S| |S 09 
(Me- | (°C) | 
ters) | | 
30°47'N. 25°25'W. | 10| — | 37.44) -0.0, 
1000 | — | 35.43) -0.5, 
26°25'N. 34°00’W. | 5000} — | 34.9 | -0.2, 
14°27'N. 66°13’W. | 5000 | 4.2 | 34.9 | -0.7, 
14°15’N. 71°02’W. 0 | 27.7 | 36.0] 0.4 
02°52’S. 88°20’W. 20 | 18.8 | 34.4 | -0.3, 
500 | 7.5 | 34.6 | -0.5, 
1000 | 4.5 | 34.6 | -0.5; 
13°25’S. 149°30’W. 0 | 27.7 | 35.8 | -0.y 
1000 | 4.3 | 34.6 | -0.9, 
4000 | 1.4 | 34.7 | -0.6, 
07°45’S. 120°30’E. 10 | 29.5 | 34.5| 0.5, 
10°06’N. 52°45’E. 20 | 29.3 | 35.5 | 0.0, 
500 | 10.7 | 35.0 | -0.2, 
2000 | 2.9 | 34.8 | -0.9, 
35°41'N. 21°50’E. 50 | 16.5 | 38.8] 0.6, 
2000 | 13.6 | 38.7| 0.8 
00°13’S. 18°26’W. 25 | 22.6 | 35.6| 0.0, 
500 | 7.1 | 34.7 | —0.5, 
2500 | 3.0 | 34.9 | —0.3, 
7500 | 1.4 | 34.8 | -0.4 
28°05’N. 60°49’W. 25 | 28.5 | 36.8] 0.3, 
500 | 16.3 | 36.3 | —0.15 
2000 | 3.7 | 35.0 | —0.55 
43°04’N. 19°40'W. | 2000 | 3.6 | 35.0 | -0.6, 
President Channel, | 171 | 9.4 | 31.4 | -1.% 
Puget Sound | 
Middle Channel, | 157 | 7.4 | 33.1 | -0.% 
Puget Sound 
Encenada, Mexico | ~0/| 19 | 33.5 | —0.8 
St. Tomas, Mexico | ~0O | 15 33 | -0.7%1 
Monterey Bay,| ~0/ 11 | 33.8 | —0.8; 
California 
La Jolla, California | ~0| 15 | 34.0 | —0.75 
Bermuda ~0 | 18-29) -- 0.6 





of the oceans occur in the waters close to the 
mouth of a river, as is the case for the President 
Channel location. If water of fully marine com- 
position is considered, then the maximum vari- 
ation of O¥ in ocean water is 1.7%, which 
would represent about 10°C variation in tem- 
perature determination. Clearly, it is necessary 
to determine the O" content of the water in 
which a marine shell grew if this particular shell 
is to be used in correlation of temperature with 
the O* content of organically precipitated shell. 
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On further examination of the data an at- 
tempt was made to establish a correlation be- 
tween the O content and the salinity of the 
water. Such a correlation is possible because 
processes which tend to increase or decrease the 
(0) of the water also tend to increase or decrease 
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Ficure 1.—Pxot oF SALINITY vs. O!8 CONTENT FOR 
OcEAN WATERS 





the salinity. The plot of salinity against relative 
08 content is shown in Figure 1. The points 
scatter badly, but a regular trend is evident, 
and a straight-line relationship is assumed be- 
tween salinity and O" content of the water. The 
relationship is only approximate as was to be 
expected. Extrapolation of the line to zero salin- 
ity gives —13.5%p» as the relative O*8 content 
of fresh water. 

The approximate linear relationship between 
salinity and O'8 abundances suggests that it 
might be more rigorously followed if limited to 
surface waters. Since surface salinity and sur- 
face temperature show some indication of being 
related linearly, we have plotted in Figures 2 
and 3 salinity vs. O%8 content and temperature 
vs. 08 content respectively for the waters of 
Table 1 coming from a depth of 25 meters or 
less. The plots indicate rather definite linear 
relations, but the data are too meager to justify 
any definite conclusion. We shall continue this 
investigation as we secure more water samples, 
particularly from higher latitudes. 

The mean isotopic composition of the oceans 
will be close to that of the water of mean sa- 
linity, and we have selected this as having 
—0.47%p relative to our standard gas and 
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Ficure 2.—P.ot oF SALINITY vs. O!8 CONTENT FOR 
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Ficure 3.—PioT oF TEMPERATURE vs. O!8 Con- 
TENT FOR SURFACE OCEAN WATERS 


corresponding to 34.8%» salinity. This gives 
—13.0%p» as the isotopic composition of fresh 
water relative to that of sea water of 34.8%» 
mean salinity. This difference is larger than that 
expected for a fractionation equal to the ratio 
of vapor pressures, and substantially greater than 
that for any fresh-water samples investigated 
by us. Lake Michigan is representative and 
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analyzes as —6.2%p» relative to our mean sea 
water. 

This large variation of isotopic composition 
of sea water decreases the reliability of tempera- 
ture measurements by this method unless a 
sample of water is analyzed at the same time 
and a correction for the isotopic composition of 
the water is made on the carbon dioxide analy- 
ses. This error, some 2.5 times as large as ex- 
pected, is disappointing, especially when paleo- 
temperatures are considered. The melting ice 
of high latitudes must have a lower O"* content 
than has fresh water of the temperate zone. The 
whole problem of water evaporation from the 
oceans and its precipitation over the land and 
oceans is very complex, and our data are not 
sufficiently extensive to justify a more detailed 
discussion. 

A correlation between salinity and O" con- 
tent is useful in that salinity data which are 
quite readily available for present-day oceans 
can give an approximate correction for the iso- 
topic abundance if water samples are not avail- 
able for mass spectrometric analyses. We hope 
the distribution of O" between phosphate and 
water will give another thermometer with a dif- 
ferent temperature coefficient, thus making 
possible the elimination of the water phase.' 


Isoropic TEMPERATURE SCALE 


Marine Shells Selected and their Growth 
Temperatures 


One problem in establishing a temperature 
scale using marine calcareous-shelled animals 
was to get shells which grew at known tempera- 
tures. To grow certain marine animals under 
controlled temperatures in the laboratory is a 
lengthy and difficult project. Nevertheless, con- 
stant-temperature baths were set up at three 
Pacific Coast marine laboratories,? and these 


1M. Steinberg is developing this possibility in 
our laboratory. 

? Oceanographic Laboratory, University of Wash- 
ington, Friday Harbor, Washington. We are in- 
debted to Dr. Thomas G. Thompson and Dr. Emery 
Swan for assisting in this work, and also to Dr. C. 
Barnes for help in locating suitable areas for collect- 
ing. At the Hopkins Marine Station we are indebted 
to Dr. Rolf Bolin for help in looking after one of the 
experimental tanks. At Scripps Institution of Ocean- 
ograph, La Jolla, California, we are indebted to Dr. 
Carl Hubbs who assisted us in setting up one of the 
experimental tanks and in making collections in 
lower California. 
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were stocked with calcareous organisms with 
notched or drilled shells. It was thought that if 
repair occurred, the new shell, when analyzed, 
would yield reliable data which could be used 
to fix points on the temperature scale. Not all 
of the material is in, and some of the attempts 
were unsuccessful. However one of the experi- 
mental batches of animals has been analyzed 
and the results are presented later. Con- 
sequently, shell-bearing marine animals were 
collected from locations which have a minimum 
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annual variation of temperature. The species 
of animals chosen were attached or sedentary 
forms which live in relatively confined areas 
so that the point at which collection was made 
would be representative of the condition under 
which the animals grew their shells. Animals 
were collected from marine locations which 
had records of maximum constancy and which 
had a minimum of contamination with fresh- 
water drainage. Collections were made in the 
cooler waters of Puget Sound, intermediate 


conditions of Monterey Bay, and the warmer | 


waters of certain areas along the coast of lower 
California in Mexico. In addition, shells from 
the still warmer waters of Bermuda were an- 
alyzed. 

Particular care was taken to estimate care- 
fully the temperature at which the particular 
fragment of shell used for analyses grew. For 
example, the samples analyzed from Santa | 
Tomas, Mexico, included one large red abalone, | 
Haliotis rufescens, and one large black abalone, 
Haliotis cracherodii, and two keyhole limpets, 
Fissurella volcano. The Haliotis were collected 


August 25, 1949, and thus the outer few milli- 2 


meters of the rim of the shell grew within the 
months May to August. Dr. C. L. Hubbs 
measured temperatures at several stations in 
this region during September 1948, May, July, 
and August 1949, and also observed the daily 
variation on August 18, 1949. From these 
measurements we conclude that 15°C is a 
reasonable average temperature. The variations 
in temperature during these months is about 
the same as the daily variation from 14.4°C 
to 17.0°C. Thus a probable error of +1°C 
seems a reasonable estimate. A complete cross 
section of the shell of the keyhole limpets was 


taken so that it was necessary to estimate the | 


average annual temperature. The temperature 
in this location goes as low as 11°C and as high 
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TABLE 2.—DESCRIPTION OF SHELL SAMPLES USED FOR DETERMINATION OF TEMPERATURE SCALE 
— Specimen a 5 oa . 
Location Part of specimen used Temp. | a 3 
No. | Name °C “a 
Middle Channel, Puget 1 Mussel Volsella modiolus | Cross section of complete | 7.4 —0.9; 
Sound 48°22'30"N., half shell 
122°58'45”W., 130-157 
meters, July 22, 1949. 2 | Brachiopod Terebratulina | Complete half shell 
transversa 
President Channel, Puget 3 | Mussel Volsella modiolus | Cross section of complete | 9.4 | --1.9% 
Sound 48°41'50’N., half shell 
122°59’30”W., 171 me- 
ters July 16, 1949 4 | Brachiopod Terebratulina | Cross section of complete 
transversa half shell 
Hopkins Marine Station, 5 | Red abalone 4Haliotis | At least a year’s growth | 12.6 | —0.85 
Pacific Grove, Calif. rufescens. At low tide of shell 
level May 14, 1949 
6 | Giant keyhole limpet | Cross section of complete | 11.5 
Megathura crenulata 50 shell 
m. off-shore, 3 m. 
Aug. 15, 1949 
Santa Tomas, Lower Cali- 7 | Red abalone dHaliotis | Last 4 cm of large shell | 15 | —0.77 
fornia, Mexico Aug. 27, rufescens 
1949 3-6 meters 
8 | Black abalone Haliotis | Last $ cm of large shell | 15 
cracherodit 
9 | Limpets Fissurella volcano | Cross section of 2 com- | 14 
10 plete shells 
Between Encenada and | 11 | Red abalone MHaliotis | Edge of shell 20 —0.8 
Morro Points, Lower rufescens 
California, Mexico. 3 
meters, Aug. 25, 1949 
Bermuda, Tank in labo- | 12 | Snail (conch) Strombus | Regenerated shell 31 0.64 
ratory gigas 
13 | Bivalve, Pinna sp. Regenerated shell (part) | 29 














| 








as 17°C. The average temperature of this 
sample of the shell was thus estimated to be 
14°C within an error of less than 3°C and 
perhaps as close as +1°C. Similar considera- 
tions of temperature data made available to us 
led to the assignment of other temperatures, 
and similar uncertainties exist in the estimate. 

Table 2 describes the location at which the 
marine animals were found, the estimated tem- 
perature of the fragment of shell used for analy- 


ses, and the relative O" of the water collected 
with the sample. 

The carbon dioxide was prepared from the 
marine shells by reacting the calcium carbonate 
with approximately 100 per cent ortho phos- 
phoric acid in an evacuated vessel. The carbon 
dioxide was trapped in a liquid nitrogen cooled 
trap and pumped free of the more volatile 
impurities, then passed through traps covered 
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by a solid CO,-acetone bath to freeze out such 
less volatile impurities as water. 

The samples of the shells described in Table 
2 were then cleaned of obvious impurities and 


isotopes of the calcium carbonate during their 
deposition in their complex shell structure, 
But not enough information has been obtained 
to justify this conclusion regarding these groups 


TaBLeE 3.—Isotopic ANALYSES OF ORGANIC CALCIUM CARBONATE SAMPLES 


Each isotopic analysis is a mean of at least four sample preparations. The measurements were all within 
the accuracy of +0.2%p». 

















Numberof spec) Qhverved |, Apatzis (No, | Analysis (With | water correction | "(No heat | With heat 
reatment) treatment) 

1 7.4 1.4; 0.36 0.4 | 1.8, | 0.85 
2 7.4 1.3, 0.1; 0.4, = | 1.75 0.5; 
3 9.4 0.5, —1.2; 1.4, | 2.0; 0.25 
4 9.4 —0.2 —1.14 1.4, | 1.26 | 0.3; 
5 12.6 0.50 —0.51 0.33 | 0.8 | 0.1 
6 11.5 1.06 —0.3; 0.35 | 1.4, | 0.0; 
7 15 0.% —0.9; 0.30 | 1.22 ~0.6; 
8 15 0.95 —0.8) 0.30 | 1.2, | 0.5 
9 | 14 0.60 | 0.45 0.30 0.% | -O.% 
10 | 14 0.46 —0.5o 0.30 | 0.76 —0.25 
11 =6|~— (20 0.2; 1.8 0.3, | 0.6. | 1.45 
12 | 31 —0.5s 2.2; —1.1, | 1.6, | 3.3, 
13 | 29 — —1.75 —1.1, | — |} 2.86 














the carbon dioxide prepared. Table 3 gives the of invertebrates. In any case, the shells here 





results of analyses when this method of prep- 
aration of carbon dioxide is used. The speci- 
men numbers correspond to those in Table 2. 
The results of the sixth column are plotted in 
Figure 4. 

The data show a correlation between tem- 
perature at which the marine animals grew 
and the relative abundance of O" in the carbon 
dioxide prepared from the shells. But the devia- 
tions from a linear or nearly linear relationship 
are much larger than the expected experimental 
error. The difference in analyses for shells col- 
lected simultaneously from the same location 
eliminated the possibility that the poor correla- 
tion between temperature and relative O"8 con- 
tent was due to incorrect estimation of the 
temperature at which the animals grew. 

The possibility that the animals did not lay 
down their shells in equilibrium with the sur- 
rounding water seemed to be unlikely because 
the relative O contents of calcium carbonate 
showed a definite correlation with temperature. 
Some preliminary results obtained during the 
investigation indicate that perhaps the echi- 
noderms and corals fractionate the oxygen 


used in the determination of the temperature 
scale were those which are deposited external 
to the animals’ soft parts. 

Comparison of results of McCrea with those 
of Table 3 indicated that the possible errors in 
the present work were of the order of magnitude 
of 1% 9 of the mass 46 to mass 44 ratio. Since 
impurities of mass 46 present to the extent 
of three parts per million would introduce an 
error equivalent to 1%p in the ratio, the usual 
chemical methods for analysis of impurities 
would be impractical in this case. Therefore, 
the criterion for valid analyses became the con- 
sistent correlation between temperature at 
which the shells were laid down and the meas- 
ured change in relative abundance of mass 46 
in the carbon dioxide. 


Extraction of Carbon Dioxide, Problems of 
Impurities, and Results 


The method of collecting the ion currents 
in the mass spectrometer gave measurements 
of the abundance ratio of mass 46 to the sum 
of masses 43, 44, 45, and 47. Consequently, 
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although the ratio of mass 46 to mass 44 would 
be relatively insensitive to small concentrations 
of impurities of masses 43, 45, and 47, these 
had to be considered. Purification of samples 





30 


25r 


°C 


SF SAMPLES _ WITHOUT 
HEAT TREATMENT 








0 N ! ! 1 ! ! 

-30  -20 =10 o «10 20 
&(0") %o 

Ficure 4.—Piot or O!8 CONCENTRATION vs. TEM- 


PERATURE: CALCIUM CARBONATE SAMPLES 
WitHout Heat TREATMENT 





could then either decrease or increase the 
value of the 46 to 44 ratio depending on whether 
purification eliminated impurities of mass 46 
or of masses 43, 44, 45, and 47. 

Many molecules or radicals formed in the 
ionizing chamber of the mass spectrometer, 
both organic and inorganic can contribute to the 
ions of mass 43 to 47. Some of these can be 
C,H,OH, CS, CH;COOH, BCI, NO,, N,O, and 
propane. The first attempts were focussed on 
the possible impurities in the carbon dioxide 
gas. Several carbon dioxide samples from differ- 
ent marine shells were mixed with small volumes 
of bromine gas and irradiated by strong light, 
with the hope that bromination would shift 
the mass range of the impurities beyond the 
43 to 47 range. Some samples of carbon dioxide 
were mixed with Van Slyke’s (1940, p. 509) 
oxidizing solution or passed over hot CuO. 
Organic impurities would then be converted 
to carbon dioxide, and if present in trace 
amounts should not affect the 46 to 44 ratio 
of the carbon dioxide sample as a whole. Sam- 
ples were also mixed with traces of ammonia 
and water to hydrolize compounds like boron 
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trichloride, and passed over hot copper to 
reduce possible traces of nitrogen oxides. 

In all these cases no appreciable effect on the 
isotopic analyses was observed. 

Samples of carbon dioxide from the shells 
were also mixed with fluorine. Only where 
there was evidence of incomplete separation 
of the carbon dioxide and fluorine and sub- 
sequent contact with glass was there a change 
in O"¥ analysis. These changes were erratic for 
any one source of gas samples and varied from 
—0.3% to —200%p» of the 46 to 44 ratio. This 
was obviously due to SiF* of mass 47 adding 
to ion beam 44 and thus showing a decrease in 
the ratio of 46 to 44. 

A more successful approach to the purifica- 
tion problem involved treatment of the pow- 
dered calcium carbonate before preparation of 
the carbon dioxide. The powder was heated at 
400°C for 60 minutes with a slow stream of 
helium passing over the powder. The helium 
was first purified by passing it over copper oxide 
at 750°C and then through a liquid nitrogen 
cooled trap. With this treatment the relative 
O' contents of the calcium carbonate from 
shells was decreased in most cases, and in a 
reproducible manner. 

Experiments to determine the length of time 
necessary for heating the calcium carbonate 
powder to attain constancy of O¥ analysis with 
further heating showed that heating at 400°C 
for 15-30 minutes was necessary to “purify” 
the samples tried. A safety factor of two was 
adopted, and the samples were heated for an 
hour. Since the equilibrium pressure of carbon 
dioxide over calcium carbonate is 10-* atmos- 
pheres at 400°C, it was necessary to test this 
treatment for possible loss of carbon dioxide 
and fractionation of the oxygen isotopes. Two 
samples of calcium carbonate, one normal in 
O* and the other enriched by 100%p, were 
precipitated inorganically (using sodium car- 
bonate and calcium chloride). Samples of each 
were analyzed without the heat treatment and 
with heat treatment. The analyses of each 
variety of calcium carbonate checked, showing 
that heat treatment did not change the analyses 
of such samples. Clearly, it was the organically 
precipitated calcium carbonate alone which 
was affected by the heat treatment. To shorten 
the time for the experiments the powdered 
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samples were heated for 15 minutes at 475°C 
with the attainment of the identical results 
observed for the longer heating times. We do 
not understand the mechanism of purification 
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of the marine shell samples by the heating proc- 
ess. In most cases it is obvious that charring 
of organic matter occurs. 

Calcium carbonate samples from specimens 
listed in Table 2 were heated before the carbon 
dioxide was prepared, and then analyzed. The 
results are given in Table 3, last column, and 
plotted in Figure 5. 

Comparisons of Figures 4 and 5 indicate 
that the heat treatment of the organically 
precipitated calcium carbonate lowered the ap- 
parent O* contents of the carbon dioxide pre- 
pared from all these samples. Apparently the 
marine shells contained variable amounts of 
some impurity, causing the inconsistency in 
the relationship between temperature and O¥ 
content when unheated samples were used. 

Least-squares calculations of the data for 
the heat-treated samples in Table 3 gave a 
straight-line relationship 


# = 11.88 — 5,915(0%), 


where ¢ equals the temperature in °C, 5(O") is 
the per mil difference between the ratios of 


masses 46 and 44 for the sample and standard 
gas. The constant 11.88 is characteristic of our 
standard gas, while the constant 5.91 should 
be independent of this gas. 


DISCUSSION 


Our constant, 5.91, for the slope of the curve 
differs from that secured by McCrea (5.4). The 
deviation is in the direction expected for in- 
complete equilibrium between solid calcium 
carbonate and dissolved carbonate ion as found 
by him. We do not believe this is the cause of 
the disagreement since a considerable variety 
of organisms all give analyses lying so closely 
on our curve.? The measurements of these small 
quantities are very difficult and, though we do 
not understand the reasons for the disagree- 
ment, it is only 10 per cent in the slope of the 
curve. Rather than regarding this as disagree- 
ment, we believe it shows that these living or- 
ganisms approximate equilibrium conditions 
rather closely. Analyses at present are more 
reproducible than they were when McCrea 
completed his work, and this may have signifi- 
cance in regard to the discrepancy. 

The possible effects of the difference in crys- 
tal structure of calcite and aragonite have not 
been explored. It is hardly to be expected that 
any large effects will be present, but, as our 
techniques improve, they should be investi- 
gated. Not only may the isotopic composition 
of aragonite and calcite crystals grown under 
similar conditions differ, but our heating pro- 
cedure may convert aragonite to calcite and 
thus modify fractionation effects in the liber- 
ation of carbon dioxide by phosphoric acid. 

The presence of a physiological effect in the 
case of certain groups of animals such as the 
echinoderms and corals, and plants such as a 
coralline algae, has seemed probable at some 
points in our work, but we do not feel it is 
certainly established. These specimens are diffi- 
cult to investigate because of the large amounts 
of organic matter present, but, since we have 
established a rather secure basis for our work 
on the mollusks and brachiopods, it now seems 
feasible to investigate these groups. The echino- 


3 As this is written an analysis of a sample of 
regenerated shell from Haliotis cracherodit grown in 
a bath at a controlled temperature of 18.9 + 0.1°C. 
shows agreement with our curve within 0.3°C. 











ar 
Wi 
liv 


thi 


or 
ple 


of ; 


Mc 








ndard 
of our 
hould 


curve 


). The 


or in- 
cium 
found 
use of 
ariety 
losely 
small 
we do 
agree- 
of the 
agree- 
ng or- 
litions 
more 
cCrea 
ignifi- 


1 Crys- 
ve not 
d that 
aS OUT 
iesti- 
sition 
under 
g pro- 
e and 
liber- 
‘id. 

in the 
as the 
h as a 
; some 
] it is 
e diffi- 
nounts 
e have 
r work 
seems 
chino- 


nple « of 
‘own in 
. 0.1°C. 
e. 





derms and corals are not particularly useful 
in relation to the paleotemperature scale be- 
cause of extensive infiltration of the fossils. The 
Foraminifera should be investigated for, even 
if not useful over long periods of time, they 
may be of interest for the investigation of more 
recent temperatures. 

The variation of O¥ content in marine and 
fresh waters, together with precise methods of 
measurement of the small quantities involved, 
makes possible an investigation of meteorologi- 
cal effects of considerable interest. The fraction- 
ation of the oxygen isotopes in evaporation and 
precipitation of water should give some infor- 
mation in regard to these processes as they 
occur in nature. We intend to continue our 
studies on the isotopic composition of water 
since the problem bears so directly on the 
measurements of paleotemperatures. 


CoNCLUSIONS 


The variation in isotopic composition of 
marine waters is appreciable in terms of the 
isotopic temperature scale. When corrections 
are made for the isotopic composition of sea 
water in which mollusks and brachiopods have 
lived, a linear relationship between temperature 
and the O¥ concentration exists. The slope of 
the curve agrees approximately with calcu- 
lated values and with that secured when in- 
organically precipitated calcium carbonate sam- 
ples crystallized at known temperatures are 
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used. Organically precipitated calcium carbon- 
ate gives reproducible analyses if it is first 
heated to 475°C for 15 minutes and then 
reacted with phosphoric acid to prepare carbon 
dioxide, the gas used in the mass spectrometer. 
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ADDENDUM 


We have delayed publication of this paper while it was in galley proof because we discovered 
disturbing variations in our temperature determination. We traced this to the use of heated copper 
oxide to remove impurities from helium gas passed over our samples during the heating process. 
Small amounts of oxygen from the copper oxide were present in the helium, and the carbon dioxide 
formed apparently exchanged in some cases with the calcium carbonate of the samples. Removal 
of all oxygen by means of a charcoal trap at liquid nitrogen temperature corrected these difficulties. 
Our tentative corrected temperature equation becomes t = 14.8-5.415. The slope now agrees with 
McCrea’s value. We will publish details of this temperature scale shortly, 





